foe] 
be 
© 
ro) 
™ 


PROCEEDINGS OF 
E AMERICAN SOCIETY 


F CIVIL ENGINEERS» 


MARCH 1981 

= 


3 $ 
< 
< 


d 


voL.107 NO.GT3. MARCH 1981 


RECA. 


qj 


= 


OF CIVIL ENGINEERS | 


American Society 


Ay 


AMERICAN ai 
CIVIL EN ENGINEERS 


 Irvan F. Mendenhall = Martin G. Abeg Richard W. 
Past President John A. Focht Paul R. Munger 
Joseph S. William R. Neumen 
President Elect GEOTECHNICAL ENGINEERING DIVISION 
James R Sims Executive Committee 
Vice Presidents = Robert Schuster, Chairman 
D. Bay Lyman R. Gillis Wilt T. Chairman = 
Francis Albert A. Grant Harvey E. 
Floyd A. Bishop —_Willliam R. Neuman 
L. Geary Byrd Leonard S. Oberman Publications Committee 


Warren J. Bake Ladd 
Don C. Banks Poul V. lade 
Richard W. Karn William H. Taylor James M. J. Langfelder 
Leon D. Luck Stafford E. Thornton Chandra S. Brahama Felipe A. Len- Rios 
Arthur R. McDaniel Robert E. Whiteside John T. Christian Gholamreza Mesri 
Richard S. Woodruff G. W. Clough J. Murphy 


Sergio Gonzalez-Karg 


Feeser D. Parkhurst William F. Marcuson Ill, Chairman 
E. Humphrey, Jr. 


John A. Focht, Jr. R. 


Herbert H. Einstein Thom L. Neff 
G. Franklin Michael W. O'Neill 

Julie E. Gibouleeu, to the H. Gray Jean H. Prevost 


Director 


q Bobby Hardin “og Adel Saada 
Louis L. Meier, Washington Counsel/ Assistant 


Cornelius J. Higgins Surendra K. Saxena 
Secretary William H. Robert L. Schiffman 
‘ William H. Wisely, Executive Director Pen, Robert D. Holtz _ Woodland G. Schockley 


Michael N. Salgo, Treasurer M.Idriss = Marshall L. Silver 
Elmer B. Isaak, Assistant L.Hoirwin Glen S. Tarbox 

 STAFFDIRECTORS ‘4 


> at Reuben H. Karol > G. R. Thiers 
34 D. D. Treadwell 
Donald A. Buzzell, Managing Director for D. Kovacs R. Uiirich 
_ Education and Professional Affairs __ ‘3 Leland M. Kraft | J. Lawrence Von Thun t 
Robert A. Crist, Jr., Managing Director Raymond J. Krizek R.N. Yong 
Publications and Technical Affairs E. T. Selig, Exec. Comm. Contact 


Alexander Korwek, Managing Director for 
Finance and Administrative Services PUBLICATION SERVICES DEPARTMENT 


Alexandra Bellow, Director, Human Resources Po David Dresia, Director, Publications — . 
David Dresia, Director, Publications Production and Marketing 
Production and Marketing chni Prof Pp 
Barker D. Herr, Director, Membership. ond 
Richard A. Jeffers, Controller Richard R. Torrens, Manager 
Carl E. Nelson, Director, Field Services = Joseph P. Cerami, Chief Copy acai 
Don P. Reynolds, Director, Policy, Planning Linda Ellington, Copy Editor ‘a 
and Public Affairs i _ Thea C. Feldman, Copy Editor 
Bruce Rickerson, Director, Le islative S Services Meryl Mandle, Copy Editor 
James M. Shea, Director, Public — “Joshua Spieler, Copy Editor 
Communications Shiela Menaker, Production Co-esdineter 
Albert W. Turchick, Director, Technical C. Scheblein, Draftsman 
Services 
George K. Wadlin, Director, Education 
aR. Lawrence Whipple, Director, Engineering | 
_ Management ervices 


L 


PERMISSION TO PHOTOCOPY JOURNAL PAPERS 


ont Permission to photocopy for personal or internal reference beyond the limits 4 4 


Lo im Sections 107 and 108 of the U.S. Copyright Law is granted by the American 


Society of Civil Engineers for libraries and other users registered with the Copyright d : 
_ Clearance Center, 21 Congress Street, Salem, Mass. 01970, provided the appropriate a i 
_ fee is paid to the CCC for all articles bearing the CCC code. Requests for special | 7 
_ permission or bulk copying should be addressed to the Manager of Technical and — > a 
Professional Publications, American Society of Civil Engineers. 


CONTENTS 


vi Vibration Tests of Full-Scale Earth Dam 


erage Study of Dynamic Response of Earth Dem 


Samuel P. Clemence and Albert O. 
_ Geotechnical Considerations for ( Construction i in Saudi ‘Arabia a 


by T. William Lambe, W. Allen Marr, and — es Silva 


office is at 345 East 47th Street, New York, N.Y. 10017. Address all ASCE correspondence 

to the Editorial and General Offices at 345 East 47th Street, New York, N.Y. 10017. ao 

Allow six weeks for change of address to become effective. Subscription price to members 

is $16.00. Nonmember subscriptions available; prices obtainable on request. Second-class | 

~postage paid at New York, N.Y. and at additional mailing offices. GT. | eer 

_ The Society is not responsible for and statement made or opinion expressed in its F 

crea Win cha of of tie 


= 

| 
| a 

| 

| 


DISCUSSION 
D 


sor 160 


Strength and Deformability of Highly Fractured Rock,” 
M. Raphael and Richard E. Goodman (Nov., 1979). 


Shear of Soils," by Edward G. ‘Prater 


_ by Yoshiharu | Moriwaki and Maurice S. Power . 367 ye 


Self of San Bay Mud," 


Wayne Clough and Gordon M. (Jan. » 1980). 


4 
J. C. P. Dalton 


Comparison Used in Tank- Failure Study," * by Roy A. Bell on 


_ by Robert E. Bigh 


Loading on ‘Single Piles in Sand," Sin- Fatt Chan 
T 


homas H. Hanna (Feb., > 


by N. K. Jain and Gopal Ranjan . 


. 


_ “Discussion period closed for this paper. Any other discussion received during this 


closure 


" Microgravity Surveys for Evaluation of Elevation Changes Due to 
| Reservoir Impoundment Waal 
65 
q 
d 


16 6096 VIBRATION TESTS OF FULL-SCALE EARTH DAM 


WORDS: Dams (earth); Dynemic characteristics; Denies stability; ow 

Dynamic tests; Earth fills; Excitation; Experimental data; Measuring 
instruments; Natural frequency; Pressure waves; Vibration 
Vibration response; Vibration tests = 

_ABSTRACT: Extensive full-scale tests, ambient (unique in pth 
dam research) and forced vibration as well as popper tests were conducted on the a 
Santa Felicia earth dam in Southern California. For the forced vibration tests, the dam | a a 
was excited into resonance in both the upstream-downstream and longitudinal — : 
directions. During the ambient test, the naturally occurring vibrations of the dam were 
‘measured. The popper test involved the use of pressure waves originating from a 
controlled, submerged release of gas. During the tests, three-dimensional measurements _ 
of the motions of approximately 25 stations along the crest and seven stations on the o¢ 
downstream slope were recorded and then analyzed in both time and frequency - 
domains. Modes of vibrations and associated natural frequencies as well as damping 


REFERENCE: “Abdel- Ghattar, Mined M., and Scott, Ronald F “Vibration Toots of 
‘Full- Scale Earth Dam,” Journal of the Geotechnical Engineering Division, ASCE, Vol. | 


_ 107, No. GT3, Proc. Paper 16096, March, 1981, pp. 241-269 pe i eis 

16097 DYNAMIC RESPONSE OF EARTH DAM paleo 


KEY WORDS: Damping; Dams (earth); Dynamic characteristics; Dynamic 
‘response; Dynamic tests; Earth fills; Earthquakes; Experimental data; 
Mathematical models; Natural frequency; Shear strain; 


ABSTRACT: Full- dynamic characteristics an earth 
fron determined from imposed dynamic forces are compared with those estimated 


from earthquake responses. The measured frequencies and modes of upstream-— 
downstream vibrations are compared with those predicted by an existing two- 
dimensional shear-beam model. Shear strains and shear moduli were estimated from _ 
_ the full-scale tests and were plotted along with those estimated from the earthquake — 
response analysis. The same procedure was followed for damping versus shear strain at 
different levels. Then the combined results were interpolated and extrapolated. The 

_ dynamic properties of the dam’s constituent materials estimated from low-strain full- bei 

_ scale tests are consistent with those determined from relatively large strains induced by © 
the two earthquakes. Also the tests revealed substantial change in the dynamic 
properties of the dam; the behavior is typical of a yielding dynamic system. ect ica 
_ REFERENCE: Abdel-Ghaffar, Ahmed M., and Scott, Ronald F., “Comparative Study _ 
of Dynamic Response of Earth Dam,” Journal of the Geotechnical Engineering 


- Division, ASCE, Vol. 107, No. GT3, Proc. Payer 16097, March, 1981, PP. es 286 Bist 


16100 PROBABILITIC EVALUATION OF LOADS Wen JAD 


Offshore structures; Pore pressure; rt 


i _ ABSTRACT: The loads that effect geotechnical designs are reviewed, characterizing 
the uncertainties in each load component. Probabilistic models for evaluating cach 

' | component are examined. Whenever data are available, the level of variabilities and 
uncertainties associated with the effect of each component on the foundation is : 

» assessed. Satisfactory performance of a geotechnical system relies, in part, on a food 7 

4 understanding of the characteristics of the loads and environments to which the an 

| is subjected. Uncertainties do exist in these factors. For dynamic loads due to 


a _ earthquake excitation and wave action, both their magnitude and frequencies ae 


| occurrences over the expected duration of the system could not be known. The long- 


term stability of soil would be dependent on stochastic fluctuation of 


pore _ pressure due to seasonal variation and changes in the 


| REFERENCE: Tang, _ Wilson H., “Probabilistic Evaluation of Loads,” Journal of the — 
'! Geotechnical Engineering Division, ASCE, Vol. 107, No. G GT3 Proc. Paper 16100, i 
March, 1981, PP. 287- 
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16106 DESIGN CONSIDERATIONS FOR COLLAPSIBLE SOILS 


KEY WORDS: Bearing capacity; Collapsible soils; Design; Design practices; a 

ABSTRACT: Metastable or soils are as any unsaturated soil that 

_ goes through a radical rearrangement of particles and great loss of volume upon 
‘worng, with or without additional loading. The existing types of collapsible soils and wa 
_ theories developed to explain collapse phenomenon are reviewed, as are the methods 
of recognizing collapsible soils in the field and laboratory. Prediction techniques for the — 


| moat of collapse are presented and summarized; current foundation design 


+ os _methods on collapsible soils are explored. Finally, present methods for the stabilization 
| of collapsible soils and possible future methods of treatment are reviewed. ca (<n a. 
REFERENCE: Clemence, Samuel P., and Finbarr, Albert O., “Design Considerations 
for Collapsible Soils,” Journal of the Geotechnical Engineering Division, ASCE, Vol. ie 
107, No. GT3, Proc. Paper 16106, March, 1981, pp. 305-317 vine 
din¥ BHA leat cts) Fo nel 
GEOTECHNICAL CONSIDERATIONS IN SAUDI ARABIA 
| KEY WORDS: Construction materials; Construction planning; Construction Sv 
procedure; Corrosion; Foundations; Geological surveys; Geotechnical 
engineering; Rocks; Saudi Arabia; Soil mechanics 
design and construction are intimately related to the environment at particular 
locations. Coastal areas are generally characterized by particularly soft and loose 
_ corrosive soils. Limestones in the interior are characterized by the potential for the 2 
presence of dry cavities and solution channels. Igneous and metamorphic rocks of good 1 
Bering quality characterize the high mountains overlooking the west coast; 
superimposed areas of sand dunes and salt flats could present special sig and 
_ problems. The availability of construction materials depends on bedrock geology and 
_ varies from relatively good in the western areas to poor in the eastern region. Offshore — 
conditions in western coastal areas are characterized by varied subsurface conditions at 
relatively short distances, affecting the design of offshore facilities. =~ + 
REFERENCE: Oweis, Issa, and Bowman, John, “Geotechnical Considerations for 
- Construction in Saudi Arabia,” Journal of the Geotechnical Engineering Division, 
ASCE, Vol. 107, No. GT3, Proc. Paper March, 1981, pp. 319-338 


16107 GEOTECHNICAL ASPECTS OF FACILITY SAFETY 
WORDS: Earth dams; Earth dam seepage; Facilities; Geotechnical 
engineering; Pore water pressure; Risk prediction (accidents); Safety; Safety — abaot 
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outlining, discussing, and illustrating various components of a comprehensive amy 
program. When dealing with important complex geotechnical facilities, experience | 
clearly shows that the geotechnical profession cannot feasibly create a facility which 
_ completely fulfills all performance criteria for the entire life of the facility. Rathez, = 


ate An. approach i is presented to to ensure the safety of facility 


_ only logical approach is to engineer a facility for its design life, i.e., create a reasonable 
design, construct the facility according to the design, and then execute an appropriate 

“ongoing safety _ program. The comprehensive Safety program consists of nine 

experience f These components are described, and Several are illus illustrated, using 

REFERENCE: Lambe, ‘William, Marr, W. Allen, and Silva, Franciso, “Safety of 

_ Constructed Facility: Geotechnical Aspects,” Journal of the Geotechnical Engineering — 

_ Division, ASCE, Vol. 107, No. GT3, — — 16107, March, 1981, pp. 339-352 


cee 


| 
€ 
« 


JOURNAL OF 
ENGINEERING DIVISION 


IBRATION Tests OF F ULL-SCALE EARTH Damo 


y Ahmed M. Abdel- 1-Ghaffar,’ A. M. ASCE and Ronald F. Scott, M. ASCE 


Full- scale experimentally determined dynamic characteristics are required both 


to improve the techniques by which such properties ¢ are calculated for purposes at 
of analysis and design, and to ‘permit interpretation of the measured response Me 
of earth dams in the event of a strong earthquake. In addition, as information _ 
on the dynamic characteristics of dams and their earthquake response accumu- - 
ie lates, it should be possible to develop better procedures for determining the 
appropriate levels of earthquake loading for design. 
‘This paper concerns full-scale dynamic tests, including (unique in 
— earth dam research) : and forced vibration a as well as popper tests (12), that 
_ were carried out on Santa Felicia Dam in southern California subsequent to q 
study of the dam’ s to two earthquakes (3 and 4). A 


to the t two o earthquakes and 4). 

_ Some notable improvements in these tests, added since most previous earth — 
: _ dam tests (6, 10, 11, and 13) were made, are: (1) Ambient vibration measurements; _ a 

= (2) the popper test (previously proved very efficient in testing arch and gravity 
_ dams); (3) longitudinal shaking of the dam, and determination of the dynamic — 
properties in the longitudinal direction not only. from the forced vibration tests _ 
but also from the ambient tests (the data obtained have great potential for al 
future research); (4) measurement of the modal three-dimensional Tesponse 
of the dam; and (5) availability of usable earthquake response records to com- 
e pare with the full-scale low-level induced response measurements. A detailed de- | 


Prof. of Civ. . Engrg., California Inst. of Tech., Pasadena, Calif. = Diet sod ke 
Note.—Discussion open until August 1, 1981. Separate discussions should be e submitted — 
J for the individual papers in this symposium. To extend the closing date one month, 
Sed a written request must be filed with the Manager of Technical and Professional Publications, — 
_ ASCE. Manuscript was submitted for review for possible publication on April 17, 1980. ; 
- This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings ' ' 
7 ¢ the American Society of Civil Engineers, ©ASCE, Vol. 107, No. GT3, March, 1981. 


4 Asst. Prof., , Civ. Engrg. Dept., Princeton Univ., Princeton, N.J. 08544. 
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a companion paper (2) compares the dynamic properties determined herein from 
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Santa 1 and 2) is a rolled- fill built 
y in 1955. It is 273 ft (83 m) high above its lowest foundation and 200 ft (61 


m) above the original stream bed; it is 450 ft (137 m) long across the valley 
at the base and 1,275 ft (389 m) long at the crest. The dam has a central 
. impervious core and pervious shells upstream and downstream. The core and 4 
_ shell materials are basically alluvial, consisting of clays, sands, gravel, and — 
boulders. A more complete description of the dam is given in —_ | and be 
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FIG. 1.—Cross Section of Santa Felicia Dam Showing be aaa Stations or on 
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other references (8, 9, 10, and 11). It basically consists of: (1) Motion ‘sensing 
including conditioners (SC-1 also by Inc.), and 
i recorder (7418A by Hewlett-Packard), and instrumentation tape recorders (3460A 
by Hewlett-Packard); (3) force generating systems or shaking machines capable 
Be producing a force up to 10,000 Ibs (44,500 N) [designed at California Institute xu 
of Technology (Caltech) and manufactured by Kinemetrics, Inc.]; and (4) data _ 
processing instruments including an electronic analog- digital converter (DDS-1103 — 
Dynamics Corp. 


Ser 


50 ft (15 m) (Fig. 1); these stations were designated by El-E12 for the eastern 
half and W1-W12 for the western half, with the center station designated 0. 
There were two arrays of measurement stations on the | downstream Slope (Figs. : 

~ and 2), each consisting of seven stations [at intervals of ” approx 50 ft (1S 
-m) along the slope]; one array was along the midpoint cross section, for 
measurement of symmetric modes, while the other was even with station E4, > 
_ for measurement of antisymmetric modes. One vibration generator was mounted — 


at station E2, and the other was placed at station W2 [the two stations were — 
200 ft (61 m) apart, Fig. 2}. The aforementioned recording instruments were = 

in the region of the crest between the shakers (close to 


"Experimental Procedure and Results. —The dam was excited into resonance 
in various modes in the upstream- -downstream (U-D) direction by the coupled 
_ pair of mechanical shakers. Symmetric and antisymmetric vibrations were 
obtained by synchronizing the two shakers to run in-phase and 180° out- of-phase, 
4 Bh gyno bee the aid of a precise and stable control system. The measurements | 
taken (in the frequency range 1. 0 Hz-6.0 Hz) on the dam included, as a first . 
step, the determination of the natural frequencies (via frequency sweeps), — 
4 preliminary modal damping ratios and mode shapes along the crest and the } 
downstream face. Eight seismometers, located at selected measurement stations “a 
to cover the entire dam crest and downstream face, were used to pick up c 
the steady-state vibrations of the dam; they were oriented in the direction of 
the exciting harmonic force. In addition to the measurement of 17 modes in 
‘the U-D direction, eight modes “of vibration of the dam’ s three-dimensional 
(3-D) response were measured using six seismometers (three were placed at 
the reference station W1, where they remained throughout the tests while the 
other three were atten from station to station). ii 3 ‘summarizes the instrumen- — 


A 
so _ The essential features of the instrumentation are substantially the same as _ 
| 
i 
i 
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} FIG. 4.—Symmetric Shaking: (a) Hy ayy Curves; (b) Normalized Response Curves | 


_ Resonance Curves and Modes of eaabiee —Figs. 4, 5, and 6 show some 
= of the 1 resonance curves  dcteemined for three selected stations. i 
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q used in the shakers; the force level (or the weight combinations) for each segment 
" of the curves is indicated on the figures. In plotting these curves, the magnitude 


4 in these curves resulted when changes were made in the weight combinations a 


_ of the double-amplitude of the response was normalized [divided by the largest a 


- response obtained for each set of weights, as shown in Figs. 4(a), 5, and 6]. _ 
_ Another alternative way of plotting the results is shown in Fig. 4(b), where = 


4 
pattern is the same as that assumed in the shear-beam theories for a triangular __ 
2 elastic wedge in a rectangular canyon (5 and 11). Peaks designated by R are _ 

believed to be either rocking modes (7 and 11), as shown schematically in Fig. _ 
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FIG. 7.—Schematic Diagram Showing First R Mode of Earth Dam 


tes 


the horizontal sinusoidal forces from the shaking machines, each peak of which | 
is between 2,000 lbs (8,900 N) and 5,000 Ibs (22,250 N) and has a moment > 
arm of about 20 in. (51 cm) which produces a maximum moment at the crest _ 
of 80,000 1b-in.-200,000 Ib-in. (2,035 Nm-22,606 Nm). These rocking modes 
have never shown up in previous full-scale tests of earth dams (6, 10, 11, — 
and 13), but the finite element idealization developed by Chopra, et al. (7) 7 
and by Martin (11) has predicted this kind of mode in the U-D direction. wen?) 
The best estimates of the mode shapes of the crest vibrations at peak amplitudes 
s from the seismometer records are shown in Figs. 8 and 9 for the 
3-D response, and Figs. 10 and 11 forthe U-D response. 


Damping Values.—Modal damping values for all the measured U-D modes 


" were determined by the logarithmic decrement method, applied to recordings 
of the damped free vibration of the dam (immediately after discontinuing the _ 


power to the vibration generating system). This method was used because the oat. . 


width of the modal peaks of the resonance curves were generally not suitable _ 
_ for the determination of modal damping values by the half-power band width > 


i 
— 7 7 for the first measured symmetric mode (R1), or U-D motion associated with 1 
longitudinal (L) or vertical (V) vibrations. Peaks designated by AR are believed _ 
to be rocking antisymmetric modes. 
| Rocking modes could be a result of the vibrational moments generated by | 
| 
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method. However, pan Pron for a few modes where the spectral or modal q 
eae were well defined at the 0.707 points were calculated . Table 1 ‘nai 
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FORCED VIBRATION vests: ON SANTA FELICIA. 
UPSTREAM-DOWNSTREAM SHAKING (TWO SHAKERS, PHASE) 


2.270 


yations were 


ties from the resonance curves” 

= a. The force level used to excite this massive structure (it ethos about | 
.. 6,000 rely was very low, particularly at low frequencies. The development a 

w, large vibration generator capable of producing v very large 


pe | 
UPSTREAM - DOWNSTREAM SHAKING (TWO SHAKERS, 0°-PHASE) 
Wenn 
FIG. 10.—Measured S Mode Shapes Determined during Frequenc 
a U-D Direction): (a) Some Shear Modes; (b) Rocking Modes — 
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| 


TESTS 

FORCED VIBRATION TESTS ON SANTA FELICIA EARTH— DAM 


SHAKING (TWO 180" OUT-OF - PHASE) UPSTREAM-OOWNSTREAM SHAKING (TWO SHAKERS, 180° OuT - OF - 


id 


2.650 HZ 


gorse 
db 


= 


WA 


[preferably 1 up to 100,000 Ib (445,000 N) at of 0.3 Hz-2. 0 
for full-scale testing of dams, nuclear reactors, and other massive structures — 
} The peaks of some modes interfered with the response at other modes 
.! of the structure. This difficulty would have been largely overcome by measuring 
the phase between the excitation and the response because measurement of _ 
: the phase allows the response to be separated into its in-phase and 90°-out-of- 
phase components, and if the latter is plotted versus frequency, each resonance 
of the structure becomes much more sharply defined. Unfortunately, this 
_ force- -response phase was not measured due to instrumentation difficulties. 
> _¢. Finally, the resonant frequencies of the structure were often too closely 
spaced (see Table 1), sometimes overlapping, to — identification of every — 
2. The mode shapes of Figs. 8-11 do not correspond to ideal versions s of 
“the sinusoidal mode shapes because of | the varying width of the valley in which — 
“the dam is built and the nonhomogeneous nature of the construction materials. 
_ 3, The magnitude of the interaction between the shaker concrete block and a 
a surrounding soil in the dam crest was such that modal amplitude measurements 
_in the vicinity of the two shakers were distorted. Consequently, the modal Gc 
anaes of Figs. 8-11 at shaker Stations E2 and W2 were modified (reduced) — 


x 

| 

a 

| i FIG. 11.—Measured AS Mode Shapes Determined during Frequency Sweep Test _ 


by com comparing amplitudes measured in vicinity of the (during 
3, a tests—see Ref. 1) to those from later measurements on the downstream — 
and at the other adjacent stations on the crest, where the influence of dl 
block-soil interaction was negligible. ‘Measurements taken at distances greater — 
than 20 ft (6 m)-30 ft (9 m) from t the shakers did not appear to be ‘be influenced — 
4, An interesting feature of the data in Table 1 is the closeness of | cach 
2 S resonant frequency to a corresponding AS one. For example, the following 
.3 modes are very closely spaced in the frequency domain: RI and AS1, S2 and 
~ ARI, R2 and AS2, S3 and AR2, R5 and AS3, S4 and AS4, R6 and AS6, and 
finally S7 and AS7. similarity i is found between the U- D 


__ 5. The measurements did not permit determination of the actual mechanism 
of energy dissipation in the dam. Some of the energy anpertes to the dam 
by shaking was dissipated by radiation into the foundation, abutments, and 


water in the and some was consumed by hysteretic in the 


TABLE 1.—Measured (Average) Resonant {in Hertz) and Estimated 
Ratios from the U-D -D Shaking ot the Santa Felicia Dam 


Force per | Average from 


in Hertz i pounds vibration 


| 

= | = | are wr 

= 


The damping valees (of Table 1) were estimated by assuming that the eo ae 
. is basically linear. In order to estimate the damping of a nonlinear dam- a 
: - structure, the exciting force should be plotted against the structural Gieplncement 


the area inside the resulting hysteretic loop i is proportional to the energy dissipated 
- phase was not measured i inthe forced vibration test. 

_ 6. Finally, it is seen from Table 1 that the damping values estimated from 
ims resonance curves are higher and so less reliable (because of the interference 
of the in-phase and 90°-out-of-phase modal components mentioned previous) 

than those determined from the logarithmic decrement method (seen | to vary Es 
from 1.8%-3. 2% of critical damping). In addition, the damping ratios do not _ 

increase with frequency (they rather decrease a little). This may suggest that . 


Vibration ae 
During” the L shaking (last of testing) one of the = 
; aeieee so the other one was used alone; also, time was limited due 


TABLE 2. Frequencies (in Hertz), Observed i in L Direction of Santa Felicia” 


Dam (L Shaking Using One Shaker) Great dae 


in Hertz Designation Total f force, in pounds 
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VIBRATION TESTS 


‘to n of the test permit, so the measurements of 
in earth dam research) were not extensive. Nonetheless, the resonance — 

_ and the preliminary modes (of about half the dam crest) were useful. Table 

2 summarizes all resonant requencies of major peaks recorded at selected stations; 

_ because only eight seismometers were used during t this shaking, it was difficult 

toco completely determine several modes corresponding to the resonant frec requencies 
of Table 2. In addition, it was not possible to distinguish between S and AS 

_ mode shapes. However, Fig. 12 shows estimations of some of the modes (up ; 
to about 3.0 Hz) obtained during the L frequency sweeps (the solid lines connecting 

the data points of Fig. 12 are estimates of the modal configuration, while the _ a 
dashed lines | represent possible extrapolations); the local ‘magnification effect 
of the soil surrounding the shaker block is also shown. ‘Rey ais 


_ From Table 2 and Fig. 12 the fcllowing observations can can be made: 
There ae many closely spaced resonant frequencies in i 
And, some L resonant frequencies are very close (even identical) to some of © 
“the U-D frequencies. Examples of these are: L2(1.675) and $1(1.635); L3(1. 850) 
and R1(1.850); L(2.275) and R2(2.270) and AS2(2.30); L6(2.625) and R3(2. 60) : 
and AR2(2.650); L7(2.875) and R4(2.920); and L8(3.05) and AS3(3. 100). 
_ cies separated by 0.05 Hz or less are defined as “‘close.’”) This proximity may — 
= a strong coupling between the U-D and L vibrations, or it may suggest | 
that due to both the eccentricity of the single shaker (it was not located on _ 
the L-axis of the dam) and the fact that the dam is not symmetrical, the U-D 
modes containing significant L motions were excited. 
2, Some of the incomplete L mode shapes are similar to the L components 
of the 3-D motions recorded during the U-D shaking; examples of these are: 
«LAL. 675) and $1(1.635), L3(1.850) and ASI(1. 875), L4(1.950) and $2(2.100), 
L(2.275) and AS2(2. 300), L6(2.625) and $3(2.82), L9(3.150) and AS3(3.100), and — 
finally L12(3.775) and both $4(3.540) and AS4(3.600). 
It is apparent from these observations that it was very difficult to tell, from 
5 4 the data available, whether a guts 7958 was purely L, associated with a U- da 
mode, ora3-Dmode. (herd aro 
| Fourier Amplitude Spectra (FAS). _The occurring of the 
dam (U-D as well as 3-D motions) caused by strong winds (of speeds ranging 
# from 10 mph-40 mph) and by the spilling of the reservoir, were measured during 
4 the ambient tests. These ambient tests were conducted over 7 days, which — 
allowed the stationarity of the records to be examined. Samples of the Fourier 
amplitude spectra (FAS) of the recorded motions (at different stations and on © 
_ different days) are shown in Figs. 13, 14, and 15. These figures have been a 
chosen as best exemplifying the basic trends seen in all the analyses and 
- comparisons (for more detail see Ref. 1); the frequency resolution in all the | 
FAS was 0.0073 Hz which proved to be high resolution. Phase, cross power, = 
and cross correlation spectra (obtained by the digital signal processor) between - 
_ different stations were used to determine the sign (in-phase or 180°-out-of: ——— 
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frequencies for the dam i in the orthogonal directions; although the dif fi 
of the frequency values in this table are small, they are distinct. = — 7 
also contains the estimated damping rs ratios for the D ‘Grection. 


and ‘the w water flowing over the spillway. As will be covered in the con ‘companion — 

_ paper, the natural frequencies and modes of vibration are influenced by the 
‘point of application of the loading (as well as by its magnitude). Also, changes _ 
in frequency due to structural nonlinearity are quite noticeable when comparing — 

- ambient with forced vibrations, but from one ambient level of excitation to 
another these changes are relatively very small (less than 10%). 
2. The ambient vibration mode shapes were determined in the same way 
ss were the forced vibration mode shapes, and these modes were almost identical — 
_ to those resulting from the forced vibration tests, so only the naturai frequencies — 
(which did vary with the level of psy sepia of the ambient modes s are presented 


orthogonal directions (U- D, L, and V); i.e., during the ‘ambient 4 


= 
| _ observations can be made: PRatt 
1. Most of the FAS contain peaks 7 
some exceptions. One possible explanation is 
a were excited in varying degrees from one test (or station) to the other because 
4 q 
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AIG. 15.—FAS of L and V Motions Recorded Simultaneously a Same Station — 


« excitations the dam sometimes exhibited motion which dominated all oe 
_ orthogonal directions (Fig. 13). The V component in most cases shows a more 


complicated spectral pattern than the two horizontal 
several major in to those common to the horizontal Girections 
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tive value of 1.61 Hz (given in Table 3). This is pore sg to the value 


1.64 Hz obtained using forced vibration (see Table 


_ §. In some instances, lower U-D modes are very distinct while higher U-D 


modes are not clearly | discernible from other nearby (Closely spaced) call 


‘TABLE 3. iia Frequencies (in Hertz) and Estimated Damping Ratios (for U-D © 


Direction) from FAS of Ambient Vibrations in Three Orthogonal Directions (Santa 


Average 
“measured 


frequency, 
Hertz > 


AS6 
R6 


«Direction 


| Mode — 

designation 
according 

to forced vi- 


Average 
measured 


Average 
measured 
damping 
ratio, as a 


possible 


frequency, 


= 


V Direction 


Average 
measured 


= 
i FAS of ambient vibration show an exceptionally clear picture of the U-D “shear” — 
: s first mode (S1) (Fig. 15). Although the values (from different days) for the | 
q 
elicia 
percentage! bration tests | in Hertz 
3.0 
4.22 ALO 3.04 


in the ambient vibration | tests, there was both strong modal coupling 
‘the fundamental frequency (Sl) ‘difference between ambient 
(1.614 Hz) and forced vibrations (1.635 Hz) might suggest that the forces generated ee 
= wind [from 10 mph-40 mph (16 km/h-64 km/h)] on the dam’ s downstream “ Ta 
_ face were larger than those generated by the shaking machines; in addition, - — 
distribution of the forces was completely. different in the two (For 
detail, see Ref. | and the companion paper.) 
_ 7. Forthe rocking modes, several natural frequencies have ‘aint ie ae 
- peaks in the FAS of the — excitations (Fig. 14), unlike those of the 


Damping Values from FAS. —Estimated damping values from the FAS of the 


of half-point width was possible. (The width of the resonance peak can be 
influenced by several factors including: the nonstationarity | of wind loads, a 
presence of adjacent closely spaced peaks, the smoothing process of the FAS 4 
and the structural nonlinearities.) The modal damping values given in Table 

3 have values ranging from 2.3%-3.6% for the first symmetric U-D shear mode. oa 
a Even for an earth dam, these are relatively large values of damping considering e " 
5 - that they were determined from ambient level vibrations. In general, the results 

- show that the damping values | do: not change greatly with a change i in the — 


4 
4 


“increases; i.e., » higher modes have lower damping values. 


The test of the dam re: to small forces was accom- 
plished by use of pressure waves, impinging upon the upstream face of the -/ 
: dam and generated from a controlled, submerged release of gas under pressure ~ f 
in the reservoir water. The purposes of this test were: (1) To determine the x 
of the popper method of testing (Ref. 12) o on a large earth 


4 characteristics. Three sets of tests were performed; in 1 popper - test I the source 
' “of energy release was placed 10 ft-20 ft (3 m-6 m) away from the dam in <a 
q - the central region. In popper test II, the source was moved along a line in — ; : 
: the central region of the dam perpendicular to the longitudinal axis, and energy 
was released at various points from 10 ft-60 ft (3 m-18 m) away. In popper 
test Ill, the boat which released the popper traversed half of the dam length ae 
(starting from the east abutment), and energy was released at eo 4 
During all three tests, some seismometers were oriented. (at 


m amount of information about 


_ as seen in Figs. 13(a) and 15(a). In other instances, the opposite is true. In “4 ag 

j addition, in some instances, lower U-D modes are very distinct from those Bs 

s > lower L and V modes, and in other instances lower U-D modes are not clearly alt a 

@ 

| 

ang are shown in iabdie 5. Vespite some difficulties in measuring Dand ths, 
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-FIG. 16.—(a) Schematical Diagram Showing Excitation of Dam’s 4th S Shear Mode, 
_ $4 (Popper Opposite Dam’s Center Line); (b) FAS of Velocity Proportional Response 
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the modal configurations. Also, by loading ‘mostly the central 
a: stations WI and El) of the dam’s upstream face (tests I and Il), 
it was hoped that generally symmetric U-D modes would be excited. For the 
tests, it was intended to excite the first symmetric shear mode exclusive of 

all other modes; instead, , the fourth symmetric shear mode, 84, of of frequency 5 

3.27 Hz [f, = 3.54 Hz from the forced vibration tests and f,= = 3.665 Hz a 
(average) from the ambient vibration tests] was excited, almost exclusive - 
all other modes, as shown in Fig. 16. Theoretically, to excite any mode exclusively, 

the distance from the source of the P-wave to the dam must be adjusted such a 
that the pressure distribution would approximate the mode shape to be excited 
(10). It seems that the ‘resulting wave-pressure distribution on the upstream 
face approximated the S4 mode; the area of the upstream face near the center 

A line was loaded more than, and in an opposite direction to, the area off to 
either side of stations E6 and W6. However, neighboring, closely spaced modes 
_ were excited as well because of the proximity of their resonant frequencies 
to that of S4 as well as because hnipte setuid the pressure on the 


upstream face. 


18.—Sample Traces from Made Simultaneously during 
Pepper Test Il at Stations E1, £7, a and E9 in in U- ‘D Direction 
_ Additional modes were excited during popper test II; note the v ~~ defined 
mode RS high peak and the high contribution of modes AS4, AS3, R4 and 
$5 [Fig. 17(@) compared to Fig. 16(b)]. One possible explanation of the strong 
appearance of mode R5 is the movement of the popper during this test; possibly, e 
previous apparent overlapping of with $4 in Fig. 16(6) was 
eliminated to some extent in Fig. 17(a). ps 
In popper test III several modes were excited with each release of bite Ps 
particularly the AS mode ASI which was obtained by placing the pressure 
wave source opposite the quarter point of the dam’s crest [see Fig. 17(b); _ 
note that the results of this figure were smoothed with | one more cycle of — 
_a Hanning window (1/4, 1/2, 1/4 weights) than those of test II (which onal 
only one smoothing cycle)]. The recorded time traces of the popper tests were 
not clear enough to allow meaningful visual interpretation; only the S4 ror 
perhaps the R5) mode shape, » frequency, an and damping value could be determined © 
of | the U-D 


a 
| 
| 

- \ > 


: | from the FAS curves. Note that the frequency resolution was taken to be 0. _. e a 
Hz for the three tests, and the filtering frequency of the recorded — was 
>... In conclusion, the popper test could be an effective dynamic test of alee, 
large earth dams if more - experimental ‘studies were done, emphasizing the 
- relationships among the pressure-wave distribution on the upstream face, the 
_ distance of the source (or popper) from the structure, the slope and roughness — 
TABLE 4.—Resonant Frequencies (in Hertz) and Estimated | Damping Ra Ratios from 
_ Popper T Tests (U-D Direction) Santa Felicia Earth Dam 2! + 


Estimated Damping Ratios 
From FAS, as a percentage 


Average value x 
(8) 


be 


4 
>. 


of the upstream face, and the 1e length of the dam crest. However, this first 
trial of the popper test technique « on an earth dam | provided 


useful 
anda beginning forfuturetests, 
| ow From the information outlined in this paper and from extensive analysis of 


e full-scale results, the pmay b be drawn. 


| 
Average Measured Frequencies | 
| 
| | Test ill |Average| Test | Test Ill 
195 | 192 | — | 1.94 260 | — | 260 
| 2.06 | 2.02 2.04 | 207 | 343 | 2.75 
213 | 214 | — | 2 | — 
AS2 236 | 231 | 235 | | 1.73 | 330 14252 | 
| 249 | 252 | 2st | 150 | 235 193) 
~AR2| 266 | 260 | 266 | 266) — | — | 
(3.25 3.29 | 3.27 3.27 — | 2.83 2.83 
| 3.46 | 3.36 | 3.44 | 342 | 120 — | 12000 
— 350 | 359 | 355 | — | 240 | 2.40 
3.94 | 389 | 400 | 399) — | — | 
| 4.09 | 4.09 | 4.07 | 4:08 | 1.00, 
$7 «| 452 | 451 | 444 | 449 | 110 | 180 | 145 
AST | 4.66 | 462 | 461 | 463 | 0.90 | — 09000 


| 


= ‘Ati is that ambient vibration tons should on 
important earth and concrete dams in areas susceptible to earthquakes in order 
to roughly estimate the structures’ dynamic properties at this low level of excitation _ 


prior to any high-level seismicevent. = 


: 4 3. The results obtained via the popper test proved to be very informative. 


on Values of the observed resonant t frequencies: from the forced vibration 


= * For the geometry of the dam studied, the shear modes contain significant 


4 L as well as V deformations, so mathematical models to represent the prototype 
realistically in both L and V directions are 


6. The L vibration data obtained from both and 


7. Further work is —— to better understand the energy dissipation mecha- A 
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COMPARATIVE Stupy OF Dynamic RESPONSE 


a 
By Abmed M. Abdel- Ghaffar,' A. M. ASCE and Ronald F. Scott,’ M. ASCE His 


which data ‘for ‘this one conducted ‘not just to "investigate the 
_ Santa Felicia earth dam itself, but rather to gather experimental ‘data concerning ~ 
its low-level dynamic characteristics which could then be used to compare with = 
_ the relatively high-level characteristics resulting from the analysis of the dam’ s 
___—‘ Structrual performance during two actual earthquake conditions (3,4). In addition, , 
: the data was gathered to be used to test and develop various analytical an and = 
: f ‘numerical ‘methods for computing the natural frequencies and mode ee of r 
_ earth dams in general, for the purpose of predicting earthquake responses. = 


" a This paper, the second of a two paper sequence describing the results and 


from the ambient, and popper vibration tests. (in the companion 
a with those predicted by an existing two-dimensional (2-D) shear-beam 
_ theory (5,6,9). Furthermore, the dynamic characteristics estimated from the 


measured response of the dam to the dynamic tests are compared with those — 


_ Also in the present paper, the shear strains and shear 1 moduli are estimated ~ 
from the full-scale tests and are plotted along g with those estimated from the — 
earthquake response analysis (4). The same procedure was followed for damping _ 
_ versus shear strain at various levels. Then the combined results are interpolated 
_and extrapolated to provide a general guide to the material properties of rolled-fill — 
earth dams for both earthquake and geotechnical engineers. Additional details 


TASS. Prof., Civ. Engrg. Dept., Princeton Univ., Princeton, N.J. 08544. 

?Prof. of Civ. Engrg., California Inst. of Tech., Pasadena, Calif. 
_ -Note.—Discussion open until August 1, 1981. Separate discussions should be submitted 
‘agit for the individual papers in this symposium. To extend the closing date one month, a 

3 a written request must be filed with the Manager of Technical and Professional Publications, & 
ASCE. Manuscript was submitted for review for possible publication on March 17, 1980. * 
Ss paper is part of the Journal of the Geotechnical Engineering Division, aa _ 

of the American Society of Civil Engineers, ©ASCE, Vol. 107, No. GT3, March, 1981. 
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| of the study presented herein can be found in report (1). 4 
_ Comparison Between Results of | U-D Shaking And Those of 2-D Shear- Beam ile 
' Theory.—To estimate the shear-wave velocity within the dam material from 
the observed resonant frequencies and to check the values of the measured i 
7 resonant frequencies, a 2-D homogeneous medium shear-beam theory was used. 


The natural frequencies and modes of shear vibration determined aaead this -* 


are given by (56,89) bg 


y= = the (n, r)th natural frequency; v, = the poe wave velocity -p 
within the dam where G = the shear of the dam 
and p = the mass mayrsntde A= the yee of the dam; B,, n = A, 2, 3, 


For symmetric modes, n = 3,4 5, 

antisymmetric (AS) modes, .. and r = 2, 6, 


a this 2-D theory of a triangular voter ina rectangular canyon, the trapezoidal 


equal to the average of the: crest ‘length and the length of ‘the base, i.e., 

‘1 = 0.5 x (1,275 + 450) = 913 ft (278.3 m). Substituting this value of /, and — 
values of Ah = 237 ft (72 m) and f,, = w,,/2m = 1.64 Hz (the measured — 
frequency of the first U- D shear mode from the forced vibration test) in Eq. - 
A Bives a value of , equal t to » 957 is close t to the 


that there is a variation as the depth: changes, but that 

a range of v, between 850 fps (259 m/s) and about 1,070 fps (326 m/s) is 

representative of the dam material. Aine, by substituting the value of the measured 
AS shear frequency = = 1.875 Hz in Eq. 1, the resulting 


vy, is equal to to ‘959 fps (292.5 m/s) which is very close to that from ‘the first 
The values of higher frequencies (>f,. therefore 
by substituting the value of v, = 958 fps (292 m/s) in Eq. 1; these computed 
natural frequencies are listed in Table 1 where they are compared with those > 
observed from the U-D forced vibrations. Only the frequencies of the observed 
_ shear modes (taken from the companion paper, , Ref. 2) have been shown in } 
Table 1 (‘‘shear’’ modes are those which are similar in configuration to the ~ 
— On me basis of the information contained in Table | as well as in the resonance 
curves, the plots of U-D mode shapes (both shear and rocking modes), and a 


the plots of 3-D components of shear modes (all in the companion paper) the the 


- following conclusions may be drawn: 


| 


1 
J 
| 
| 


by Eq. 1 (the are than the computed 
j “- requencies by from 3%-17%, but overall it is fair to say that the correspondence | 
‘. between the observations and Eq. | is good over the entire frequency range 
_ evaluated. The computations ; were based on the assumptions — of an empty 
h reservoir, a rigid foundation for the dam, a rectangular canyon 913 ft long, C 
an average value for the shear wave ‘velocity of 958 fps (292 m/s), and only d 
elastic shear (low-strain level) deformations in the horizontal direction (U-D). __ 
_ 2. Because of the selected deployment of the two shakers, it was possible 
to excite pure first (1,1) modal response (response characterized by all points © 
in the dam moving in the same phase). Several t trials were ¢ made in in the field > 
to measure modes (2, 1) and (2,2) which have similar configurations to the first — 
S mode (1,1) and the first AS mode (1,2) along the crest and have a node 
along a vertical (V) section [like modes (2,3) and (2,4)]. Unfortunately, these 
‘ - modes were never excited or measured during the forced vibration tests. Nor 

e these modes ever been measured and reported vaste the literature (6,8 9). 
TABLE 1 _—Comparison between Observed Resonant Frequencies (from U-D Sheking) 
_and Those Computed by 2-D Shear-Beam Theory (Eq. 1): Santa Felicia Earth Dam 


AMTISYRMETRIC VIBRATION 
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2273 
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‘een a one in phase with the « exciting force, and one one at 90° out-of-phase J 
—. to the exciting force, and if the dam were to behave as a linear ‘structure, — 
_ then all of the response of the resonating mode would appear in the epson 
at 90° out-of-phase to the exciting force. However, as a result of damping, 
all of | the other modes of the dam would respond at phase angles other than 
0° or 180° to the exciting force, and therefore they would appear to some 
extent in the component of response at 90° out-of-phase, depending on the 
_ amount of damping, the spacing of the natural frequencies, and the location 
of the exciting forceonthedam. 
4. The shear- beam theory (low levels of in ressonably good 


/ _ Instrumentation which would accurately (and less laboriously) measure the phase § 
: between the excitation and the response would be of great help in exploring 
modes, indicating the inadequacy of the theory for earthquake response (high 
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EARTHQUAKE RECORDS 2-D Shear Beam | 


Earthquake, Earthquake, | First Frequency 


Observed Frequency | Measured | Measured 
frequency, | frequency, - frequency. 
in Hertz | e in Hertz 


AR2 
- 


ASI, 6) 


15 
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AS7(1, 10) 
AR4 
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= of strain) computations. In order to dhiiie Seinemnent between the 
c, computed and all of the observed frequencies, it appears that a representation 
= the dam more sophisticated than that provided by shear-beam (or wedge) — 
theory is needed, and possibly a more accurate knowledge of the physical 
properties of the soil in the dam is also required. Perhaps more accurate dynamic 
_ properties, using a 2-D shear-beam theory with nonuniform soil properties = 
the dam depth, may be attained by assuming for instance, that the shear modulus" 
G varies with depth (say, G « y'” or G « y*/*); this has not been tried 
. _ 5. From the shapes of the first four measured S (and AS shear modes at 
a amplitudes, obtained from the simultaneously recorded three components — 
_ (U-D, L, and V shown in Ref. 2), it is seen that the U-D modal motions are > » 
_ accompanied by contributions from the longitudinal (L) as well as the V motions; : 
these contributions range between 10% and 25%. For S vibration, the L and | 
V motions along the crest are similar to the U-D configuration, while along 
the downstream slope there is some disagreement. For AS vibration, the 
configuration of the L component is similar to that of the U-D but the V component 


is completely different along the ¢ crest. And again, there are ; disagreements among 


nation 
1.66 1.850 | Ri |Rocking) —_ sho 
— | — | 2600 | R3 |Rocking} — 
2.88 15 | 27 2.840 | — 
2, «3.33 3.550 |S42,3)|Sher | — | | 
“= 34 = 3.600 | AS4(2,4) | Shear 
3.93 2, 3.68 | — | — | 4100 | 8 Shear 
- 1, 4450 R6 | Rocking | ot | = i 
i 
f=: that: although the dam is costly, it would — a 


 2-DShearBeam | AMBIENT VIBRATION ME MEASUREMENTS _ «ped Average Measured 
g Theory (Based on Frequencies From 
Observed First Measured Frequency From Fourier Amplitude 
d measured; 

Frequency, 1.635 Hz) Fourier Amplitude Spectra, in Hertz | Spectra, in He Hertz 

Frequency | quency, 


Hertz 


Average 
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Between Results of Full- Scale Dynamic Tests and Dam’s Earthquake 
_ Responses. — Tables 2, 3, 4 contain summaries of comparisons made between 
and V directions, respectively) — 
in the dynamic | tests and those the spectral analysis of the 
dam’s two earthquake responses (3), 
_ The time history as well as spectral analyses of the dam’s earthquake responses _ 
(3) indicated that the structure responded primarily in its fundamental U-D shear | 
- mode with apparent natural frequency of 1. 45 Hz (0.7 secs). However, the 
_ response in the L direction showed significant contribution from higher as well © 
as lower modes; the V component indicated a very complex response. Further- 


dam undergoes rather small relative displacements of its crest at a moderate 
intensity of shaking; the peak accelerations recorded during the 1971 earthquake — 
were 0.22 g on the abutment and 0.21 g on the crest. | we. Perea fi 
_ = U-D Direction.—Based on the results of Table 2 the following observations 
- relating to the U-D motion from the dynamic tests and from the earthquakes _ 


' isting Shear-Beam Theories (Santa Felicia 
— wo ow) sit} mab sd 
13/ | 14/ | 15/| 16/ | 17/] 24/ | 27/ | inHertz | Test | Test | Test | er 
nr 78 | 78 | 78 | 78 | 78] 78 78 | “7 days” | m | age 
_ (13) | (14) | (15) | (16) | (17) (18) (19) | (20) (2%) | (22) (24) 
1,1 | 1.635 |1.620] 1.621 |1.586] 1.619 |1.603) 1.620 | 1.626 | 1.614 [1.513 | 1.503 | 1.512 
— 1.753} 1.779 |1.792| 1.773 |1.745] 1.787 | 1.780 | 1.773 |1.727 | 1.719 | 1.677 708 
1.884 ]1.924) 1.835 ]1.916] 1.881 | 1.926 | 1.900 1.892 | 1.802 | 1.831 | 1.862 
— | _—  |2.025} 2.020 2.030] 1.978 | — | 2.035 | 2.004 | 2.022 |1.954] 1.922) — | 1938 
2208 2.130 | 2.130 |2.097] 2.044 |2.050 2.119 | 2.097 | 2.062 | 2.027 | 2.027 | 2.037 
— 12.224] 2.244 12.168] 2.161 12.225] 2.237 | 2.229 2.133 | 2.141 | — | 2.137 
q 1,4 | 2.610 |2.368] 2.350 |2.395| 2.333 [2.377] 2.367 | 2.320 | 2.359 [2.355 | 2.307 | 2.354 | 2.339 
| 2.611 ]2.549] 2.540 2.568 | 2.500 | 2.564 [2.489] 2.522 2520] 2510 
| |2.680} 2.693 |2.680| 2.695 |2.632| 2.699 | 2.646 | 2.675 |2.659| 2.604 | 2.659 | 2.641 
3.048 |2.900] 2.809 |2.783) — |2.842| 2.869 | 2.859 | 2.846 |2.778 | 2.805 | 2.735 | 2.773 
— 2.942 |2.903] — 3.038} 2.970 | 2.989 | 2.975 |2.954| 2.920 | 2.897 | 2.924 
4,6 | 3.510 |3.091] 3.230 [3.128] — [3.190] 3.214 | 3.135 3.093 | 3.057 | 3.062 | 3.071 
3.595 [3.157] 3.348 ]3.230) — |3.289) 3.105 | — | 3.226 | 3.186 | 3.157 | 3.157 | 3.167 
| | — | 3.606 [3.474] — 3.539 | 3.582 | 3.536 |3.246 | 3.289 | 3.270 | 3.268 
| 3.891 | — | 3.745 3.641] — ]3.515] 3.708 | 3.669 3.665 | 3.457 | 3.362 | 3.435 | 3.418 
3.987 | — | 3.943 [3.853] — |3.847] 3.866 | 3.830 | 3.868 |3.569| 3.496 | 3.589 | 3.551 
— 24 4131 | — | 4.146 [4.136] — |4.124) 4.181 13.206 | 3.675 | 3.735 | 3.705 
— |} — | — | 4.243) 4.206 | — | 4.221 |3.935 | 3.890 | 3.999 | 3.961 
1.8 | 447% | — | 4.360 |4.282| — |4.339) 4.325 4.327 | 4.093 | 4.093 | 4.065 | 4.084 
26 | 4752 | — | 4672 — |4.433| 4.476 | — | 4445 | 4.230] 4.256 | 4.292 | 4.259 
1.9 | 49m | — | 4596 — | —| — | — | 4555 | 4.400] 4366) 
4.695 |4.656] — | 4.636) 4.598 4.646 | 4.517] 4.507 | 4.444 4.489 
1,10 | 5.470 | 4.780 | 4.749] — | 4.804) 4.790 | — | 4.781 | 4.662 | 4.616 | 4.607 4.628 
—| 4973 | — 4.973 | 4.745| 4.693 | 4.706 | 4.715 
| — | 5.140 5.140 | 4.851} 4.831 | 4.867 | 4.850 
the: indicated that a modern rolled-fill earth 
— 
— 


1981 


shear various full-scale tests that estimated from 
the response of the dam during the two earthquakes, it is seen that the natural = 
period was longer during the earthquakes (/, = 1.45 Hz or T, = 0.69 sec) = 


2 
a a) | 


San Fer- 
nando “fornia 


quake "Tests lIENT VIBRATION MEASUREMENTS 

Frequencies From ¥ frequen- 


Spectra, in Hertz 


Possible 
Possible 


Possible 


than it was during the tests (forced: f, = 1.635 Hz or T, = 0.61 sec; ambient: i 


SL a. ¥ 614 Hz or 7, = 0.62 sec; popper: f, 1.509 Hz or 7, = 0.66 sec). 


TABLE 3.—Comparison between Resonant Frequencies from Earthquake Records 
Full-Scale Dynamic Tests (Santa Felicia Earth Dam): L Direction 
§ 
oe 
| = Hertz | Hertz | Remarks | 3/13 | 3/14} 3/15 | 3/16 | 3/27 days” 
| 1.27 | 1.50 | 1.42 | 1.47 | 1.41 | 1.43 
1.707 166 | 168 | 1.67 | 1.63 | 1.66 | 1.58 | 1.65 | 1.64 
1.86 | 1.85 1.88 | 1.88 | 1.87 | 1.81 1.80] 1.85 
| 2.03 | 1.98 | 1.95 | 2.04 ]1.92 | 1.99 
2.15 | 245 | 216 | LS | 2.23 | 2.15 | 2.18 | 2.16 216 | 2.18 
| 228 2.36 | 2.36 | 2.33 | 2.24) 2.27] 232 
x, |,74 | 2.44 | 2.48 | 2.49 | 2.40 | 2.45 | 2.45 
— 2.64 2.63 | Lo | 2.65 | 2.66 | 2.62 | 2.67 | 2.60 264 ; 
{2a 2.87 | 2.93 285° 2.94 | 2.81 | 2.88 
| — | 3.06 | 3.08 3.01] — | 298] 303 
| 3.22 | 43.23 | Possible | 3.21 | 3.21 | 3.24] — | 3.26] 3.23 
3.35 | LIO | 335/333} — | — |338] 332 
30 | — | 34 Li 3.49 | 3.49 | 3.46 | — | 3.50 | 3.49 
3.85 | 3.71 | 3.78 Li2 | 3.75 | 3.80 | 3.82 | 3.70 | 3.81 | 3.78 
| 4.06 | 4.04 | 4.07 | — | 4.03 | 4.05 
4.20 | 14.18 4.17 | 4.22 | 4.16 4.22 | 419 
437 | 437 | 433) — [437] 436 
4.70 | Lid | — | 464} 470) — | 466] 467 
488 | — | 488 | LIS — — 488 
i i _ — 


— 


TABLE 4.—Comparieon between Resonant Frequencies from. Earthquake Records 
and Ambient Vibration Tests (Santa Felicia Earth Dam): V Direction vedo ae 
AMBIENT VIBRATION MEASUREMENTS 


frequency, 


3/14/ | 3/18/| 3/16/ | 3/27/ | inHertz 
78 | 78 | 78 | 78 | 78 days” 


3. 92 


=F 
— 
| 1.24] 1.23 | 125 | 126 
147 | 148 | 143 | 144 | 1460 
| 1.63 | 1.66 | 1.63 | 163) 1.64 
198 1.96 ’ 193 | 204 | 198 
we 2.19 | 2.22 | 239 | 223 | 229 | 226 
14264 | 263 | 265 | 265 | 263 | 264 
3.22 | 3.22 | 3.24 | 3.19 | 3.19 
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q The lengthening of the first S shear period during earthquakes by factors of a 
_ 1.13 (as compared to forced tests), 1.12 (ambient), and 1.04 (popper) corresponds» 
q : to a decrease in overall stiffness s by factors of 1.28, 1.25, and 1.08, respectively. — a 
This behavior is typical of a a yielding dynamic system such as the soil c 
an earth dam which is known to be extremely nonlinear. However, in vesewing 
57 this variability of the dam’s dynamic characteristics one has to be very careful 4 
about concluding whether these changes in frequency are caused solely by a 
change in the stiffness due to different excitation levels of the dam material — 
or by other factors | as well. For example, this » variability could be partially \ 
due to an increase in the v vibrating mass of the dam caused by an increase 
in the level of the reservoir water, or it could be partly due to the | nature 4 
of the dynamic excitations which varied greatly in their time history characteris- - . 
tics, spatial distributions, and intensities. = 
a During the full- scale measurements, the water level i in the | reservoir was at 


‘the reservoir level was s quite high but not at its m maximum level. The virtual = 


additional mass could have had some effect, but it is believed that such effect 
. would be very slight and probably unnoticeable (because of the massiveness 
of the structure and the flatness of the upstream face, Ref. 7). Thus, the data _ 
presented here seems to imply that changes in dynamic characteristics observed - bie 
rece the tests were ‘mainly the result of the loss in stiffness due to the nonlinear _ 
behavior of the dam materials and less significantly due to the nature of the 
_ Observation 2.—The frequency difference between the first S shear mode 
t's the U-D direction from the ambient tests (1.614 Hz) and that from the forced 
_ vibration tests (1.635 Hz) might suggest that the forces generated by the wind 
e (from 10 mph-40 mph) on the dam’s downstream face were larger than those 
generated by the shaking machines. (During some forced vibration tests, the an 
wind prevented measurement of the steady-state response of the dam, but during 
; most of the testing this was not a factor). An estimate of the wind forces 
A acting on the dam can be made by using the dynamic pressure expression (usually — 
for calcalating wind loads on structures) p = 0. 00256 C, (in pounds 
~per square ft; 1 psf = = 47.9 Pa) in which C, = all = the shape factor; and V = 
the wind velocity, in ‘miles per hour. Now, assume C, = 1.0 (a convenient 
assumption for the rectangular horizontal projection area of the downstream | 
_ face) and assume an average velocity of 20 mph; then, the foregoing expression 
F gives p = 1.024 psf. The exposed effective area (horizontal projection) of of the 
_ downstream face can be estimated to be 100 ft (depth) times 800 ft (average 
length); this value produces a total wind force of F,, = 1.024 x 100 x :00 
... one assumes that the wind was not really perpendicular to the downstream — 
face, then the aforementioned value can be reduced by, say, , 40%; this gives 7 


a total wind force of 29,152 lb, , which is approximately five ; times larger than 
the average steady- state sinusoidal force. Thus, the total amount of energy 
_ imparted to the dam by the wind, very likely, was greater than that produced 
Finally, comparison also indicates that, in exciting the first S shear mode 
al g the force generated by the popper (at 1.509 Hz) was more than that generated 


during the ambient and mechanical- ~shaking tes Unfortunately, there 


4 
were 
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no dynamic pressure transducers mounted on the upstream face to provide 


the order of magnitude of the popper forces. 
_ Observation 3.—During the two earthquakes, very | few high modes (shear — 
or rocking) were excited in the U- D direction (the response was primarily from — 
the first mode), whereas, contrarily, the higher modes were excited efficiently 
and clearly by the three dynamic tests. In addition, the identification of some 
of the higher modes excited during the earthquakes was not reliable. As = 
result, it is difficult t to mak make site reliable comparison based on a ' frequency-to- to- 


Observation 4. — Looking at the of the first Ks mods, Al ASI, 

81 Hz “1971 earthquake,” 1.86 Hz “1976 earthquake,”’ 1.875 Hz ‘‘forced, 

— 1.900 Hz “ambient,” and 1. 842 Hz “popper”), it is seen that there is a maximum > 

z= pone of only 5% which may suggest that the ASI mode may not have 


e ground shaking at the dam causing S modes to be excited more effectively; 
and (2) the location of the crest accelerograph in the central region (about = 
50 ft away from the crest midpoint), i.e., in a position more likely to pick — 
7 up the dominant s motions since this region represents nodes of AS modes. 
Observation 5. —Some “of the lower U-D rocking 1 modes were 1 not excited 
during the shaking of the two earthquakes; examples of these modes are ARI, 
R2, AR2, and R4. Apparently the first R mode, Rl, was excited during oe 
1976 earthquake (with frequency 1.66 Hz); this frequency was increased by 
_ 3% (1.708 Hz) during the popper tests, by 7% (1.773 Hz) during the ambient 
tests at and by 11% (1. -850 Hz) during the forced Vibration tests. ‘The | third 


2.47 Hz compared to frequencies of 2.510 Hz (2% increase) from the popper 
_ tests, 2.564 Hz (4% increase) from the ambient tests, and 2.600 Hz (5% increase) = 
from the forced tests. Finally, the resonant frequencies of the rocking modes 
ARI, R2, and AR2 resulting from the full- scale tests vary slightly from one 
| _ type of test to another; the variation ranges from 1%-5% as indicated by Table 2. ¢ 
Observation 6. ——The frequencies of S and AS shear modes S2, AS2, and 7 
ad -AS3 also vary only slightly from the earthquake responses to the full-scale 
tests. However, the third S shear mode S3 was excited more strongly during . 
the small earthquake of 1976 (frequency 2. m1 Hz) than during the large earthquake © 
of 1971 (frequency 2.88 Hz); the frequencies from the tests were 2.840 Hz, 
2.844 Hz, and 2.773 Hz from forced, ambient, and popper tests, respectively. 
- This just demonstrates that greater shaking does not always excite all modes 
ia Observation 7. —It is interesting to compare the frequencies of the fourth ie 
Ss shear mode, S4, resulting from the ‘different: tests and the two earthquakes 
cages the mode associated with this frequency was : excited very. effectively 
_ by the popper test (2). The frequency determined from the 1971 earthquake 
response was 3.470 Hz, then it increased to 3.610 Hz from the 1976 earthquake 
response. The popper tests provided a frequency of 3.268 Hz while the fi requencies = 
from the forced and ambient vibration tests were 3.550 Hz and 3.536 i -F 
“respectively. Note that the popper test excited the S4 mode more strongly than 


the strong earthquake. This observation is a surprising one since the 


eu 
| 

| 
| 

urn WO Ca uakes are presented in Kel. . Keasons 
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dam responded ‘mainly i in its shear mode during the tw two earth- 
quakes, while during the popper test the dam was excited, exclusive of all Ey 


en 
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FIG. 1 Shear Moduli and Corresponding Dynamic Shear Strains (Evaluat-— 
ed from Hysteretic Response of Two Earthquakes as well as from Full-Scale Dynamic : 
Fernando of Fe 197) (M, 
Forced Vibration Test meme) dor > 
hae 


FIG. 2. —Estimated Damping Factors and Dynamic Shear ‘Strains 
_ (Evaluated from Hysteretic Response of Two Earthquakes as well as from Full-Scale 


fl Observation 8.—The dynamic shear moduli and damping factors of the dam’s 


material were estimated from the measured responses (the hysteretic responses _ 


of the first S shear mode) to the two earthquakes. Both dynamic properties — 


| 


| 
| 
| 
~ 
| 
| 
| 
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"were ‘determined as functions of the induced dynamic strains (for more details — 


see Ref. 4). The relationships between each dynamic property and the induced 
dynamic strain are shown by the semilog plots of Figs. 1 and 2 for the first 


20 sec of the 1971 earthquake and the first 6 sec of the 1976 earthquake. ee ; 


_ In order to reveal any change in the dynamic properties of the dam, the 
eo shear moduli and damping factors from the full-scale tests are also 
shown in Figs. 1 and 2. In addition, Table 5 shows some of these dynamic 


. properties associated with various levels | of response. Both the figures and the . 
t 


able indicate that the dynamic properties of the dam’s constituent materials s i 


TABLE 5. —Comparison between Dynamic Properties (of Earth Dam) Estimated from _ 
rthquake Records and Those Determined from Full-Scale Dynamic Tests ha —_ 


‘Long- 

Source Estimated Estimated acceler- | Maximum Direction 
of data damping ation at displace- 

levels of frequency. a percent- crest, in 
excitation) in Morte age of g millimeters 
1971 San J Ranges from 
Fernando 
earthquake 


fornia 
earthquake 
M,=47 


Ranges from 
vibration 3.6 x 10°- 
tests (wind 43x 

and spill- 
ing of res- 


ervoir) 

‘Hydrody- 

mamically 
generated 


force 


(popper 
(tests) 


are consistent with those determined — 
ae the relatively | larger strains induced by the two earthquakes. part = 
i By interpolating and extrapolating the data of Figs. 1 and 2, Figs. 3 and — 


4, accordingly, are Proposed to provide _Teasonable estimates “of the dynamic 


_ for any study of the dam’s U-D earthquake-response characteristics. In addition, 
the variation of material properties with depth should be taken into account © 
_ for any realistic dynamics study; the variation of the shear modulus of the _ 
dam material with | depth is proposed in Fig. which i is based on a an in- in-situ 


q 
| 
< 
_ 
¢ 
1976 South- } Ranges from | Ranges from | Ranges from 0.85 1.86 
ern Cali- 10°-| 0.10 x 15-50 | 
Forced | 1.635 | 4.10 x 10° | 0.716 x 107° 0.0037 (0.00336. 
Wa C10 ncasulre a we a uc alla WO Ca UUGAC a 
a (3). Figs. 3, 4, and 5 should provide a good | al 
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of rolled- fill, essentially cohesionless material, with or without a relatively thin 


"observed during the two. corthanshes (at 1.35 Hz, 1.70 Hz, and 1.86 Hz), it 
S 6s iss difficult to make any reliable comparison with the values of the resonant 
eat 


n 


ANTA FELICIA 


(percent) 
FG. 3. Shear Moduli and Corresponding Dynamic Shear Strains 
for any Earth Dam Composed Predominantly of Rolled-Fill, , Essentially Cohesionless_ 
Material, with or without Relatively Thin Core pee 


percent ) 


santa FELICIA 


ARTH 


DAMPING RATIO( 


SHEAR STRAIN (percent) 


Material, with or without Relatively Thin Con 

frequencies in the L direction resulting from the forced vibration tests and 
those resulting from various ambient vibration measurements, as indicated by 

| Table 3. As was considered previously in the oma paper (2), there are | 


| _ properties in the dynamic analysis Of any carth dam Composed predominantly | =m 
. 
— 
| 7 frequencies are very close (even identical) to some of the U-D resonant frequencies | r 
“which may suggest a strong coupling between the U-D and L vibrations. 


first L mode was apparently excited more during. the small 
a, earthquake (where f = 1.27 Hz) than during the large (1971) earthquake (where 
i.35 Hz); the ratio is about 6% which corresponds to a decrease in eal _ 
bya factor of 1. 12. The frequency increased tol. 425 Hz when determined 


to an average value of 1.446 Hz when determined from the ambient tests (a 
14% increase over the 1976 earthquake). These increases of the first L frequency — 
correspond to decreases in overall stiffness by factors of 1.25 and 1.30 for 

the forced and ambient | respectively. This behavior is similar to that 


= 4, where while the comparison between the results of the two ie 
is quite consistent, the comparison between the two earthquake results = 


HEAR MODULUS x x 10 
y! 
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« results of the two earthquakes showed well-separated frequencies. = = 
_ Finally, there is also substantial indication of modal coupling and interference = 
between the three orthogonal directions; there were many with 
identical values m 
Concuusions, ano RECOMMENDATIONS 

appears s from all the data available on the dam’s responses, ‘including the 

1. Full-scale tests at levels of excitation much lower than those experienced 
the earthquakes produced s in the dynamic 


§. —Suggested Relati Depth 
a those of the ambient vibration measurements is not consistent at all. The results m a 


ae 


“of the da the dam. The behavior is typical of yielding dynamic systems stevvabuling soil, 
which is known to be extremely nonlinear, but has rarely been examined a at 
- _ 2. Dynamic properties of the dam’s constituent materials, , such as shear mc moduli ‘ 

and damping factors, estimated from low-strain full-scale tests and determined 
_ from the relatively larger strains induced by the two earthquakes are consistent 

_ with established shear moduli-induced strain and damping-induced strain curves, 

- Previously proposed ¢ on “eS basis of laboratory te tests. s. The lower levels. of excitation 


covered the lower of induced strains" in 1 these curves, ves, while the strong 
x ground motion during the earthquakes resulted in dynamic properties covering — 
the intermediate and higher-induced strains,§ 
The losses in stiffness due to high strains have been attributed, for the — 
“part, to the nonlinear behavior of the dam materials. An exact determination | 
7 the causes of these nonlinearities is almost impossible due to insufficient 
data on the influence of the foundation, the reservoir and the pore pressure 
- changes. These influences are complex and difficult to analyze accurately. Finally, 
the variations of these dynamic properties for this specific dam have been 
_ generalized to aid in developing analytical techniques to predict accurately the 
low-strain as well as high-strain level earthquake response for most modern — 
gi The correspondence between resonant frequencies from full-scale tests — 
and those estimated from the spectral analysis of the two earthquake records 
is reasonably good over the first few frequencies, but higher modes could not 
be reliably matched. . A comparison of these first few frequencies revealed that 
the natural frequencies of the dam increased significantly during the full- “scale 
_ tests in all three orthogonal directions (U-D, L, and V). Clearly, the dynamic 
_ characteristics of real structures are not actually constant as implied by linear 


4. _ During U- D shaking of the dam, the ‘Shear ‘modes indicated significant 


weaken the cohesion of the dam’s L ‘vibrations are also ‘important 
_ because they could cause transverse cracks across the axis. Thus, mathematical _ 
_(i.¢., analytical or finite element) models are needed to represent the prototype re 
in both and V Guections. These models should take into account 


good agreement with predominantly shear-like modes, do not predict all the “ 
U-D motion-type modes, indicating the inadequacy of these theories for compre- bs 
hensive earthquake response computations. The nonshear U-D modes are believed ie 
to be either rocking modes or indications | of U- D | motion associated with | . 4 7 

. 6. It appears that precise field tests and more theoretical analyses over a 

wide range of strains are required to establish the nature of the dynamic response 
characteristics of earth dams. Despite the complexity of several factors influencing 

the Tesponse of the dam, it is } believed that full- scale test t studies will yield 


wt Concurrently, it is s also desirable tor reassess the installation and deployment a 


| 
| 
| 
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reliable and informative results for interpretation of the dynamic behavior of — 
these structures during earthquakes. In order to get reliable information on 

‘i any dam’s 3-D behavior during an earthquake, more than two sets of strong- -motion | 
ncn should be deployed on and in the vicinity of the structure. Proper 
_ placement will yield information about the response of the dam, the nature 


of different modes of Rrenagys, Ti coupling of these modes and the nature 


A 
a. A minimum set of four three-component olla should be located — 
; ai on the crest and the downstream slope. One instrument should be placed os 
m onthe midpoint or the central region of the crest with a second on the midpoint — 
= lower than this point) of the U- D face, anda similar set of two accelerographs 
should be located on any given “cross section between the crest midpoint — 
a either of the two abutments. The records from these instruments would — 
F- help to identify the different modes of vibrations, 
2 ___b. Asset of two instruments should be placed on the right and left abutments _ 
7 (crest level) of the dam in order to correlate the ground motions at the two — 
sides and to evaluate any phase differences as well as to compare the — 
motion at the sides with that of the downstream foundation. lobarel ory. 
¢. A set of at least two instruments should be located in the U-D area. 4 
This set would provide information about the input ground motion and the 


of ground motion underneath the tory, 


q from the beginning of ground motion, after simultaneous triggering. Digital 
recording greatly eases data analysis. = 
ra Finally, it is believed that the analysis and the data presented here will 


dynamic characteristics of ‘earth dams and will ‘promote further research 


ams 
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The following are suggestions for appropriate locations of the instruments; 
ls f circumstan 
| 
_ U-D, L and V directions. In addition, records from the instruments should © 4 : 
a =. be time tied. Continuous digital tape recording with a several second storage ff 
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| ENGINEERING DIVISION 


(© jeonlev nt yousqstaaib (2) to 
PROBABILISTIC EVALUATION OF LOADS*® 


s + a geotechnical design, engineers have been primarily concerned with the 
_ determination of soil properties and the prediction of soil performances. Although 
also a major input parameter in design, loads are in most instances already — 
, prescribed to the geotechnical engineers. The same is also true regardin 
‘ ground-water level and other environmental f actors. Geotechnical engineers have 
_ long recognized variability in soil properties as exhibited from laboratory test © 
_ results and observed field performances. Significant amounts of uncertainty 
a exist in the loads and environmental factors. The uncertainties may arise 
_ from spatial as well as temporal fl fluctuation of load and environmental factors, 
and from idealized assumptions ne necessary y for ev: evaluating load effects. In view ; 
& of such uncertainties in both soil and load characteristics, safety of a geotechnical 
_ design cannot be guaranteed. Although a large proportion of the uncertainties _ 
in the load and environmental factors is due to inherent natural variability, 
which i is perhaps beyond the engineer’ s control, nevertheless, 1 an in understanding cam 


= factors to reflect the relative uncertainty levels of the individual loads 

a load combination, consistency in reliabilities could be maintained in the design. ; 

This paper reviews the loads that affect geotechnical designs, 
_ the uncertainties present in each load component. ‘Probabilistic ‘models for 
evaluating each component 


the level of variabilities and uncertainties associated wit with the effect of ‘each 


component on the foundation will be assessed. ead) 


_ ie Presented at the April 2-6, 1979 ASCE Convention and Exhibition, held at Boston, 
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Note. —Discussion open ‘until “August 1981. Toe extend th the date one month, 
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ASCE. Manuscript was submitted for review for possible publication on January 1, 1980. > 
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Gravitational load ae - the weight of structural materials such as concrete, _ 
steel, or form work is commonly referred to as dead load, whose ) magnitede 
usually does not change during the life of the geotechnical system. In general, — 

_ dead load could be estimated relatively accurately. Possible sources of uncertain- an 
ties include: (1) Variability in material density due to inhomogeneity and tolerance. 
of manufacturer; (2) discrepancy in dimensions from the | design values; (3). 
uncertainty in the final choice of building materials; a and (4) variability in the a 
nonstructural components such as decorative architectural forms. A coefficient 
of variation (COV) of 10% has been suggested (20) for the variation of dead = - 
load in the design of a typical floor in the superstructure. The total dead load 
applied to the foundation is generally composed of the sum of contributions: 
from many component dead loads. According to the law of large nu: numbers, 
the variability in the total dead load will be less ‘than the 10% value sug suggested - 
for the design of the superstructure. > example, the COV of the total load — 

Y from n floors can be shown to be (e. see Ref. 1). 


in 1 which 8, =the COV o of | the | load from inc individual fl floor; and p = - the cor correlation. 

coefficient between any two floors. Suppose p = 0.3 and mn = 20, then 8, 
i only 58% of the value of 5,. -In fact, ‘in n those foundations where dead load 


as in the evaluation of overburden pressure in settlement anabyel. The vertical 


and lateral earth pressures will be treated separately water soil induced 


Most st loads acting on a foundation do vary with | time. The simplest edits : 

_ is the sustained live load, caused by the weight from regular human occupancy ; 
(as opposed to “occasional large gatherings which will be treated later as San ; 
extraordinary load), furniture, and equipment. The regular occupancy load is 
not expected to vary much throughout the duration of a given occupancy. Such 
load, however, could change significantly between occupancies. For practical os 
purposes, the regular occupancy load could be assumed to ) stay ¢ constant during 
each occupancy. An example of the changes in regular occupancy load over 
the life of a structure is shown in Fig. 1. The times at which occupancy changes © 
do not necessarily take place at regular intervals. The Poisson process has 

@ suggested to model the occurrences of such occupancy changes, with 

_ parameter v as the average rate at which occupancy changes. This in turn implies — 

that the duration of each occupancy follows the exponential | distribution. 

_ mean number of load changes over a structure’s life of T years is simply vT. 
as bearing capacity) and settlement sand, 


| 
Za that the variability of the resultant dead load is expected to decrease 
| 
| 


| 
in which Fy the CDF of the load each regular occupancy; i= 
_ the number of changes of occupancy; and the values of X, between occupancies 
‘ are assumed to be identically and independently distributed. If exactly n 7 
of ee are gourd to have occurred, the CDF of X.,,,, is given nby 


‘The Statistics of regular occupancy load X, have been collected by ore a 


ma baer auapt sey in the United Kingdom and by Culver © in . the United 


FIG. 1.—Changes of Regular Load 


fier area of 2, 000 sq 1 ft (186. lm m?), the COV of X,is 0.5. 
= of the lifetime maximum regular occupancy load could be evaluated for & a. | given 
Ms distribution of X, and a given expected number of occupancy changes vT. 
Assuming X, to be lognormally distributed, and for eight changes of occupancy, 7 
7 the mean and cov of X,,, for a floor 2 area of 2,000 sq ft (186. 1 m”) could os 
4 m’) and 0.30, respectively. 
The total ‘regular “occupancy loads. from the individual floors are required 
_ forthe design of foundations supporting a multistory structure. Since the maximum 
_ ‘Regular occupancy loads of each floor may not necessarily occur at the same 
- time, | it would be somewhat conservative if the maximum regular « occupancy 
loads i in each floor are combined to give the maximum total regular « occupancy 
load for foundation design. Since the times at which changes of « occupancy > 
pret are random, the total regular occupancy load at the foundations will 
. appear as shown in Fig. 2, where load changes occur much more frequently 
si especially where there are many floors. The load process now is correlated 
al since at given load change, the load at only one floor changes whereas the 
Toads i in the other floors remain as before. Wen (34) has studied the Statistics 
of the maximum of this correlated load process over the lifeti 


_ foundation performance is governed by the maximum load during the lifetime Bie : 

of ‘the: foundation. T herefore, the maximum regular occupancy load over the : 

. - foundation’s lifetime T is of interest here. The cumulative distribution function 7 

(CDF) of the maximum regular occupancy load X,,,can be derived from 

aq 
the standard deviation of X, will decrease as floor area increases. For a typical _ 

| 
ail 


= 


For a 50-story structure with vT = 7.2 for each floor, the mean i tae maximum — “ 
lifetime regular occupancy load is about 23 psf, compared to 20.5 psf for the © 
- mean maximum of a single floor. A correlation coefficient of 0.3 was assumed 
as suggested in Ref. 19 to account for the dependency of loads between successive ds 
floors in the foregoing analysis due to possible common usage. 
SS On the other hand, if settlement due to consolidation is to be evaluated for ff 
a foundation, the loads considered would not be the lifetime maximum, which = 


 walgined average of the ¢ regular occu occupancy load will be needed. T he « 


in which X,, = the ith regular occupancy otiaa ond a,= 
factor depending on the duration of X , and the time of occurrence of X,, 


relative to the entire consolidation process; and is the summation from 


= ae n where n = the number of ‘occupancy changes occurring | over the 


cov of the effective load could tobe 


in which = the COV of the regular load to the 
a same for each load period and statistical independence is assumed between 
i - each value of X,,. The coefficient B can be shown to be less than unity and ad 
bi its value will depend on the number of changes of occupancy n. As an example, 
for n = 3, B = 0.6 (if the ratios of the values of a, are not larger than say, a 
_ 4) indicating that the variability i in the effective load for settlement ca calculation is ; 


will be smaller than that of each regular occupancy load. In fact, as the time 

required for consolidation increases, n will be larger and B will decrease - 
‘i accordingly, thus the variability in the effective load in a consolidation settlement ae 


analysis would be relatively small compared to the variability i in the soil parameters =. - 


| 
| 
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MG PROBABILISTIC EVALUATION 
a heiite classified as extraordinary loads include forces due to snow, wind, 
and storm waves, earthquake and extreme occupancy loads. Ths pao charac- 


a _ the time, such loads do not occur at all (or normally their tert are negligible 

"such as light breezes or ordinary waves). When such loads do occur, the duration 

_ as well as the load magnitude are uncertain. In view of the short duration of | 

_ load application, relative to dead load or regular ‘occupancy loads covered 
previously, the contribution of such loads to consolidation settlement will be 
negligible. ‘The extraordinary loads are, , however, critical for problems related | 
to § stability and transient 1 motions. Some common features 0 of and statistical ‘2 
concepts relevant design against extraordinary loads are “considered in the 
Lifetime Maximum Extraordinary Load. —Assuming the occurrence of 
extraordinary load follows a Poisson process with mean rate of occurrence 

v, and the magnitudes of the _ between occurrences are statistically indepen- 


ab 


Wind 


- dent, then the CDF of the lifetime maximum of an extraordinary load rs on 
_ be shown to be the same expression as that in Eq. 2, in which X, is interpreted 

_ as the maximum load magnitude X in each occurrence. "Alternatively, if the 


- statistics of an extraordinary load, such as wind speed and wave height, are 7 


—Extraordinary Loads 


_ available in terms of the annual maximum values, then the CDF of the lifetime 
-maximum load is the same expression as that in Eq. 3 in which X, = 
_ is the random annual maximum load assumed to be statistically edepentet 
; between years; and mn = T is the foundation lifetime. Once the CDF has been 
evaluated for | a of x or X,, the mean, variance, and COV 


dingly. 
‘When vT or T is large, could be shown to approach o one of the limiting 
distribution, such as the type I extreme value (Gumbel) distribution or the type a 
II extreme value (Fisher-Tippet) distribution or the Weibull distribution, = = =«_—> 
Return Period of Design.—Conventionally, a 100-yr wave, or a 50-yr wind 


velocity or an earthquake acceleration with a return period of 500 yr may be be 


q 
4 
| 
| 
 extraor inary 10a that wou appue to a foundation system. ese esign 


wes’. 


The return period T is simply the average time until the nares —_ is exceeded. 
It could be shown that thy) cme al 


in which = the | probability that ‘the “design load x exceeded in a given od 

year. For example, if 50 yr is the design return period of an extraordinary — 
load, the probability that such design load will be exceeded in a year is 1/50 
or 0.02. Therefore, the probability that the , design load will be exceeded during wp 
the foundation system’s expected lifetime, say of 30 yr, could be evaluated i; 

_ using a binomial distribution as 1 — (1 — 0. 02)” or 0.45. It is interesting to 

observe that the probability of exceedence in 50 yr (i.e., within the tales 2 
return period) is 0.636, implying that the design load in this case is not necessarily _ 
exceeded at the 50th year or once every 50 yr. In fact the design load could 

be exceeded in 7) year including e even the first few years of operation,  =—> 


FIG. 4- 4.—Earthquake Acceleration Time History 


Dynamic Amplification Factor.—A inspection of the characteristics 
of an extraordinary load when it occurs may reveal fluctuations as shown in — Z | 
“Fig. 4 (for the case of earthquake acceleration). The average period between - 
‘peaks of such fluctuating load cycles could range from a fraction of second 4 


(as in earthquake excitation) to a fraction of a minute (as in ocean waves) = 
or perhaps a longer period (as in snow accumulation). Depending on the vibrational rel 
characteristics of a structure-foundation system relative to the frequency of a 


fluctuations of the applied loads, the dynamic system response in terms of 


= . Pseudostatic analysis: An example of this is the seismic coefficient adopted — 
in quasistatic analysis of embankment “stability, which could be a \ prescribed — 
_ value or determined through a dynamic visco-elastic response analysis (e.g., 
see Ref. 25). Wu and Kraft (35) suggested the COV of the prescribed seismic 
coefficient could exceed 20%. The uncertainty of the seismic coefficient deter-_ 
mined from the latter — is is expected to be smaller, but it should incorporate 


— 


_ | 


in parameters s such as soil damping, mass shear wave ve 
velocity, and most of all in the visco-elastic model itself for representing the — 
_ 2. Time History Approach: The effect of the uncertainty in the load time 
history input could be analyzed | by ‘Tepeating the dynamic analysis to various 


‘to the s same maximum  sccstorations. Moreover, using a Monte-Carlo simulation - 4 
‘procedure, uncertainties in the parameters of the deterministic model such as rs a 
_ the stress-strain relationship, damping characteristics, and density of the soil 
= might also be incorporated in assessing the overall variability in 
a dynamic stress. Tremendous amount of computer time is generally required. 
_ 3. Random Vibration Analysis: Christian (4) summarized the concept of random — 
_ vibrations as well as the state of the art of its application to oo 
engineering. Through random vibration analysis the statistics of the Guaie 
stress response could be derived imanowins oft 
Nevertheless, the statistics of the stress would account for only 
fk uncertainties in the inherent variability of the time sequence of load fluctuations. a 
_ Uncertainties in the system parameters would also affect the dynamic stress. - 
Since the dynamic response could be analytically expressed in terms of the | 
system parameters, first-order analysis of uncertainties (based on the Taylor 
Series approximation) might be applied to incorporate the system uncertainties 
_ in evaluating the overall uncertainties of the dynamic stress response. An example — : 
Of this approach to wind loadings could be found in Ref. 22. a nua mu te oi} to 
In the following sections, each of the extraordinary loads that would be of 
j ; interest to geotechnical engineers, is examined with respect to its sources of 
uncertainties and probability models that have been suggested. 
a _ Wind Load.—Wind load generally acts laterally on a structure, thus creating _ 
: an overturning moment on the foundation. As one side of the foundation receives 
downward load from the overturning moment, “uplift forces will « occur on the © 
other part of the foundation. Pulling of piles and caissons could be potential 


modes of failure in addition to stability failures under wind action. Moreover, — 


if the wind velocity is much larger in one direction, say from the west as 7 

: compared to thet from the east, then differential stress may be induced on 7 
_ two opposite ends of the foundation, which may lead to differential settlement 


or tilt if the wind environment persists over a long time period. 


— 
_ The annual maximum wind velocity has been suggested to follow the extreme 


‘should be used throughout the United States. However, it has also been shown 
(27) that in a well-behaved wind climate, such as the continental interior which 
, is unaffected by hurricane disturbances as in the coastal areas, type I extreme 


value distribution. Thom (30) proposed that the type II extreme value distributions 


value distribution may be more appropriate. In areas where tornadoes and = | 
4 = are prevalent, though at infrequent intervals, the effects of these 1 
extreme winds should be superimposed on the ordinary wind effects. The } 
distribution of the combined maximum wind velocity (ordinary wind plus 
hurricanes and tornadoes) has been derived (32), which could deviate 
pressure is proportional to the square of the wind velocity. 
-— mathematica simplification could be achieved by adopting the type II instead i : 


i 
- 
i 
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‘of ome I distribution for on wind velocity, since the corresponding v wind pressure > 
can be derived to follow also the type II distribution. The typical COV of | 
the maximum mean hourly wind velocity would range from 0.09-0.12. The cov = 
of the corresponding maximum wind pressure will be approximately twice esl 
pa V’), thus ranging from 0.18-0.24, with the lower value for the suburban 
area and the higher value for the inner city. These statistics, however, assume — 
only static wind load. For a structure-foundation system with large natural period 
of vibration, static wind load may be satisfactory; otherwise, the dynamic response 
"4 of the system to wind forces should be considered. Combining the pebity 
distribution of wind velocity and direction and the response characteristics 
obtained from wind tunnel tests, Davenport (7) established probability contours : 
of specific response (stress or displacement). Nevertheless, uncertainties in the _ 
response and structural parameters were not considered in his evaluation. If — 
the effects of the uncertainties in the gust response factor and the turbulence _ 
nature of the wind environment are incorporated, the value of the COV of 
wind forces could range from 0. 27-0.52 for city environment 0.25-0.48 
for suburban areas (22). An independent evaluation (21) suggests s that the COV 
of the lifetime (50 yr) maximum wind pressure would typically range from 
7 0.37-0.41 (assuming a specific COV of 16% for the lifetime maximum wind : 
which appears to agree well with the previous study. 
Earthquake Load.— Earthquakes have been regarded as a major threat to a 
4 dam safety. For slope stability consideration the maximum lateral force due _ 
to earthquake motion could be modelled quasistatically as a function of mass _ 
of the structure, the maximum ground horizontal acceleration during earthquake = 
excitations, and the dynamic amplification factor. Several factors need to be t 
_ considered in evaluating the statistics of the maximum ground acceleration; — 
these include: (1) Identification of all the potential sources of earthquakes (such 
aS existing faults) in the vicinity of the site; (2) the probability distribution 
(Of earthquake magnitude; (3) frequency of occurrence of earthquakes; (4) the | 
direction of fault slip, if any; (5) the distance from each source of earthquake; | 
- and (6) attenuation characteristics of the region. Analytical models for estimating — 
_ the maximum acceleration at a site have been proposed by Cornell (5), Der- 
Kiureghian and Ang (9), considering separately different types of a 
sources (i.e., fault lines, aerial sources and point sources), each having the 
magnitude distribution and activity rate assessed based on available earthquake © 
7 q records. The relative frequency of earthquake magnitude in a region is assumed 
to follow the Richter’s Law. They also used attenuation equations to relate 
site ground acceleration intensity t to © earthquake magnitude and focal distance. ‘ 


a concentrated at a point whereas DerKiureghian and Ang’ s s model assumed ‘hat 
_ the slip length for a large earthquake may be several hundred kilometers, and 
thus, the largest contribution to the acceleration intensity at the site might be — 
caused by the slip location closest to the site. Although the latter model appears — 
_ to be more realistic, difficulty exists in defining the magnitude of an ome : 
in such a case. Because of lack of data for establishing the attenuation equation J 
for sites close to the sources, prediction of acceleration level is expected to 
be very crude. The maximum ground acceleration for San Francisco and San 
Juan, Puerto Rico have been evaluated (8); the results are reproduced in Fig. 

5. The The probability that a @ ground of 0.4 g will be exceeded ina 


: : 
° 
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q is 10 times in a San site ‘than i in a Based 

i on this annual maximum acceleration distribution, the distribution of the lifetime pat yg 

J 
acceleration could be derived accordingly. To convert the ground 
acceleration to the maximum acceleration in the structure-foundation “system, a 
__ the response characteristics of the system should also be considered. Portillo & 
& (20) estimated the COV in the dynamic amplification factor is typically around | ; 
0.33 for reinforced concrete structures. For soil material whose dynamic stress- _ 

a _ strain relationships are less well understood and where soil parameters are subject 7 
to spatial variation and sampling uncertainties, the COV of this factor is expected _ 
to be somewhat higher. Nevertheless, since the level of uncertainty associated © 
_ with the attenuation «avi ranges from 0. 5-0.85 in terms of COV (10,11 17), 


a 


acceleration would perhaps come from the ‘seismicity ‘factors than fr 
Structural systems response. 
In those soil conditions where substantial strength loss could result from 
pore pressure development during earthquake excitation, liquefaction analysis 
_ is required. For the evaluation of liquefaction potential, the number of significant ra q 


i equivalent stress cycles during earthquake excitation is also important (26). The 
= of equivalent cycles, N.,, depends on: (1) Duration of an earthquake; : 


(2) magnitude of an earthquake; “G) equivalent stress level expressed in terms 
of percentage of the peak stress; and (4) the soil strength curve at the specific 


| 
=, 
| 
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. site characterizing the red. reduction in 1 soil strength . with increased number of constant 
- load cycles. Based on the data in Ref. 16, the relationship between N ., earthquake 


magnitude, equivalent stress levels, and soil strength curves has been enslyned 


statistically (13). _ Nonlinear regression analyses: were performed to obtain the 
mean and variance of N,, as a function of magnitude M. The sensitivity of a. 
the prediction equation is found to be minimum when the equivalent stress 
level is adopted as 75% of the peak stress. The corresponding mean value 


and the scatter about this prediction curve is represented by a constant standard — 
= deviation of 5.4 cycles. Combining Eq. 7 with the probability distributions of 
the annual maximum earthquake magnitude and acceleration, the annual —- 


ity of the liquefaction “ a given site could be estimated (14) — 


in which = the of liquefaction at t 


j site assuming the load parameters, namely the value of the maximum acceleration , ; 
and the number of ‘equivalent cycles, , are exactly known; 


f = the joint probability density function (PDF) of and 
Ne on The liquef: action risk over a given expected lifetime could be then evaluated 
fi rom the annual risk. The foregoing approach determines the risk of liquefaction 
at a single point (which is perhaps the weakest) in a soil stratum. However, 
7 large volume of sand needs to liquify to have an adverse effect on a structure — 
a the ground surface. Zones of liquefaction ar are nia suman in a recent probabilis- 
Wave Load.—With the increase of offshore activities, geotechnical engineers 
- confronted with a new kind of load and environmental factors. The forces 
from ocean waves have been recognized as a major contribution to the loads 
oe the stability and displacement of offshore platforms. Although eco- 
nomic lives of most offshore structures are expected to be | approx 20 yr, the — y 
100-yr wave has been recommended as a characteristic load for ‘subsequent 5 
deterministic analysis. Realistically, the maximum wave height during a oa 
E period is a random variable. Based on a probabilistic stability analysis for a 
_ typical Condeep gravity platform located in the North Sea (15), and neglecting — 
| Caadiiaten effect on the soil from the platform weight and previous storms, — 
the probability of stability failure assuming a random 20-yr maximum wave = 
. ee: is approximately equal to the probability of stability failure assuming ‘ 
a deterministic 100-yr wave height. Thus, the 100-yr design wave height may 
perhaps be a consistently (reliability-wise) chosen characteristic load for stability 
= The long-term wave height distribution has been modelled by the > Weibull 
aa (2) or alternatively by the Gumbel distribution (23). The general 
lack of sufficient field measured data and possible errors associated with the F 
_ eememes distribution model will give rise to uncertainties in addition to the 


inherent variability. Calibrated wave and those hindcasted from 


1 i 


storms are often used as supplemental data. In th this case, calibration and | ca 


error would need be considered. ‘For a wave height and 


or p leherstery model tests or sophisticated analytical studies. Varying levels of 
are associated with each of these example of the i 


es North Sea based on 4 yr of wave height data, is summarized in Table 1. “ 
. sensitivity analysis of the laboratory model test data reveals that the estimated — 
uncertainties in wave height and period would result in a COV of 11% in the 

_ wave forces; whereas discrepancy between laboratory model tests and field 

_ performance could be subject to a COV of 20% based on the assumptions 

_ that the uncertainty level should be somewhat less than that observed | for smaller 

structural members i in steel jacket structures (30%-40%). 
‘The undrained shear strength of a soil c could be reduced 


* Baer wave action of the storm prior to the occurrence of the largest - 


of Uncertainty Analysis of Wave Load for Platform 
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Wave height: insufficient data bs 


Wave height: distribution model error 
‘Wave height: Total UF of 


Wave forces: due to uncertainties in wave height and period (heal evil 


wave. The reduction in strength is related to the accumulated cyclic strain which | 
_ in turn depends on the height of the individual waves. The distribution of wave — 
heights in short-term sea state (such as during a storm) has been suggested — 
to f follow a a Rayleigh distribution, from which the | e statistics of the reduction — 
in undrained strength could be evaluated. Toe 
In problems where the dynamic system response is critical, wave , energy 
"spectra such as Jonswap or modified Pierson- Moskowitz-Kitaigorodskii (P-M-K) © 
have been suggested to describe the stochastic nature of waves. Using these _ 
7 spectra as input, random vibration analysis may be performed to obtain statistics 
of the dynamic response. The problem of which spectrum ‘to is still 
— Geaeeenne On the basis of 41 wave spectra from six storm records, Sartre 


(23) suggested that the ahongype peak of the samsttiod spectrum appears to fit 


defining the spectrum for ‘the ‘severe storms, €.8.5 that v with a 20-yr return period ( 


y 
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. _ Each load component considered previously acts separately on a foundation, 
whic occasionally « could be subject to stresses resulting from combined actions na 
of two or more load components. The critical load combinations that a foundation ‘a 
may encounter sometime during its lifetime are identified as follows, together 
with a brief review of the corresponding models for probabilistic treatment. _ 
- Dead Load + - Maximum Regular Live Load.—Since dead load is assumed 
to act permanently, and regular live load only changes with each turnover of a 
occupancy, it is obvious that the combined action of these two loads will achieve 
a maximum during the occupancy with the largest live load. Let the maximum a. 


which X,,, = the maximum regular occupancy leads and Y= = the & 
tatic load. "Assuming X,, and Y to be statistically independent, the CDF of | d = 
could be shown (e.g.,see Ref. l)tobe 


i in which | le may be determined from Eqs. 2 or 3; ; and fy) = 

Alternatively, if only the mean and standard deviation of Z are required, 

they could be evaluated directly from Eq. 9 in terms of those of X,, and : 


= « 


Dead Load + Regular Live Load + 1 Extraordinary Load. —When one of «ff 

_ the extraordinary load (e.g., wind, wave, snow, earthquake, or cnsecotinnsy 

~ live load) achieves its maximum value over the lifetime of the foundations, 

it is unlikely that it occurs during the occupancy with the maximum regular — 

live load. / Assuming that there is only one extraordinary ¢ and live load over | 

the period of one occupancy could be first established as follows. Let 7 


in which X = the regular live load ‘during a Sada ounminan: and E= the 


| maximum value of an ‘extraordinary load over the period of the given occupancy. _ 
Similar to the tog case, the CDF and statistics of Z could be shown to © 


the: 


be evaluated according to the procedures covered in previous sections. Assuming — 

Poisson process for the occurrence of occupancy changes, Eq. 2 can be applied 
a to saath the CDF of the the maximum c - extraordinary | and live load over 


a 


in which v = the mean rate of ‘occupancy change; end T = the lifetime of 
the foundations. The maximum combined load over the foundation’s lifetime — 
ist then the sum of Z,,,, and the permanent static load Y. Replacing X,, by 
TOYA Eqs. 10-12 could be then used for determining the corresponding CDF 
Dead Load + Regular Live Load + 2 Extraordinary Loads. —Let R» denote 
the maximum load value due to two independent extraordinary load sources 
over the foundation’s lifetime. The value of R will be higher than that of E 
_ in the previous case for two reasons. First, the second extraordinary load itself Pe 
could have a maximum lifetime load higher than E. Secondly, the two extraordi- 
_ ‘fary loads could occur + simultaneously and the « combined load could exceed 
_E. The probability of simultaneous occurrence of two extraordinary loads will 
_ increase with the frequency of occurrences of each load as well as the duration — 7 
of each load . Suppose the loads follow Poisson processes with mean rates of 
occurrence v, and V2, and the mean durations of each extraordinary load when 
it occurs are Ma, and respectively. Wen (33) suggests an approximate 
xpression for the CDF of over a time t as 


from j 
in which F,, and Fea = 


for evaluation ‘of the CDF of the lifetime n maximum load. 
. The aforementioned formulation could be extended to the case cm more than 


4 


> encountered by geotechnical « engineers is that 
! by the soil itself, such as the overburden pressure on a given soil layer in | 
a settlement analysis. The average density over a soil stratum would be the a 
relevant parameter here instead of the density at a Specific hiyoms Because - 


a - of reduction in the variability depends on the volume of soil under consideration 

_ and the spatial correlation characteristics of the soil density. In general, for | 

: t natural deposits or compacted embankment, a much larger correli ition is observed 
in in the horizontal plane than along the vertical direction. Assuming perfect q 
‘correlation | among soil density in the horizontal plane and an exponential decay 7 
“Sensis ov model vertically with correlation distance b, the COV of the — 


Ba over a soil layer of thickness H can be derived $F ainvisia, 


the S Of extraordinary loads £,, £,, an 
- (31) suggested that the variability in the average property (e.g., density) over § 


rrr 8, = the COV of densiiien between the individual soil | samples. For 
Sale if 5, = 10% and b= 0.3 m (0.98 ft), 8, for a l-m (3.28- -ft) and 
10-m (32.8-ft) thick soil layer could be shown to be 6. 5% and 2.5%, respectively. 7 
4 Therefore, the variability in the average density over a large stratum of soil 
bi g expected to be very small, except for the case where the soil stratum is 
— Lateral earth pressure, encountered in retaining structures, are aie 
: ‘empirically calculated in terms of the average density ) and a stress coefficient _ 
such as K,, K,, or K,, which is in turn a function of the friction angle 4. ¥ 
The COV of could reach 10%-20% (28). The uncertainty in the coefficient a 
oe K due to the variability in could be evaluated using the first-order approximation | = 
formula. For example, with K, = (1 sin o)/(1 + sin >), the mean | and be 
oy 


— 


6.—Histogram for Earth on Braced Cut 


a) Dre; } 
which yields of K. 3)8, after simplification, an illustration, 

_ for a design situation where the mean friction angle is 30° and a COV of 10%, 
_ the corresponding COV of the active pressure coefficient K, is approx 
Additional error could also exist in the empirical model itself, assumptions adopted — 
in the Static of residual s stresses and | values of wall friction. 


34 
: ‘amean of 0.4 anda COV of 0.23. After removing the contrib 
iz 


known as negative skin friction) could be generated on the piles or omescn 
z= to relative downward movement of the soil adjacent to the pile surfaces. 

_ force is function of soil parameters such as friction angle and cohesion — 
coefficient. ‘Again, spatial variability, measurement, and sampling error cece 7 

 ealian ‘and model i uncertainties would give rise to variability in the drag 

‘a — which could be conveniently analyzed through first-order approximate 

Pore pressure is another major parameter in geotechnical analysis. Even in 

: simplest case where pore pressure results from hydrostatic condition defined © i 
by the ground-water level, uncertainties in pore pressure at a point in the stratum 
‘still arises due to seasonal changes, rainstorms, construction, or dewatering 
activities in the adjacent neighborhood. Spatial fluctuation in pore pressure could 
also occur in addition to temporal fluctuation. However, foundation performance __ 
is generally governed by the average pore pressure over an area. Thus, the 
effect of spatial fluctuation will not be as crucial as the temporal fluctuation. 
If piezometric measurements are available, long-term predictions could be made 
from the recorded seasonal cycles and short-term stochastic fluctuations (f rly 


rainstorms). Alternatively, if piezometric data are not available, which is generally 

the case, , hydrologic m odels considering precipitations and ‘infiltration : rates may . 
be used to predict ground- -water levels. From the annual maximum water level, 7 
statistics of the lifetime maximum ground-water level could be evaluated sol 
procedures suggested earlier. According to the measurements recorded by the | 
Illinois State Water Survey, due to natural causes, the annual maximum ground — 

water is typically 5 ft 53 m) above the mean ground- “water level with a 

of 35%. Over a period of 20 yr, the | high-water level could be expected 

af be more than 8 ft (2. 44 | m) above the mean ground-water level. On the other 


hand, if low ground-water level is of interest, the CDF of its lifetime minimum > 


Frigg = 1 = 
in which X, = the minimum ground-water level, 
n ~ Movement of water through soil stratum would give rise to seepage forces” 
which is a function of the hydraulic gradient which in turn depends on the 
7 relative permeability over the soil stratum. It is generally recognized that - 
error in measuring permeability could be plus or minus one order of magnitude. * 


The ratio of horizontal and vertical permeability, that is important in the 


‘~ construction of flow nets, could be thus subject to large uncertainty. In addition, 
: due to inhomogeneous deposits of sand and clay, and erratic location of sand 
lenses, further uncertainties are associated with the prediction of seepage force — 


; In the case of slopes, based on Refs. 3 and 24, the COV of the pore ¥ pesmuee 


was estimated (36) to vary, ‘between 0.05 and 0.1. With various assumption 


| from K, (estimated at 12%) the uncertainty in the model is estimated to be te 
i about 0.19. Nevertheless, data are still generally lacking for evaluating such > ; 
a hard lave Overiain Dv 4 consolidating atum force pm mon 
Ss q 


_ flow conditions, they also estimated that the COV of the pore pressure ‘gia : 
the sliding surface for a slope at Selset to be 6%, in addition to the uncertainty — oo 


with adequate statistical data, the probability of encountering a certain magnitude - 
of load or load combination over a foundation’s lifetime could be evaluated. rs 
Such information is necessary in rational planning of geotechnical system and q 
the establishment of load values for design purposes. ~~ oe ee. 
The ounia reported here was supported in “part by the ‘National Science _ 
 Senhaien (NSF) under grant No. PFR 78-19680. Any opinions, findings, = 
1 conclusions or recommendations expressed herein are those of the writer and 
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F(x) cumulative distribution function; ol 
thickness of soillayer; 


number of floors or occupancy changes; 

probability of exceedence, wind pressure; 

maximum of two extraordinary loads combined; Ta 


T = expected service life; 


coefficient; 


density of soil; watt W 
mean rate of occurrence; 
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Metastable o or r collapsible soils are defined as any ‘unsaturated soil that 
through a radical rearrangement of particles and great loss of volume upon > 

wetting with or without additional loading (2). The existence of these soils —~ 
throughout the world and difficulties with building on them have long been 
_ recognized. The reason for the lack of information on many y of these soil deposits . 
is that they are located in predominantly arid regions with limited economic — 
development. Recent advances in irrigation have made it possible to open up — 
many of these regions to farming and industrial development. These developments . 

include industrial and urban complexes, and provide opportunities for the use - 
Of water in 1 large quantities. The consequent problems resulting from excessive ae 


sive state-of-the-art ; papers on n collapsing ‘soils a are by Northey (18), Sultan nl 
| and Dudley (9). Bara (2), in 1976, has also presented a brief summary of present 
theories and treatment of collapsible soils. Since 1970, the major research has 
. been devoted to determining the mechanism of collapse. Other studies have 
been devoted to predictive methods, treatment methods, and case histories 
(3,13,21). The most recent study is a comprehensive review of poe orn 

if ‘ investigations of loess in North America by Lutenegger et al. (16). _ oa 
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q The most extensive deposits of collapsible soils are Aeolian or or deposited 
san 


ds and silts (loess). In addition, alluvial flood plains, fans, mud flows, colluvial 
deposits, residual soils, and volcanic tuffs may produce collapsible soils. In 
most cases, the deposits are characterized by loose structures of bulky-shaped Pp 


UUIA 


- grains often in the form of silt to fine sand size ———_ 
_ Aeolian Deposits.—These deposits consist of materials transported by 
which form dunes, loess, loessial-type deposits, aeolic beaches, and large volcanic 
dust deposits. They consist of cohesioniess or slightly cohesive soils and may | 
have a low relative density. The natural structure of these soils may contain 2 
a clay cement binder such as the loessial soils. These deposits are characteristic = 
of arid regions where the water table is at great depth below the ground surface. ~ 


Upon saturation, the clay cement binder loses its strength. Under these conditions, 
the material suffers structural collapse. The fine aeolian sediments are sometimes — 
found protected from weathering in nature ss de a ‘surface crust of impermeable ; 


— 


to 
_ Clay. This crust protects the entrance of large quantities of water, thus preserving 
the natural structure. For this reason, even in climates of medium rainfall, 
_ fine aeolian deposits are modified either 1 not at all or only partially. HOS aes 
: _ Large areas of the earth’s surface, particularly in the midwestern and western = 
- United States, and parts of Asia and southern Africa, are covered with loess 
that is susceptible to a large decrease in bulk volume when saturated. Loess, 

_ a yellow to reddish-brown soil covers approx 17% of the United States, — : 
17% of Europe (including the low countries), parts of France, Germany, and 7 
J Eastern Europe, 15% of Russia and Siberia, and large areas of China. Loess “ 

is found in New Zealand and the plains region of Argentina and Uruguay, | 
principally between latitudes 30° and 40°. Loess appears to be ened by 
wind-borne traveling over glacial outwash causing of of the 


Water- Laid laid deposits ‘consist primarily of loose water. 


| fi 

wos 

a 
CLAY BOND TYPE OF STRUCTURE 

| 
Oaquin Valuey Oo auiornia, €.g., Cases have been oDserved Of subsidences 


: 
up to 15 ft. in water-laid deposits. These were originally deposited 
7 by flash floods or mud flows derived from small watersheds subject to cloudbursts re 
at infrequent intervals. These deposits dry out and never again become saturated _ 
tae until the arrival of another flow. Flows consist of poorly consolidated materials a a 
4 that contain considerable. clay. The amount of clay in a flow has an ‘important — 
bearing on the soils’ behavior. Bull (6) has shown that the maximum subsidence a 
"occurs when the clay amounts to about 12% of the solids. Below 5% there : 
is little subsidence, and above 30%, the clay swells. In between, there are 
- cases in which soils with larger pore spaces settled more than those at the 
same void ratio, but with smaller pores. ight 
Residual Soil.. —Residual soils are the product of weathering, i.e., the disin- 
- tegration and “mechanical alteration of the ‘components < of parent rocks. The 
particles of residual materials may vary in size from large fragments to gravel, — 
a ~ sand, silt, colloids, and in some cases, organic matter. The collapsible grain 
Le "structure has developed as a result of leaching of soluble and colloidal material. 
a This leaching out of the soluble and fine material results in a high void ratio 
and unstable structure. Low densities may be found in the upper part of the 
an due to eluviation. Generally, residual soil profiles are encountered in 
q thicknesses ranging from a few inches to several feet, depending on the dimese 
3 


and physiographical environment of the region. Residual decomposed granites “ 
in South Africa and northern Rhodesia, e.g., have been found tol have a total 
collapse amounting to about a 7%-10% increase | in density (5). ef fiok 

- Other soil types that exhibit collapse uf upon “wetting are those derived from > ¥ 
volcanic tuff (17), gypsum (25), loose sands cemented by soluble salts, dispersive _ 


4 - Russian studies (11) have shown that there are four main types of ‘weting 
x that can trigger the collapse of soils; 


1. Local, shallow wetting of a random nature caused by water sources e 


on 


“J 


_ pipelines or to uncontrolled drainage of surface water from construction. Water 
from local wetting does not usually penetrate to great depths and there is normally — 
th rise in the ground-water level. In this case, settlement occurs mainly in 
“the upper layer of soil below the wetted zone. suugradou 
ie Intense, deep, local wetting of soil caused by discharge of industrial effluents _ 
or irrigation. If the flow rate is sufficient to cause a continuous rise in the _ 
ground-water level, then the entire zone of a soil may be saturated x 


is extremely 1 uneven n and ‘dangerous. ‘It may include the | ‘whole thickness of e 
the soil layer and collapse may occur either under the load from existing —. : 
orunderthe weight of the soilitself. 
3. Slow and relatively uniform rise of the ground-water level under the influence Pa - 
of water sources outside the collapsible soil area. The settlement is normally — 
uniform and gradual. aa) larg ¢ on) dowe 

‘ 4, Gradual slow increase of the m moisture content of a thick layer of — 

za resulting from condensation of steam and accumulation of moisture caused 7. a, 


1 
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partial. Correspondingly, the settlement i is incomplete, and slowly ir increases 


_ with the wetness of the soil. Collapse may be triggered by water alone, or 
by saturation and loading acting together gon!) 
According to Barden et al. (3), appreciable collapse of a soil Tequ.res the 
following three conditions: (1) An open, potentially unstable, partly saturated — 
structure; (2) a high enough value of an applied or existing stress component _ 
_ to develop a metastable condition; and @) a strong soil bonding or cementing 
to stabilize intergranular contacts, with a a reduction which, we tting, 
All cases so far studied have shown that these soils have a hartge structure 
of bulky-shaped grains with the grains held in place by some bonding material 
or force. When the support is removed, the grains are able to slide over onto > 
one another, moving into vacant spaces. 
In cases where the soil consists of sand with a fine silt binder, the temporary _ 
"strength is due to capillary tension. As the soil dries below the shrinkage limit, — 
_ capillary action places the water under tension. This increases the apparent | 
_ Strength of the soil. However, addition of water reduces the strength, and if 
_ the soil is porous, a rapid decrease in volume may occur. According to Ancecs 
and Donald (1), for uniform spherical grains in an open or cubical packing, ; 
the maximum added pressure due to the capillary | effect occurs at a moisture — 
- content of about 32%. For the densest packing of uniform spherical grains, © 
_ the maximum tension occurs at 10% moisture. Investigations conducted with 
: various collapsible soils indicate that the peak effective stress values occur 
at moisture contents less than saturation and above 10% water content. i 
Knight (15) has shown that sands and silts loaded | in a the unsaturated state 


4 line. It has been suggested by Blight 4) that the \ uniqueness of the virgin \ saturated | 
; : compression line holds for swelling soils, although Burland (7) has expressed | 
_ the view that the soils with high clay content, the effect of stress history on 
\ soil structure will become significant. Holtz and Hilf (12) explained the collapse — 
on the basis that saturation releases confining stresses, and the sample fails 
in shear, thus undergoing additional settlement. Holtz and Hilf (12) also obtained 
_ compression curves showing that collapse occurs in soils compacted to high 
_ densities when dry of optimum moisture content. 
_ The majority of collapsing soils involve the action of clay particles in the 
J bonds between the bulky sand grains. A number of structural arrangements 
of the clay Particles a are possible depending on the history of the soil (see 
Fig. 1). Clays that are formed by authigenesis possibly form a parallel plate” 
 onionskin effect around the quartz particles. Gradual evaporation of the pore X 
‘water can cause the clay plates to retreat with the water into the menisci at 7 
interparticle contacts. Knight (15), using an electron microscope, found that a 
_ under such conditions the clay grains cluster around the junctions in a random 7 
flocculated arrangement, giving a buttress type of support to the bulky grains. _ 
The nature of the clay bonding is complex and it is not yet clear how much © 
is due to how much to capillary effects. However, 


| 


209 
Ss it has been observed that in most om, the ae the water content, , the greater 
Cementing agents as iron oxide, carbonate, oF the welding of 
_ the grain contacts, provide ‘Strength for many collapsible soils. This comenting — 
action is often the main agent in loessial soils. The rate at which a bre 
agent loses its effectiveness depends on the degree of contamination of the 


Soils have 


“Soils have not 
generally 
observed to 
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An increase in load could add to this effect, as an increase pressure on 
‘the soils will increase their rate of solution, which could produce a delayed ca 
_ increase in consolidation. However, whatever the physical basis of the bond 
- strength, all collapsible soils are weakened by the addition of water. Collapse 


= strength is more immediate in the case where the grains are held together 
by capillary suction, slow in the of and much slower 


in the case of clay b 

ee RECOGNITION OF oF Prosiem In AND 

geotechnical engineer needs to be to identity the 


ia ~ could collapse and to determine the amount of collapse that may occur. Soil — 
deposits most to are: Loose fills; ; 2) altered wind-blown 


| 


sands; (3) hillwash of renalaaeal and (4) decomposed | granite or other 

Gibbs (10) has proposed the use of the natural dry density and liquid | limit 
as criteria for predicting collapse. His method is based on the premise that 

_ a soil, which has sufficient void spaces to hold its liquid limit moisture at — 
saturation, is susceptible to collapse upon wetting (see Fig. 2). Soil densities 
g plot above the line shown in Fig. 2 are in a loose condition, and ‘when 5B 


- fully saturated will have a moisture content greater than the liquid limit. ue a 
\ very simple test which can be performed in the field is the ‘‘sausage”’ 
_ test A hand-size sample of the ‘soil to be tested is broken into two pieces, 
_ and each is trimmed until they are approximately equal in volume. One of 
them is then wetted and molded in the hands to form a damp ball. The two 
- volumes are then compared again. If the wetted tell ‘is epi smaller, then 
ne 

at ate. 

wil 


= 
OG. 3. —Typical Collapse Potential Test Result 
A more satisfactory test can be done with a consolidation test, as suggested 


by Knight (15). A sample of the soil is cut to fit in the consolidometer | ring, 

_ and loads are applied progressively until about 200 kPa is reached. At the — 

of this loading, the specimen is flooded with water and left for a day, and 

_ the consolidation test is then carried on to its maximum loading limit. The 

: "resulting curve is shown in Fig. 3. The atts potential is then defined as: 


us 


— 


in which Ae, = change in void ratio upon wetting; and e, = natural void ratio. 
a he collapse potential may also be defined as: 


| 
q | 
ig 


AHe yin height. upon wetting; and Ho = 
_ collapse potential is only a guide to the collapse which may be ¢ encountered. 
Jennings and Knight (13) have suggested some values for collapse potential — 
which are shown in Table 
(Carcuation of Couarse 


_ A method to predict the collapse settlement of a soil for pe purposes 


_ has been proposed by Jennings and Knight (13), using the results of a double — 
_ consolidation test. The use of this test will will even not tonly a qualitative ¢ determination a 


edjustment of me. curve 


soi/ af nat. me 


sooked ofter 24 bedding or /kPa 


ent » Consolidation Test and / Adjustments for Normally Consoli 
_ of the possibilities of collapse, but also quantitative information to allow pe : 
to be made of the magnitude of the collapse. 
_ Two undisturbed similar samples, preferably block samples cut by hand (those — 
taken | using a thin- walled tube are usually disturbed), are tested. Both samples — 


are carefully rye into consolidometer rings, and placed in a consolidometer _ 


i! 

| 

: _is inundated by flooding with water, while the other sample is kept at its natural 4 
left for a further 24 h. If the soaked 
a rebound effect is being experienced — 


= the soil may be of the h heaving type. Generally, most collapsing soils of ne 


4 sandy type will show little or no > change ae this 24-h scaking period. 


. iene and that the two curves do not start from the same point ports 
_ Fig. 4). The total overburden pressure, P.,, at the depth of the sample, is calculated 
| plotted on the e versus log p curves obtained from the two. tests. The g 


precompression pressure, P., is found from the soaked curve and compared he. 
Normally ‘Consolidated Case.— —For the case in which P. /P, is ome to 0.8-1.5 | 
(essentially normally consolidated), compression is considered to take place on i 
the virgin ¢ curve ‘e and the natural moisture content curve is adjusted tc to ) the ¢ e, 


shown i in Fig. rw 


= 
collepse 


FIG. 6. —Double Consolidation Test an ont Adjieeene for Overconsolidated Soil (13) ; 


_ If the loading is increased to Ap, then the unit settlement for the soil, without 
change of natural moisture content will be Re 


— 


If the applied loading remains constant and the onl increases in water content, 
(ene the unit additional settlement willbe 


If the values are used with various selected layers in the profile, then the © 
total settlements under both conditions may be calculated. It should | be noted 


ms if a is onal, then the accuracy of the estimates will be less. Mase 


g 
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Overconsolidated Case. case of an (P. /P, 
1.5), the adjustment to the curve largely follows the ordinary settlement nt 
- tion practices. The only difference between th the two cases lies in the determination 
Jennings and Knight (14) conducted several field studies and found that the _ 
e agreement between observed and predicted collapse settlements was good. There | 
. are several but infrequent cases in which tests have been conducted with samples < 
_ having an equal or nearly equal (e,),. This can be attributed to disturbance 
r - during sampling. Since there is a variety of mechanisms supporting collapsible _ 
soils, and a number of predictive techniques, a field test should be performed — 
whenever possible to check the laboratory analysis and the efficiency of the 
wetting technique used. The major drawback of field tests however, is that 
they are expensive and time consuming. Recent studies by Riggs (19), and Riggs — 


_ and Schmidt (20), have indicated that can be 


rg 


OF 


from field plate load tests, laboratory testing, or « stendard penetration tests. a 
The initial moisture content is | is a critical parameter in determining the most 


in collapsible soils | ‘is results from 
laboratory or field tests can be used to predict the amount of settlement that 
‘can be expected. In many cases, deep foundations (piles, caissons, etc.) may _ 
be required to transmit foundation loads to suitable bearing strata below the 
4 ‘collapsible’’ soil deposit. In cases where it is feasible to support the structure © q 
: on shallow foundations on or above collapsing soil, the use of continuous strip F 
_ footings may provide a more economical and safer foundation than isolated — 
footings. Differential settlement between columns can be minimized, and a more 
equal distribution of stresses may be achieved with the use of strip footings. — 


| — 
FIG. 6.—Continuous Footing Design for Collapsible Soil (26) 
_ A practical procedure for the design of a strip foundation, suggested by Zeevaert 7 
_ (26), is shown in Fig. 6. The foundation system is comprised of reaction beams - * 
__ formed by footing beams and load balancing beams in the longitudinal direction; = 
. 7 the load balancing beams are reinforced to make the system sufficiently stiff. | 


schemes | | may be improved, i ie. ., the vertical displacements may be reduced, 
and more load carrying capacity gained with the use of compensated footings — 
(26). If the footing area becomes greater than 50% of the entire area of the — 
building, then a mat foundation should be considered for the entire foundation. 
Treatment MerHoos 
; — Ons some sites, it may be feasible to apply a pretreatment technique in order 
to either stabilize or cause collapse of the soil deposit prior to construction 
of a specific structure or facility. The amount of treatment and type depends 
TABLE 2. —Methods of Collapsible Foundation Soils 


treatment desired Foundation treatment method 


— 


Current and Past Methods; 
-Moistening and compaction (conventional extra- -heavy, 
impact, or vibratory rollers). 
Overexcavation and recompaction (earth pads with or 
without stabilization by additives such as lime or ce- 
ment); Vibroflotation (free-draining soils); Rock 
umns (vibroreplacement); Displacement piles; Injec- 
cord tion of silt or lime; or Ponding or flooding (if no | 7 


Any of the aforementioned, or combinations of the | 


‘a aforementioned, where applicable; Ponding and in- 


filtration wells; or Ponding and infiltration wells with 
_ Possible Future Methods: gel 
Heat treatment to solidify. the s soils i in place; Ultrason- 
_ ics to produce vibrations that will destroy the bond- 
ing mechanisms of the metstable soil; Chemical addi- 
tives to strengthen the bonding mechanism of the 4 


_ methods of application); or Use of grout-like addi- 


tives to fill the pore spaces before solidification. 


on the depth of collapsible soil and the support requirements for the proposed 
facility. A great variety of treatment methods have been used in the past. Table 
is a summary by Bara of various "treatment methods used. 
technological advances are necessary before heat treatment, ultrasonics, or 
ey additives (other than lime or cement), are feasible, 2 
4 _ The Russians have conducted extensive studies with chemical stabilization 7 
ie techniques (21,22). The methods currently employed are: (1) Gaseous silicatization 
of sandy and loessial soils; (2) strengthening of carbonate cements by polymers; _ 
and (3) chemical strengthening of alluvial soils by clay-silicate solutions (21). 
_ The gaseous silicatization treatment involves a mixture of soil, carbon dioxide, oe 


tL 
ae) 
| Hew 
| 
om of 
and a 1um siicate solution. ests on noncardonate- ype Sandy Sous 


COLLAPSIBLE SOILS 
_ pretreated with carbon dioxide have shown strength increases of 20%- 25%. 
Recent investigations by ‘Sokolovich (22) have shown that Stabilization of 
loessial soils by treatment with ammonia is possible. In this 1 method, gaseous 
ammonia is injected via boreholes into loessial soil prone to slump-type settle- 
ments. The ammonia is absorbed by water films of the loessial soil and reacts — 
with its absorbing complex. As the result of an exchange reaction with the 
absorbed calcium, highly-dispersed calcium hydroxide i is formed. During eso 
_ tion saturation of 1 m’° of loess, 5 kg-8 kg of gaseous ammonia is adsorbed, — 
‘leading to the p precipitation n of 10 kg-17 kg of highly-dispersed calcium hydroxide. — 
Reaction of the precipitated calcium hydroxide with the silica and colloidal 


silicic acid of the soil leads to formation of a calcareous-siliceous binder which © 
“stabilizes bos Yo aninastion is . adel 


In the case of gaseous stabilization of loessial soils, the lime precipitated 
by the > exchange reaction is located on the surface of the colloidal absorbing | 
- complex. In this case, only the skeleton of the loessial soil is stabilized, and ia 
its natural structure is strengthened, whereas the porosity of the stabilized loessial 7 
soil is virtually unchanged. In contrast to ordinary strengthening methods, the — 2 
aim of which is to impart high strength to the ground, gaseous stabilization 
& of loessial soils leads “— to a decrease = tendency f for collapse- type 


4 semi-arid regions. Since many of these ar areas are b beginning to develop economically 7 
and industrially, it is ; becoming increasingly more important to understand and — 
predict the behavior of these soils. Structures or facilities (canals, embankments, 
etc.) of any nature constructed on these soils should receive careful attention. — 


fe Inm many instances, collapsible soils are are ‘undetected during initial soils investiga: — ! 
‘g tions. An increased awareness of the existence of these soils and more full-scale 
research are needed to evaluate th the > present Predictive te techniques and treatment 
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attention of 1 many engineering firms and contractors competing for a share of 
the Saudi construction program. A prior knowledge of ground conditions is ; 
necessary not only for a proper basic design but also to assure that the construction 4 
_ iscompatible with environmental geology. This paper describes ground conditions m4 
7 and their implication i in design and construction in the Kingdom of Saudi Arabia. 


~ Many of the concepts presented are applicable in many other areas of the 


To assess ground conditions i in Saudi Arabia, a general knowledge of topogra- 
J] phy, climate, and bedrock geology is essential. These factors control the 
: development of surface soils and rocks and, therefore, control ground conditions. 
Saudi Arabia, approximately 830,000 sq miles in area, occupies most of = 7 
Arabian Peninsula (see Fig. 1). The terrain is quite irregular. The principal — 
_ physical divisions are the relatively flat western coastal area, the western mountain 
zone (Asir and Hejaz), the central plateau which declines northward in elevation 
to the Plains of Jouf, and the relatively flat terrain of the Eastern Province 
bordering | on the Persian (Arabian) Gulf. Several depressions, believed to be | 
- the result of tectonic movement, are superimposed. Examples are re the Rub a 
- Khali (Empty Quarter) Desert, the Dibdibba Basin, and Wadi as Sirhan Basin. 
_ Other basins such as the Tabuk Basin that extends to the Jordanian border 
are believed to be the result of erosion, = = 


tus.) 


_ 'Vice Pres., Converse Ward Davis Dixon, 91 Roseland Ave., P.O. Box 91, Caldwell, ; 
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x ‘Note.—Discussion open until August 1, 1981. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on June 3, 1980. > 
This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings — eo 
of the American Society of Civil Engineers, © ASCE, Vel. I 107, No. GT3, March, 1981. 
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a The development of surface sediments is influenced by the arid and semi- atid 
oh hot climate in the Kingdom. In addition, this climate has significant effect on 
materials = in superstructures. Chemical reactions affecting durability of 
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NUM 
 WADIS: 1. WADI AL BATIN, 2. WADI BIRK, 3. WADI BISHAH, 4. 
-DAWASTR, WADI HADHRAMAUT. 6. WADI NISAH, 7. WADI RANYAH bed 
AR RIMAH, 9. WADI AS AH'BA, 10. “WADI AS SIRHAN, nd 
WADI TATHLITH. ‘ 
LANDFORMS: 12, ADEN HINTERLAND, 13. AL 'ARAMAH ESCARPMENT, 14. AD 
15. AD DIBDIBAH, 16. GREAT NEFUD, 17. HADHRAMAWT 
ss PLATEAU, 18. HADHRAMAWT VALLEY, 19. AL HIJAZ PLATEAU, 20. a 
7 Ke PLATEAU, 21. JABAL SHAMMAR, 22. AL JAFURAH, 23. MUSANDAM 
ss PENINSULA, 24. NAJD PEDIPLAIN, 25. OMAN MOUNTAINS, 26. AR Rus! J 
Gone AL KHALI, 27. SABKHAT MATTI, 28. AS SUMMAN PLATEAU, 29. SYRIAN 
; PLATEAU. 30. TINAMAN, 31. TUWAYQ ESCARPMENT, 32. UMM AS SAMIM, a 
EMEN HIGHLANDS. 


FIG. 1 —Mapc of Arabian Peninsula Showing Principal Geomorphic oii [Courteey. 
concrete are accelerated by higher te ‘temperatures, and extreme temperature 
fluctuation causes expansion and contraction inc concrete elements and pssceceoy 7 
_ paving. Some of the effect of the hot climate on n concrete construction (compare aed 


ie a ages climate) is shown in Fig. | 


| 
| 
| 


The. Arabian Peninsula i is bisected by the Tropic of Cancer. —— summer r months, 

_ daytime temperatures frequently exceed 40° C (extremes above 50° C have been © 
seconded). Daily variations in temperatures of 15° C are not uncommon. Winter 
_ temperatures are about 20°C on the average and in some elevated regions, _ 
like Qasim, temperatures may drop to below zero. Precipitation is quite low 4 
(some areas receive little or no rainfall for several years) and evaporation rates — 


“floods. . With little or no soil or ; vegetation | cover in hilly areas, high runoffs — 

occur, resulting i in in flooding and the development of of ‘watercourses (wadis). The 

flash floods can cause significant damage to property and | construction sites. _ 

_ Recognition of the flood potential is particularly significant in assessing the = 

need for protection of highway embankments. Some of these wadis change s 
_ course over the years, and thus a proper assessment of the need for highway 


_|COMPARED TO TEMPERATE CLIMATES 
NTHS __YEARS A 
EXTERNAL 


EXCESSIVE CHLORIDE 
CONCENTRATIONS 


_ protection in the 1 term is is required. Many small culverts are usually preferred 


_ handbooks may lead to gross underestimates if the ‘‘flash flood’’ — 
- is not recognized. The major wadi flood areas are shown i in Fig. 4. st.  ovaddying 
‘The development of salt flats (sabkhas) or salt bearing soils in Saudi Arabia — 
7 is mostly attributable to climate. These saline soils are evaporites that form 

& desert basins ¢ or wadi lakes, 0 or form i in in other areas by fluctuating groundwater. — 


a to intense evaporation of salty water drawn to the surface by capillary q 
action, or due to evaporation of tidal water that periodically inundates the —— 
_Evaporite salts are also present in local bedrock and therefore shallow groundwater 

or water drawn to tl to the surface by capillary action is | usually — The ieee 7 


... high, exceeding precipitation by a factor as high as 30. While the | 

a 

| 


and the less soluble gypsum. Saline soils could cause corrosion of steel and } 

_ deterioration of foundation concrete. Structural collapse of salt cemented soil — 

= occur upon prolonged exposure of the soilto water, ee 

i Climatic effect is also evidenced by the loss of support of pavements in 

4 coastal areas. Evaporation is prevented or retarded by impervious | asphalt © 

-_ surfaces, and capillary action increases moisture content beneath the pavement. 
This destroys capillary tension and saturates the subgrade. For subgrades sensitive 
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to saturation, , pavement instability occur. For this reason, most 
use ‘‘soaked’’ California Bearing Ratio (CBR) for pavement design. 
4 one- third of | ‘Saudi Arabia is underlain by a Precambrian 
_ Basement Complex of igneous and metamorphic rocks (Arabian Shield) (see a 
Fig. 5). The remaining area (nearly 600,000 sq miles) is underlain by oe 
rocks of various ages (Cambrian to Quaternary) deposited over the basement 
= rocks with thicknesses of several I hundreds of meters and gentle dips (to north, 
south, and east) of 1°-2°. Many of the sedimentary formations are 
_ and contain water and oil. The source of groundwater is primarily the runoff Y 
from elevated areas of the Arabian Shield where rocks are impervious and 
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ARABIA 
"where: is confined to wadi alluvium. 


= _ The western coastal plain along the Red Sea is low in elevation, in many 
areas covered by the salt bearing deposits, and is flooded as a result of rainstorms. — a 
In some locations, salt bearing soils are absent and in their place are terraces . ; 

of coralline limestone overlooking the Red Sea and, further inland, localized 
deposits of argillaceous silts and gypsum outcrops can be found. The coastal ® 
plain is usually protected or isolated by offshore bars and limestone reef chains 
which often form lagoons. Elevated limestone reefs often outcrop, and raised & 
_ beaches, s sand dunes, and d occasionally gypsiferous sands s and gravels ar are encoun- 
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—Main Syeteme of Saudi Arabi Arabia (3) 


a West of the coastal plains, offshore coral banks and reefal limestones can 
be found either outcropping or covered with mostly carbonate sands. _ 74, em 
_ The coastal plains of the Red Sea rise gently into an eastern plain (Piedmont), jf 
consisting of granular soils (gravels, boulders, cobbles, sands, and silts), overlying : 
Tertiary (limestones, shales, siltstones, and sandstones) or crystalline rocks. a 
a In some cases, especially near the fringes of the coastal plain, the granular > 
- soils were deposited during floods originating from the mountains to the a 
; ‘ and thus overlie soft coastal plain deposits. Salt beds, e.g., halite in a 
can be found in the Tertiary rocks. This Piedmont zone extends to the mountainous _ 
_ Arabian Shield. Gravels and boulders are usually found in the high zone of 
cy the slope, and finer soils prevail with increased distance down the gentle slope. 
‘The Precambrian rocks of of the Aral Arabian Shield | are } mainly of of boa 


| 


“metamorphic type, e.g., in Khamis Mushayt are found granites and gneisses: 
and in Taif are found granites, schists, amphibolites, and gneisses. The — 
form a dome- -shaped topography and often are overlain by younger volcanics 
7 (mainly basalts) or thin (less than 10 m) strata of alluvial sands and gravels, — 
> both. In addition, late Precambrian nonmetamorphosed red arkoses and 

; conglomerates can be recognized in isolated basins or folds over much of the — 
4 Uplift of the shield is believed to have occurred i in the late Precambrian 
- Age. The eroded Precambrian basement is ssions of mostly 
Z sandstones, overlain by shales, with the lower sandstones ‘extending northward 
_ x into the marine Cambrian rocks of Jordan. Erosion has continued since the 
Precambrian time producing present-day alluvial fans and wadi deposits. 
“a Sedimentary rocks (mostly limestone, but also sandstones and shales) of various 
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FIG. 5.—Generalized Geologic Map of of Arabian Peninsula (Modified) (Courtesy, 
: = cover (a few meters, resulting from either rock weathering or windblown 
P sand) or with a thicker (several meters) cover of | wadi alluvium or residual | 
soils (mostly stiff clays) that are also the product of 1 rock weathering. 
_ Sedimentary rock deposits dip gently to the east (1°-2°) and extend to a 
Gulf. These rocks have been subjected to uplift and erosion (probably in the 
Tan Age) and the more resistent or massive rock units form west-facing — 
_ escarpments such as Jabal Tuwayq near Riyadh. Sands and gravels of the 
+ Quaternary Age can be found in wadis between escarpments such as wadi Hanifa 
west and south of Riyadh and along escarpments near Riyadh. Deposits of z 
Tertiary sands and gravels and lacustrian limestone can be found along erosion * 
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SAUDI 
i. consisting of salt bearing soils, inter sabkhas, areas of sand dunes, 


@ outcrops of limestone and calcareous sandstones. The near-shore Gulf — 
"consists of shelly and clayey sands, corals, wr and sandy limestones. 


- sands and rocks in the Gulf Plain and floor is related to Quetaiey marine — 
_ Foundation Zones.—From a geotechnical viewpoint the Kingdom < can be divided 
into four major zones as follows: gravel or 
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«FIG. 6. —Subsurface Conditions, Zone A-1 


‘1, Zone A—Western coastal plain consisting of Zone A-1, generally | covered © 
by soft sabkha coastal deposits, and Zone A-2, covered by alluvial deposits 

-.. 2. Zone B—Arabian Shield extending north and east to the coastal plains 
= . Zone C—Sedimentary rocks from Shield to coastal plains of the Gulf. 


_ Foundations. —Determination of the most suitable method for foundation 
support depends on on the subsurface conditions | at a particular site, considered » 


AAR OB 
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a Zone A Western Coastal Plain). —Soil conditions in vary 
generally soft and loose surface sabkha deposits in the western portion (Zone _ ; 
; A-1) to the generally dense alluvial fans and wadi deposits extending = 
to the foothills of the Arabian Shield (Zone A-2). The sabkha deposits, varying : 
in classification from sand to clay within the same site, can be found with 
variable | ‘thickness but usually not more than 10 m and typically less than 3 
m thickness, and may be underlain by dense granular soils, reefal limestone, 
_ or loose sand and soft clay layers with coralline rock layers or ledges. The 
; soft and loose coastal soils may range in thickness from a few meters to 20 


m or greater, so that a typical could be misleading. some Sites, e.g., 
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- Figure 6 depicts conditions encountered by various borin s for a large e site 
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‘in the vicinity of Jeddah. The density of the coral formations vary erratically — 

as well as the occurrence and stratigraphy of soil strata. In a marine euviccament J 
_ such as Zone A-l, extrapolation of data to an adjacent unexplored location 
Methods used for foundation support in Zone A-l vary, depending on on the 
loads and performance requirements, e.g., preloading and conventional founda- © 
tions have been used successfully in Jeddah to support | story-2 story structures. _ a 
The time of preloading is generally less than 6 months but depends on the — 
stratigraphy — at the specific locations. As in other places, deep stabilization — 
eae may be used to densify the loose granular soils but a high shell 
content may retard the usefulness of this technique since it may not be possible _ 


to achieve the required density, 
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The use of piles or ‘other deep foundations should consider the extreme 
vations in the rock surface, e.g., in Jeddah, and the presence of hard rock — 
ledges or layers above loose soils. A properly conducted site investigation prior _ 
* design and construction should reduce the risk of surprise related to the 
_ extreme variability of subsurface conditions, and identify potential problems 
as posed by high water table. Requirements, design, and construction problems “ 
in Zone A-1 vary along the coast, e.g., the width and thickness of this — 
i Yanbu is less than in Jeddah. Piles driven onshore in Yanbu can be expected 
_ to be shorter and may bear on the dense granular stratum beneath the upper 
loose and soft stratum (and not necessarily on rock of unknown thickness and — 
extent with large variations in pile lengths, as is sometimes the case in Jeddah). - 
Zone A-2 ‘Soils dense sands and or hard clays) present no 


OR ALLUVIUM) 


‘a DENSE SAND TO GRAVELLY SAND OR 5 
a" V7 SILTY AND CLAYEY SAND (RESIDUAL SON 


SEDIMENTARY ROCK (BRECCIA. SAND- 
STONE OR LIMEST: IGNEOUS GRAVEL AND COBBLES 
- 


(GENERALLY GRANITE) 7 


¢ 
unusual design an and construction problems. Example stratifications are » depleted 
in Fig. 7. Spread or strip footings on natural soil or compacted fill may © be 
used. In some coastal locations, e.g., Jizan, soft surface soils are sometimes — ‘ 

- encountered and associated problems are compounded by the presence of soluble 
salts in the ¢ sedimentary formations. For these conditions, support by piles bearing 
in layers not affected by solution have been suggested. However, proper site 
“a and facility drainage to prevent prolonged exposure to water would be more 
; appropriate. The same applies to soluble gypsum bearing soils in other coastal J 
4 areas, e.g., Yanbu. In addition to the potential for solution, the gypsum bearing ~ 
Soils havea structure that is prone to collapse (resulting i in significant settlements) — 
‘upon on prolonged exposure to water. "Identification of such soils” is critical to 
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‘procedure to evaluate the potential for collapse under prolonged 


exposure to moisture. is _— use of the | plate load test on a ‘soil under « dry 
and saturated conditions inthe field. 
- Zone B (Arabian Shield).—Problems related to the presence of soft soil : . 
conditions are not common in this zone. A geotechnical section (see Fig. 8) 
- would show relatively shallow cover composed of in-place weathered soil, wedi 
alluvium, or windblown sands and silts over rock (mostly igneous and metamor- 
s phic) with varied degrees of weathering. In some cases rock weathering has 
7 ‘resulted in layers of residual soil (generally dense sand and gravel but occasions 4 


hard clay near the eastern border of the Arabian Shield). At some sites, e.g., __ 
: in Taif, soil cover (dense wadi sands and gravels) as deep as 7 m can be ws PS 
a _ encountered. The use of the saturated penetration test in this zone for site __ 
inenignten should be supplemented by test pits, trenching, and shallow seismic ae 


= nr surveys. Because of expediency, hollow stem augering is — n used 
_ for the borings. The classification of the materials based on cuttings is poor, 
_ especially when distinguishing residual soil from | completely weathered rock. 
For a typical site, individual footings could bear on dense natural soil (or 
compacted fill), or residual soil or rock. Depending on the magnitude of design 
loads, some special site preparations may be required to minimize differential 
settlements. Evaluating rock rippability and blasting requirements are of primary 
significance in any in this z zone. In general, rock encountered 


can be “expected. ee tg may be used effectively to evaluate rippability in = 
“conjunction with other techniques (seismic velocities, point load test, RQD, : 
etc.). As the rock is relatively i impervious, shallow groundwater conditions may 
be created, depending on site grades, as a result of irrigation, e. >in ‘Taif, 


{ 
Sea 
FIG. 9.—Cavity Development, Limestone of Arab Formation, Riyadh 
| 
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the Arabian Shield to the eee plains of the Gulf. Foundation —— 
in the zone are influenced by: (1) The presence of cavities in limestone bedrock, 
e.g., in Riyadh, Dammam, Buraydah, and other locations; (2) the presence of % 
gypsum bearing soils in basins and other poorly drained locations, e.g., Tabuk, 
Qaysumah, Al-Kharj, Riyadh, _ and other locations; (3) the presence of inland 
_ sabkhas, e.g., in wadi Sirhan, southeast of Buraydah, and other locations; and 4 
- (4) variability of foundation material within the same site, e.g., rock versus 
- Thick sequences of sedimentary rocks have been deposited i in this zone since 
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Y EXCAVATED SURFA\ 
ew SURFACE 


be expected in between the more resistant limestone formations. Some —vtntocnl 
s g., the Arab Formation, contain anhydrite beds i in which past surface dissolution _— 


the (Fig. 9) is most frequent along ‘the ‘north-south axis 
a the city of Riyadh. Figure 10 shows details of a solution channel found 
ad. s&s a site in Riyadh. These cavities are often encountered (during construction) 

_ between locations where borings showed no cavity development (the subsurface — 

section in Fig. 11 for a site in Riyadh typifies the variability of the limestone 


aon ae upper Tock su surface | in 1 Zone C is usually shallow (typically 0 Lod 


underdrains may be necessary, depending on the design particulars, 

— 
— 
— 


Numerous cavities have been encountered within 10 m the surface. 
Unfortunately, detailed bedrock verification techniques are expensive and may : 

_ not be effective, depending on particular site conditions. Usually, reliance is 
on borings and cheaper probes using wagon drill. Other 


and verification ( (bys wagon drill probes at foundation locations) during construction 
are quite important and sometimes imperative if foundation failure is to be 
avoided. Another measure useful for design of lightly loaded structures is to ag 
use properly sized and reinforced strip footings (instead of spread footings) 
_to span undetected cavities. When spread footings are employed, it is preferable — 7 
and relatively inexpensive to verify conditions below each footing during oS 
- construction. In case nae are detected, excavating of the cavity zone and 


ELEVATION IN METERS 


Xb 


(ROCK QUALITY DESIGNATION 


- conditions. As is the case in Zone B, evaluating rippability and need for blasting 
is essential and can be accomplished by site and laboratory investigation 
techniques, ©.8., Tock coring and logging, shallow seismic refraction 


a developments have also been encountered in — Dammam, and 
_ Hofuf and some cavities have even occurred near Al Batin beneath the overlying 


|. = 
— qf 
«FAG. 11.—Typical Subsurface Section, Riyadh __ 
A 
_ profiles (alluvium or alluvium over residual soil) may be encountered. In a Ad : 
7 _ site situated in downtown Riyadh, the soil profile consisted of 10 m of dense 
1 gravelly sand underlain by 10 m of stiff plastic clay over marly cavernous | 


¥ lanteteae. Figure 12 exemplifies some soil profiles i in various a Cc sites. 
4 _ Some sedimentary formations, e.g., the Russ Formation in the eastern region, _ 
1 contain thick layers of marls. Attempts are usually made to auger through the cd 
; : highly indurated layers (or yrs and if this fails, coring is then attempted. - 


q ian a Allowable bearing \ values < on ‘marl a1 are e variable. Allowable bearing values = 
= for footings on the Rus Formation have varied from 3 tons/sq ft to as ‘ 
high as 9 tons/sq ft. In a typical profile, layers of rock or rock fragments 
are embedded in a dense soil matrix. The marl is often used as compacted ; 
& to support spread footings. When this i is 3 done, op site and facility drainage S 
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Subsurface Conditions, Various Locations in Dene c 


 Gypeem bearing soil can be encountered in poorly drained areas om 
"jm Zone C. Support of footings on gypsum bearing soils can be made safe 
by extending footing depth to zones where the percentage « of gypsum can bee 


grading the site for water to drain away from the footing, concrete paving, 

— enclosing water carrying pipes in impermeable galleries, etc. The most economic 

_ solution would depend on specific site conditions. It should be mentioned that 

gypsum bearing soils may occur in localized areas and for a large site, they 

Bear be missed by borings. Therefore, proper verification during construction 
is essential. The percentage of gypsum in soils used for compacted fill should 

be controlled depending on the grading plans and the moisture environment. 
_ As explained earlier, plate load tests may be used to evaluate the potential 


| 
| 
— 
| 
for collapse of the 
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a of coastal sabkhas. The depth of inland sabkha zones varies s and the * 

current highway practice in the Kingdom is either to avoid these areas or a 

Zs the weak material by fill. Typically, bridges are pile a_i te 
filling scheme is not gatirely satisfactory as significant roadway Settlements a a 


_ Superimposed on Zone C are extensive areas of sand dunes and areas of ie 
‘ alluvial deposits. The alluvial wadi deposits (fine gravels and sands), ¢.g., in 
AbKhasj, may be several meters thick overlying bedrock and usually 


Rs = REFUSAL TO STANDARD PENETRATION: Tuan 
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DENSE SAND AND SILTY S. SHELLS: 


SANDY SILT) TIALLY CEMENTED 


WEATHERED AND AND OR HARD 


Sand dune areas may be loose, eapecially near the surface. Mobile par 
_ present maintenance problems for roads. In addition to design techniques (aligning 
_ route considering prevailing wind direction and location of sand source, ~ 

or no crown, etc. ), techniques such as cutting of dunes, fencing, planting, ,paving, 
_ or oil spraying and impregnation may be considered for dune stabilization. None = 
: $ of these techniques are problem free and as most are relatively expensive, Fry 
continuous maintenance may be cheaper and more effective. = = 
‘a Zone D (Eastern Coastal Zone).—Most of the eastern coastal plain is categorized F 
. a by the presence of salt bearing soils or sand dunes. Typically, subsurface ee 
conditions are defined by an upper zone (usually less than 3 m thick) of soft a 
and loose deposits in low-lying areas near the shore and generally medium 


dense to dense sands below. of stiff to or 


— 
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Artesian groundwater conditions also may be encountered with piezometric level : 
_ several meters above ground surface. Some soil profiles below the mudline _ 
_ in the offshore portion of Jubail are exemplified in Fig. 13. Figure 14 depicts ~ 
conditions along or near the shore. Surface limestone caprock is often encountered 
along the shore with thicknesses up to 2 m, which could present a 
during dredging. Soil conditions below the mudline can vary from loose to 
very dense depending on location. ove 23 
ene of suitable methods for foundation support in Zone D again ms sap a 


snot 
gration 
ON: Penetration Resistance blows/30cm. 
LOOSE SAND OR SOFT CLAY (SABKHA) 
MEDIUM DENSE TO DENSE. SANDS AND pion view 
sary Clays TO CLAYEY SATS. 
or 


in Zone D is generally less than in Zone A. Thus, excavating the soft or loose + 
soils and replacing by fill probably would be more cost effective in Zone D. 
Shallow foundations bearing on natural soil or compacted fill are feasible in 

the medium dense to dense sand areas. The same methods of support may 
be feasible for supporting light to medium heavy structures | on properly filled — “4 

and preloaded loose sabkha areas. Pile foundations may have to be used for 
- heavy structures depending on the sensitivity of structures to settlements. Soil | 


improvement techniques, e.g. , Vibroflotation, terraprobe, and dynamic compac- 


j 
L 
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tion, may be considered if shelly layers are prevalent. The 
of these techniques should be assessed by construction of carefully controlled — 
___ As indicated previously, loose soils may contain large volumes of shells and — 
these could preclude the effectiveness of the stabilization method selected. The i 
_ underlying denser soils may also contain large concentrations of shell that t may “@ 
affect both pile installation and support. A properly conducted site investigation A 
is therefore prudent. The presence of rock ledges and lenses that are not continuous 
should be considered in planning underground structures, dredging, and design a 
Construction Materials. —The availability and quality of construction materials 
are dictated primarily by bedrock geology and the presence of alluvial deposits. *, 
- The alluvial fans (Zone A-2) usually provide a good source of natural aggregates, 7 
while the Arabian Shield (Zone B) provides a large source of rock which is 
suitable for crushing into aggregates. The windblown sands and beach sands 
in Zone A-1 and Zone D often are unsatisfactory for use as fine aggregate a 
because of either poor gradation, excessive fines, or excessive salts. Quarry 
‘sand is often used because of improved ‘gradation. Data on physical and 
engineering properties of soils and rock at various locations are shown in Table 1. 
In general, it can be stated that within a distance of 30 km of a site, it © 
_ should not _ difficult to find a suitable source for aggregates in Zone A or 
ee B, g., Suitable sc sources in the vicinity of guerel can 1 be found at Jabal 7 


extends from the Taif mountains), and Ras at Belejh for a sand borrow area. 

_ In the Yanbu area and to a lesser extent in Jeddah, construction materials are 
abundant and there appears to ‘be no problem of finding and transporting — 
construction materials. However, resources of fresh water are nil or limited, - 
* and salts that render the materials unsuitable in natural aggregates without 
processing may be encountered. Aggregate washing may not be feasible and 
processing may have to be conducted in the dry. 

_ In Zone C and Zone D, the availability of construction material is also dictated © 

| _ by the quality of rock and wadi deposits. In several quarries, the Jubaila limestone _ 
_ formation in Riyadh provides a satisfactory source of crushed stone aggregates. 
At many locations the upper meter or so is case hardened and is of relatively 
better quality to the extent that it is used as a facade material for buildings. 

Wadi sands and gravels, e.g., in wadi Hanifa, Jubaylah, Sulbukh, Al-Ammariysh 
and others, may be satisfactory depending on, among other factors, the degree 
of their alkali reactivity and fines content. Most of the sandstones and shales 
in Tabuk do not usually meet specifications for concrete aggregates and reliance - 

is made on wadi deposits in the area, e.g., wadi Al Uthayli and wadi Al Altheeli. 
In some locations, aggregates have to be obtained from distant off-site sources 
as exemplified by the town of Taima, where the aggregates used are obtained 

from wadi Al Qa about 100km south ofthe town, 

ew Unfavorable bedrock geology is the key to the lack of "good quality A. 2 ely 


| 
_ of aggregates from limestone rocks of the Al Khobar member of the Dammam a _ 
_ Formation. Proper quality control should be exercised to avoid contaminated _ 
: - or weak aggregates. Another potential aggregate source is the limestone of the | aa 
km | west of the Jubail /Dhahran axis. Some outcrops 
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and Engineering Properties: Various Sites thet 
| |inpounds| in | mum | mum 
Sound- | per cubic den- den- 
ness | foot | cent sity sity 
7.1-23.4 9.6-13 | 23-75 
mol 
in | 2.35-255| 18-5 | 31-51 | 7-75 | 
ALKharj ne | 78 | 26-40 
120-130] 7-9 | 4690 | 
| situ soaked CBR; moisture content 13%-20%. 


q 
brackish water with total dissolved above 1,500 ppm. 
. large construction, on-site water treatment may be required depending on the = 
availability of potable water supply. In some cases, natural waters are of adequate 
quality to be used in construction anes treatment. However, potential sources — 


TABLE 2.—Some Results of and ‘Water 
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|bsboratory 


Sand (SP) 
Tap water 
Sand (SP-SM) M) 
Water 
Sand (SP) 
Water (well) 
Khamis Mushayt | Sand (SM) $s 
f 
Tap water 
Sand 
Water (irrigation 
canal 
Water (tap) 
Sandy clay (CL-SC) 
‘Sand and gravel 
(GM-SM) 
Water (tap) 
Sand (SM) 
Water (farm well) 
Sand (SP) 
Water (tap) 
Sand (SP) 
Sand (SP-SM) 
Al-Batin area Sand(SW) ; i ohm-m 


120 ohm-m 


0.73 chm-m 
1234 ohm-cm 


for water and aggregates should be identified and ak prior to ae design 
phase of a specific project. The timing for testing could be important, especially 
_ whee alkali reactivity of potential aggregates with cement paste is suspected. 
X _ The American Society of Testing and Materials (ASTM) Standards C295 (Pe 
_ trographic Examination of Aggregates for Concrete), C289 (Test for Potential ’ 
_ Reactivity of Aggregates—Chemical Method) and C227 (Mortar Bar Method) 


- are useful in 1 evaluating the potential of reactivity y of aggre aggregates. If the e relatively — - 


“| [em] | per | per 
| 
| 
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test (C289) or petrographic « examination indicate a potential for 


reactivity, then the mortar bar test (C227) should help quantify that expansion 7 


> (and thus the potential for damage) due to potential reactivity. The mortar bar 
test duration should be at least 3 months and preferably 6 months or a year, 
; = on the test results obtained. The writers are not aware of concrete 


failure in the Kingdom directly related to alkali-silica reactivity of aggregates 


with cement paste. Reported probiems in concrete are primarily related to high 

: content of salts (sulphates, chlorides) in concrete constituents or surrounding 

_ environment. It is known that dolomitic limestones in the eastern region contain 
chalcedonic cherts having the capability of producing a reaction. However, in 

7 cases known to the writers, such cherts in commercial quarries, e.g., rimrock _ 


— at Dhahran, are found in such small proportions that their effect is not harmful __ 


to concrete in existing structures. Chalcedonic cherts could also be found in 
_ dune sand. . Recognizing th the potential harmful ll effect o on concrete, it is prudent a 
to investigate concrete constituents from all unknown so sources with regard to ; 
- potential reactivity. The cost for doing such investigations is small. The specifi- _ 
_ cations should limit the salt content in concrete since marginally high = 
im aggregate coupled with marginally high salt content in water would yield — 
Corrosivity and Chemical Attacks.—Examples of sulphate attack of concrete 


are common in Saudi Arabia, especially i in the sabkha and shoreline areas where ¥ 


_ conditions are suitable for capillary action carrying salts to and around concrete 
: foundations. Design and construction safeguards include the use of sulphate — 


af 
resistant (Type V) cement, provision of a capillary barrier beneath the laden ae 
all 


(OF, in some cases, tanking the foundations, i.e., using bitumemastic coating id 
on all sides and the bottom. In most areas of the Kingdom (Zone B and Zone ~ ¥ 
© the use of Type V cement would most probably be sufficient. In several _ 
: areas of Zone B, special precautions are not required as the salt concentrations nil 
2 Corrosion of buried pipes and metallic structures is a real problem. Chemical 
tests. fo or ‘Salt content and laboratory and field electric resistivity tests and perhaps 


at a particular site. The lower the electrical resistivity, the higher the potential 
for corrosion. Soils with resistivities less than 30 ohm-meter are usually considered _ 
very corrosive and embedded metallic structures may require cathodic protection. _ 
Table 2 shows data from several sites. Corrosivity is more severe in a saturated 
or moist environment. Buried pipework, at most site sites, would require protection 
by coating and wrapping techniques, o or, in some cases, ses, cathodic protection mn 


= could be necessary. Soil acidity is not usually a corrosion factor for the sites — 
“listed i in in Table 


Le 


“ ‘Several geologic environments prevail in Saudi Arabia. Problems in design 5 ; 


and construction are intimately related to the environment at particular locations. | 

_ Coastal areas are generally characterized by particularly soft and loose corrosive 

- soils. Limestones i in the interior are characterized by the potential for the presence — 
of mostly dry cavities and solution | channels. Igneous and metamorphic <= 
of good bearing quality chesestatinn the high mountains overlooking the west 


| | 

NOlarityv measurements are neriormed [oO evaluate the notential corrosive natiri 
— | 
— 
— 


t. Superimposed are areas of sand dunes and salt flats which could present ; 
special foundation problems. Availability of construction materials is dictated — 
by bedrock geology varying from relatively good availability in the western» 

_ areas to poor in the eastern region. As is the case in many marine environments, 
offshore conditions in the Red Sea are characterized by varied subsurface — 
conditions at relatively short distances, thus affecting the design of offshore 
- facilities. Offshore conditions are relatively more uniform along the P 

_ The writers are grateful to the staff of the Foundation, Material and Survey 
Branch, ‘United States Army Corps of Engineers, Middle East Division, for — 
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J. Sorensen, who typed the manuscript. The writers are also grateful to ASCE 
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SAFETY OF A CONSTRUCTED 
By T. William Lambe,’ F. ASCE, W. Allen Marr,’ M M. ASCE, 
Civil engineers have always | given priority to the | safety of the facilities they : 
4 design and help construct. Indeed, the public, through its action groups and 
_ governmental regulatory agencies, demands a high degree of safety for certain 


7 


constructed facilities. However, differences exist in the degree of safety employed eo 
ed engineers. Structural engineers use factors of safety ranging from 1-1/2-4_ 
ve higher in designing structures of steel, “concrete, and wood to allow for 


‘unknown and unforseen factors.’’ Geotechnical engineers also employ factors 

of safety for the geotechnical aspects of constructed facilities. Unfortunately, 

economic considerations often require that the geotechnical engineer use lower | 

factors of safety than those used by the structural engineer despite the fact © 

- the properties « of soil, the geotechnical medium, vary more than those = 
of steel, concrete, and wood. The geotechnical engineer usually employs factors a 

of safety ranging from 1.1-1.5 for earth slopes (but higher safety factors for 
As a consequence of these relatively low safety factors, the geotechnical | - 
engineer necessarily employs greater supervision and review of construction 


safety -of their constructed facilities. bang, ‘gohe Radi 
With the recent dam failures in the United States [a coal slag dam at Buffalo — 
_ Creek, W. Va. (Feb. 1972); the Bouldin Dam near Montgomery, Ala. (Feb. 


: __ ‘Edmund K. Turner Prof., Dept. of Civ. Engrg., Massachusetts Inst. of Tech. , Cambridge, 
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1975); th the Teton Dam in ‘Idaho (June 1976); and the Kelly Barnes ‘Lake Dam 
at Toccoa Falls, Ga. (Nov. 1977)], and the resulting loss of life and property, f 
national attention has focused sharply on dam safety. In August 1972, Congress ) 
: passed Public Law 92-367 concerning dam safety; in a December 2, 1977 statement 5 
the President initiated the Federal Dam Safety Program. We can “expect further 
. regulations on dam safety from both state and Federal agencies. Responding = 
to the increased concern for safety of geotechnical facilities such as dams, . 
x offshore structures, and foundations for nuclear power stations, engineers | hav we 
devised and executed systems of “quality control’’ to help reduce — = 
_ The details and the effectiveness of an appropriate safety program for a 
“constructed facility depend on various factors, including; (1) Consequences of : 
failure; (2) stage of life of the facility; (3) desired degree of safety; and (4) 
a nature of the facility and the site. An engineer can tailor the i comprenensive > 
Z safety program presented herein to suit a particular situation, = 
_ This document presents an approach t to “safety of a a constructed facility” a 


“safety program.” Aided by Massachusetts Institute of Tecnology 
‘ "associates, we have, over a period of 25 yrs, developed this safety program 
and applied parts of it to approximately 20 facilities, primarily dams and natural re 
slopes. We continue to lo develop the program as 3 we gain | experience with its - 
—When dealing with important, ‘complex geotechnical facilities, 
; _ experience clearly shows that the geotechnical profession cannot feasibly create 
a facility which completely fulfills all performance criteria for the entire = 
of the facility. Rather, the only logical approach is to o “engineer a facility for 
its design life,’ ” i. e., create a reasonable design, construct the facility according — 
to the design, and then execute an appropriate ongoing safety program. 
7 Factors Dictating Approach. —For two decades, we have investigated the 
_ geotechnical engineer’s capability to predict the performance of a geotechnical 
_ facility. Students of MIT, practicing engineers, and leaders in the geotechnical _ 
profession have made predictions for a variety of situations prior to the occurrence 
of the event. The predictors did not know of the actual outcome at the time 
of their predictions. By comparing predicted performance with measured perfor- . : 
mance, we obtained an approximate measure of the profession’s capability to __ 
predict performance. Our evaluations suggest that the geotechnical ‘engineer ey 


predict with the following reliability; 
Deformation—Predicted vertical deformation +50% of ‘the 
3 deformation. Predicted horizontal deformation = +150% of measured deforma- | 7 


‘Stress—Predicted change of pore +25% of measured change 

of pore pressure. Predicted lateral stress = +50% of measured lateral stress. 4 

_ 3. Stability—Predicted factor of safety = +25% of measured factor of safety. _ 
_, 4. Flow—Predicted flow = + one — of magnitude of measured flow hea — 


Candidate Facilities. — —Experience sh shows the desirability of a long-term safety 


| 
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pes, cut 
=  (2)ex excavations; (3) fluid retention systems; and (4) off shore and waterfront 


facilities. Natural slopes, particularly those involving a reduction of stresses “ 


from cutting, can become less safe over a period of time. The drop in stability ‘ 
esate from a rise in pore water pressure and, for certain types of soil, a 
toss in shear strength. Professor Skempton’s work in England and our own 
experience in Puerto Rico | and Venezuela show natural slopes that failed years — 
after formation. oft svieado Muode edT 
Dams, like natural slopes, can weaken with time. Filling a reservoir tends 
_ to increase the pore water pressures in the dam, thereby lowering its strength. a 
_ In addition, the removal of soil fines by flowing water may worsen with time om : 
* _ Structures designed for earthquakes o1 or sea storms will normally experience | 


| use measurements of performance during the ‘early life of the facility to check 
on its safety prior to arrival of the maximum loading. This procedure, a feature 
- of the program described herein, may prevent a failure similar to that which 
occurred at Sines, Portugal where a large breakwater failed under storm conditions © 
less severe thanthe design storm. 
- Components of Program.—The Geotechnical Safety Program has the following 


1. To determine the cogue of safety of a constructed facility (geotechnical 
2. To help maintain a predetermined degree of safety of a constructed facility. . 
oe 3. To furnish the technical base to permit alteration of the facility to obtain oil 
a different, usually higher, degree of safety. 
Our comprehensive geotechnical safety program has nine 


4 Ensure that the criteria meet appropriate legal and accepted 


Identify the critical mechanisms of 
¢. Identify typical and critical sections, 
od. Review the field and laboratory data used by the designer. Inspec = } 
= the ‘site and run selected tests on the soil(s) th that greatly ; affect the design. . 
e. Examine the prediction ‘methods employed by the designer. Check 
see ‘whether the designer’s methods rest on mechanisms approximating the — 


Check the s predictions of pe performance and them 


used by the designer as well as his prediction methods and parameters, and 


7 a Prepare a Design Semmary which stat states clearly and refer and references Sees Oe loads 


d 
| q 
7 
b 
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h. Use| the methods and parameters to predict 

conditions expected during the early life of the facility, = = 

_ _ Identify major uncertainties and critical aspects of performance. 


wy b. Install instruments which can reveal _ mechanisms and values of key 
‘parameters. The engineer should observe the installation of instruments a 
ensure that they are placed at the correct locations and that they read correctly. 
Periodically check and maintain the field instruments. 
on a. Make periodic i inspections of the facility during its construction. Compare by, 
the actual construction with the design and document the significant features — 
of construction, particularly those features which vary from the original design. 
___b. Assess the significance on predicted performance of such departures — 
om c. Make frequent site visits where difficult site conditions ‘prevail. For ” 
highly complex situations, have a representative present on the site during 
ts a. Make periodic visual inspections of the project, searching for 
evidence of malfunction. “fi of 
Measure important aspects of performance. 
c. Obtain information—data on mechanisms and 
time) field measurements, events, and predict 


and performance criteria. Attempt to explain any differences shown by the 7 

Compare predicted and measured mechanisms and parameters. 

a. Determine and indicate (in terms of ‘safety fac factor or or probability o of ' failure) e 

the degree of safety of the facility. ote 440. ler we 

a —b. State and evaluate any change in degree of safety during the period 
since the last safety assessment. 
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State whether or not the existing of meets performance 
a. Design and assist in executing to the 
of safety to a satisfactory level if the safety assessment reveals an inadequate 


T 
' The preceding nine components constitute the comprehensive safety program. 


The following portion of this paper Wastrates “Design Assessment,”’ ‘‘Perfor- 
mance Evaluation,” and “‘Safety As These from 


I The preceding section listed the components of a design assessment. The 
i thoroughness of and approach | to the desi ment < depends on the nature 


andi importance of any given project. The engineer making the design assessment 
_ Should maintain independence from the organization designing the =— The 


‘Review the design documents and drawings. _ hy 
| 2. Examine the site. 


3. Consider the design with the designer and the owner ni] ah emil 
_ 4, Perform selected check or supporting tests and analyses, or both. 
6. Analyze the assessment v with the congeee and the owner. 


par 


wie silty sand 
Silty sand 
Brown sand 


* 


= Whitesand 
Black siltysand 


Silty sand 


s Value ‘Assesser’s Value 


— 
| 
£ 
i 
store 
7 
J 
| 


= 


“a 


= 


42 


4 


2 


Bulsnp ys0m jIpewes Ao; | WG PINCD moys 


A 


NOILWWHOI30 
30404 


= 
Swoss3ssy 


iS gig 
| 
te 
one 
4 
| 
| : | 
13 BL] 


water, and later clay The. foundation | of layers 
_ sands, and sands. The embankment consists of hydraulically-placed mine mine t tailings 4 
_ with a drain and an upstream face of mixed and rolled silty sands.§ 
ice _ Tables 1 and 2 summarize a design assessment for this facility. Table 1 compares = 
b parameters for permeability « and strength used by the designer with values selected 
@g the assessor. Different | results from lab | ) tests ar and different interpretations — 


of the field conditions « explain the 1¢ disparity i in values. Table 2 presents a detailed = 


evaluation of the design. The first column indicates critical aspects of performance - 
_ for the facility and the second column lists numerical criteria for performance 
as adopted in considerations among the owner, designer, and assessor — 


The fifth column gives the assessor’s evaluation of the designer’s prediction 
methods, parameters, and mechanisms, and the last column lists departures 
of predicted performance f from the owner’s s design cr criteria and indicates required 

This design assessment resulted in investigations whieh 
_ modifications in the design. Performing the assessment helps to clarify the ' 
intended function, predicted performance, and any er problems for 


the owner ond d designer, as well as for the if 


Figs. 2 and 3 show two epens of a performance ‘evaluation mid a recently > 
constructed dam. Figure 2 portrays stress conditions downstream of the dam 
m for a headwater elevation 15 ft (5 m) below the design reservoir level. As 
_ this figure indicates, effective stresses fell below the values specified in the : 


- Figure 3 shows the flow net for a partially full reservoir using pore pressures. 


that flow through the dam met the cr criterion, , predicted “4 flow for full reservoir 


The —T ormance evaluation summarized i in Fi ig 4 =e on the design and | execution 


; 


- 80 


FIG. 2.—Pertormance Evaluation f for | Effective Stresses (1 ft = 0.305 mm, 1, 
ft = 4.88 tons/m?=47.85kN/m*) bist) 30 dal 


| 
{ 
comparison Of periormance Criteria and predicted periormance, along with an 
§ —s Standing of design criteria for the facility). The third and fourth columns summarize ; 
a the designer’s predictions of performance and the basis for these predictions. : 


ELEVATION in 


- of —_ measures which corrected nearly all of the deficiencies. Continuing 


Periodic evaluations check deficiencies not yet fully corrected. 


“Approximate Sarety Assessment For Dam 


| The er engineer often must perform an approximate safety assessment due to 


J 3. A need to define the scope for a a complete AS oP 


_ 4. The project existing at an intermediate stage of completion. 
An approximate safety assessment for dam com to the existing 
conditions of the dam as follows; 
1. A safe condition, i.e., action 
An unsafe condition, i.e., ‘remedial action 
A condition, i.e., further investigation 1 needed. 
approximate assessment normally utilizes existing construction 
records, performance data, anda site inspection. The site inspection team searches 
for the presence or absence of design and construction features associated with 
good performance, or both, plus any evidence of unsatisfactory apr 
An pcan assessment t of dam safety normally does not include additional 5 
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experience permit one to list features associated with good 
and poor performance of dams. We assign a scale to these features to indicate b 
the relative importance of the various undesirable features and malfunctions. 
‘Fig. 5 gives a list which uses both a numerical scale and a descriptive scale. 
es descriptive words frequently require concise definition to avoid confusion — 
‘in the safety 1 assessment. Definitions of terms used in Fig. 5 include: 

. Property— —Characteristic quality, action, state, or condition. 

Behavior—Manner in which something acts, 

4. Safe—Freedom from danger, hazard, or 

Stable—Firmly established; fixed. 
Inspect—View closely; scrutinize. 
8 Measure—Ascertain extent, 

10. Surveillance—Close watch, constant guard, control, management. _ a A 

11. Assess—Analyze critically and judge nature, significance, 


| | 

mrs 

] | 
| | | 
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Evaluate—Express examine and judge worth or condition of. 
14, Interpret—Explain the meaning of. 
he _ A figure depicting as accurately as possible the cross section at the station = 
Z a maximum head difference and a plan of the dam and reservoir constitutes _ 
a useful means of systematically assembling and assessing information on the 
= dam. A lack - of critical | information in this figure constitutes an wars = 
which one must consider in the assessment. 


Execution of a complete safety assessment requires sufficient information 
concerning loads, - parameters, and field conditions to select, represent, and 
analyze a critical section of a facility. Such a clearly-defined situation rarely — 
. _ occurs, and the engineer must therefore accommodate incomplete knowledge — 
and limited data by making “‘conservative’’ choices. eed 
7 __ Engineers it in fields such as aircraft design, nuclear power plant design, and 
\ 4 fire hazard reduction have developed safety programs which utilize a 


in these techniques, particularly i in the design of foundations for nuclear ‘power ‘i 


of alternate designs, and to evaluate the benefits of remedial work. A numerical = 
- assessment of safety promises results superior to the normal deterministic finding ‘ 
= _ A risk analysis attempts to identify all undesirable or abnormal events and =] 
assess the probability and consequences of each One can exclude 
consideration of events of low consequence. Te ult 
To help organize a risk analysis one can arrange events in an event tree 
= or a fault tree. Figures 6 and 7 present a partial combined fault-event tree — 
_ for a waste storage dam. The owner plans to raise the dam crest, thereby = 
_ increasing reservoir storage capacity. Conversations with the owner and an initial 9 
: safety assessment of the dam revealed that the most serious event consists 


of a loss of contaminated fluid stored within the waste storage facility. Loss — - 
of fluid would result in major negative consequences—principally, the flow 
of contaminated water off property and the shutdown of the processing plant _ 
“until containment of fluid was achieved and the water supply replenished. rat © 
following five mechanisms result in loss of fluid: id: 

en ihe pro”ability of tone of f 

Flow occurs along a path of high permeability resulting in flow through ; 
a dam in excess of the capability of the perimeter collection ditch. ae le 
* Excessive rainfall nullifies freeboard and overtops the dam. 


. A sinkhole develops beneath the dam or reservoir. hos 


_ 4. A shear slide occurs through the crest of the dam. wise? 
A pipe develops through the dam by internal ¢ erosion. 


yee 


= 


= 
The occurrence of fombination of various 


arious techniques 
to determine the probability that each 1 mechanism will lead to loss of containment, — 


e.g., determining the probability of loss of containment by a large shear slide — 
in the downstream section involves estimating uncertainties in geometry, por 
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= 2 borer no sudden loss of conteinment 
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te 


essure, strength, and analytical method. Figure 7 shows the : steps involved _ 

obtaining these failure probabilities wiig an event tree. Steps in the probability 
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resides i in the e effectivences core trench to cut off seepage. 
2. Establish peobabiity distributions for permeabilities and strength using 
laboratory data, field data, and the experience of the o> 
—s Obtain pore pressures from flow nets drawn for ranges of permeability 

Compute the minimum factors of safety against a shear § slide for at a range 
a of friction angles for each set of pore pressures. We used the simplified Bishop — 


‘method of slices with | circular arc failure surfaces which lead to loss of — 
‘a . Compute des mean and standard deviation of factor of safety for each > 
“set of pore pressures using the probability distributions of friction angles from 
& 2 and an estimated analytical uncertainty of +6% error in the calculated # . 
. 6. Calculate the probability of failure using the uote f from Step 5 for each . 
x Pe of pore pressures, i.e., F.S. <1 with some probability distribution for F.S., io 
‘for which in this case we tried both normal and lognormal distributions, "with 
s 7, Combine results for each set of pore | pressures and geometry into < one 


= probability of failure using probabilities from Steps1,2,and6. | 


- __ Figure 7 7 gives the probabilities for each of the aforementioned steps, e.g., 
_ from piezometric data, we estimate that an imperfect cutoff exists with a probability _ 
of 0.9. The probability that flow net IV (which we based on sedimented gypsum 
being more than three times as permeable as compacted gypsum, and sand 
being less than 1/3 as permeable as the sedimented gypsum) represents flow 
in Step 2. A probability of 0. 23 exists ‘that friction angles | combined 
_ with pore pressures from flow net IV give a factor of safety less than or — 
: to one. Consequently, the probability of a shear slide due to the combination — 
os low friction angles, pore pressures from flow net IV, and an imperfect 7 Fe 
equals 0.23 x 0.19 x 0.9, or 0.039. Combining all paths to failure in Fig. a 
A gives a probability of failure by a large slide and loss of containment of 0.22. o i 
Internal erosion, or piping of the silty sand foundation stratum underlying 
: the reservoir, causes great concern. Substantial pipes developed in the existing 
4 dam, but loss of containment did not occur, 
The fact that pipes have already developed in the sand foundation of the 
existing dam with some pipes inducing small slumps leads us to suspect | that 
the probability of failure by piping of the enlarged facility must equal or exceed : 
the probability of a shear slide. Consequently, the fault tree shows a probability — 7 
loss of containment by piping of greater than 0.22. 
_ _In summary, the probability of loss of containment in the enpented facility ¥ 
i equals more than 0.44 during its 10-yr design life. Essentially all of the risk — 
ae stems from failure by internal erosion and a shear slide. . With this information = 
the owner readily identifies the main contributors to the risk and assesses the 
a. _As parts of this analysis suggest, we do not know the actual probabilities _ 
as precisely as the number of significant figures Fig. 6 and 7 indicate. The 
figures include the artificial precision to allow ‘Clear illustration of the the 
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_ This paper describes a co mp heaton program to help ensure that a geotechnical ‘We g 
facility meets performance criteria, particularly safety criteria. The program’s 
4 systematic approach rests on the fundamentals of geotechnical engineering—spe- 
cifically, the effective stress principle, the stress path method, Darcy’ s Law, 
4 Successful application of the program yields: tr 
ih portrayal and evaluation of actual performance. ai) 
_ 2. Amassessment of prediction capability, 
3. A reservoir of data on the effectiveness 


of re remedial work. on. 


These activities include: (1) Field i (2) laboratory tests; (3) analyses; 
6 field measurements; and (5) site inspections. The program employs ‘‘average 
2 elements,” a feature of the stress path method (first and second writers, 1979). F 
No program can guarantee ‘‘zero risk’’ of failure. The engineer must still search 
for ‘‘minor geological details” that may cause a facility to malfunction. Further, — a 
the geotechnical engineer must still deal with phenomena such as piping, dam 
_ cracking, and soil disturbance, which the profession only partially understands. 
Believing in the approach of ‘‘engineering certain constructed facilities for 
entire lives,”” have developed—or actually evolved—the program 


of the effor the effort required to to implement it. et ed 
This program, both in concept and execution, benefitted from contributions ; 
ys many — R. V. Whitman and G. B. Baecher assisted i in n the risk assessment 


_ The Borden Company and LAGOVEN, 's. Poy encouraged ae assisted in the 


_ development and application of the program described in this paper. They 


4 on the assessment of dam safety. wade 
_ We acknowledge with appreciation me: assistance we received. Se 
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ELEvATIONCHANGES 


By Dwain K. Butler’ 4 


q One of the primary y observable of | large reservoirs” 
is the deflection of the earth’s surface, presumably | due to the imposed load. _ 

2 For small reservoirs, elevation changes may be immeasurably small, while for > 
large reservoirs, elevation changes are easily observed. Also associated vith 
__ reservoir filling (as well as with subsequent seasonal changes in reservoir level) — 
x are changes in ground-water levels in permeable formations in the vicinity of © 


the 1 reservoir. The objectives of this note are Assess the potential of 


impoundment; (2) determine what “ size”’ reservoirs are re likely to produce elevation =] 
_ changes large enough for detection by microgravity surveys; and (3) orl 
the role of coupled microgravity and surface leveling or ground-water-level . 
monitoring surveys in areas with permeable bedrock. Microgravity survey nets 
in n operation. at Libby I Dam and Reservoir (Montana) and the 


m/s? An change of 10 mm at a given gravity survey 
no associated mass loss or gain in the subsurface, would produce a change 4 
in the vertical component of the gravitational attraction of about 0.003 mgal | 
or 3 wgal. Studies have shown that, with careful repeat gravity measurements 

: vat stations of a survey net with a LaCoste-Romberg Model-D gravimeter, it 
is possible to obtain solutions to the survey net with a precision and accuracy | 
in the estimate of relative gravity values of 1 pgal-2 pgal (11,12,13) with a A. 
high confidence level. The LaCoste-Romberg Modei-D gravimeter has a reading 


"Research Geophysicist, Earthquake Engrg. and Geophysics Div., Geotechnical ‘Lab. oc 
- United States Army Engineer Waterways Experiment Station, Corps of Engrs., P.O. a 
bd _ Note.—Discussion open until August 1, 1981. To extend the closing date one month, 
a written 1 request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on July 30, 1980. ’ 
This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings — 
of the American Society of Civil Engineers, ©ASCE, Vol. 107, No. GT3, March, 1981. 
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of about | ygal. Thus, using a true 

7 as the Model-D, with care, elevation changes or differences of 10 mm can 4 

be detected by repeat surveys of the same gravity survey net at different times. _ 

_ The precision and accuracy achieved in “‘routine’’ microgravity site surveys 
with the Model- D to determine relative gravity variations is 3 pgal-6 pgal 

(2,3, 11,12); thus even in this context, “elevation changes or differences of 20 

mm can be detected by repeat surveys. Observed and predicted (calculated) < 

crustal deflections associated with various reservoirs will now be surveyed in — 


order to estimate the “size” for which microgravity surveys 

— 


CTION Ossenvanons AND Carcura 


Downward crustal deflections have been ‘observed following the. impounding 
on many reservoirs. Elevation changes of 200 mm have been measured at Lake 7 


1. of Calculated Maximum Vertical Deflections 
a 
Reservoir 


= 
Calculated maximum. 
vertical deflections, 


millimeters 


‘Total 
Joad, in kilograms 7s 
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_ Lake Kariba, Rhodesia, the third largest man-made lake in the world, deflections 
4 along a 1 carefully repeated level traverse were as great as 130 mm next to the 
-_ reservoir, and agreed almost exactly with elastic deflection calculations (4,5). — 
The calculations for Lake Kariba predicted a maximum deflection directly beneath 
_ the reservoir of about 250 mm. Similar calculations have been made for fl 
_ other smaller reservoirs (1,7,9,10), and Table 1 summarizes the calculated 
_ maximum vertical deflections and total surface loads represented by the six 
reservoirs; the total ‘surface load ‘represented by atten Reservoir is included © 
_ Using the data in Table 1, the relation of total surface load and maximum 
calculated elastic deflection, or observed deflection, or both, is shown in Fig. 
1. The dashed line indicates the trend of the data. The data for the six reservoirs 
suggest | 
reservoirs have different areal extent, varying water depth ‘distributions, and, 7 
of course, for the field observations, all of the complexities of differing geologies ot 9 - 
are present. This seems to imply that the total load is the key parameter in a 


ves rION 
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determining deflection. Results of measurements of elevation: changen near Oroville 


__ Reservoir are smaller than corresponding values from the calculation (see Fig. 
= 4) . However, the elevation change data for Oroville are complicated by the ; 
occurrence of the 1975 Oroville earthquake (1, 14). b 
_ Due to Elevation Change Alone.—Considering the free-air gravity changes due | 
to the elevation changes, neglecting any other effect at this point, the maxim 
we 
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FIG. —Vertical Deflections, and Calculated, as Function of Total Surface 
a to be expected from reservoir impoundment would be about 75 pgal for an 
_ “elastic-type”’ crustal deflection. However, cases in which an earthquake occurs _ 
in the vicinity of a microgravity survey net may exhibit locally greater elevation 
and gr gravity change (8). The 75-ygal figure is an upper limit from a practical j 
point of view since the corresponding deflection would probably be directly 
under the reservoir. For the case of Lake Kariba, 200 mm represents the maximum : 
calculated deflection which would have been accessible to measurement, or 
about 60 pgal gravity change, which is well within the detectability limits of . 


the Model-D meter with ‘survey could be be detected 
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« 


MARCH 1981 


with standard gravity such as the Worden or LaCoste-Romberg 
- Model- G. On the lower end of the spectrum, a gravity change as small as 
i 3 wgal is detectable with a reasonable degree of confidence with an instrument — 


having microgal precision and accuracy. A 3-pgal gravity change implies a 10-mm 

elevation change (assuming no other interfering effects). From Fig. 1, this implies 

- that microgravity will be of possible usefulness in detecting vertical deflections 

of the crust only for reservoirs representing loads: of the order of 10” kg (1. 1 

x 10° tons or 8.1 x 10° acre-ft) or greater. 

_ Interfering Time Effects.—The assessment thus far has made the rather idealistic — 
a that there are no gi re effects, i.e., that gravity differences 


1. Even though observed deflections agree with elastic deflection anion, J 
: the response time for de deflections probably cannot be ‘‘instantaneous”’ considering 
2. Reservoir filling will not only raise the ground-water level around the 7 
Reservoir but force water into permeable basement rock. 
3. Reservoir level fluctuations will cause ground-water level fluctuations similar 
“4 4. Rainfall percolation into the subsurface can have a significant effect on 
Gravity readings by causing seasonal fluctuations in ground-water 


5. Earth-tide effects must be included in the gravity data analysi 


6. Long-term drift in gravity | 


_ 7, Background noise level. 


‘The ground-water | level changes (items 2-4) cause mass changes which cannot 
be separated from elevation changes in the overall gravity effect wines : 
independent information on either the ground-water level variations or the elevation on 
changes. Item 1, of course, implies that a long-term monitoring program \ _ 
i be involved, e.g., reading a survey net two or three times a year over a several 
year period before and after reservoir filling. The data can be corrected for : 
earth- tide variation by the use of theoretical tide programs © ora base station 
__ The subsurface 1 mass aes represented by ground- -water level fluctuation: 
are seen to be a major problem. One solution is to coordinate precise leveling 
surveys with the gravimetry. While this would not determine elevation changes — 
directly with gravimetry, it would allow the separation of elevation effects from 
mass change effects in the gravity data (15), which would be of considerable 


the possible role played by w water movement and pore pressure changes i in triggering 
~ earthquakes). ‘It would be necessary in this regard to begin the gravity survey — 
several years before impoundment in order to have knowledge of the effects 
On gravity readings caused by normal, seasonal ground-water fluctuations. a“ 
~ Another is to locate the gravity stations on of ‘‘impermeable’’ 


el fluctuations, provided by piezometers and other instruments. Knowledge 


 gource of inter ering ¢ ects. Ouowlng 1S a partia ist Of time-relate 

A 

| 

— 

od 

— 


water level at time of each allow the ‘separation 
_ Of gravity changes due to elevation changes from those due to ground-water _ 
level changes. In general, previous work has indicated 10-pgal—30-ygal changes ae 

in some areas due to seasonal ground-water level fluctuations (11) and changes — 
in the -30-pgal-40- gal range have been noted in repeat surveys over an area 

before and after heavy rains (presumably due to very shallow | ground- water 

concentrations) (unpublished data by Butler). lo, 
_ Survey Procedure.—The survey procedure itself is of crucial importance (3). : 
The same meter should be used over the complete traverse and each traverse a 


and long-term drift. errors (items 6 ‘and. 7) is to use more than one instrument 
Bs s the survey (6); however it is essential that each instrument be used pad 
1 the complete traverse and the data reduced separately prior to combination. | 
The survey grid points should be substantial concrete monuments a 


variation, so that repeat readings at stations are taken at similar points on the 
tidal curve as earlier readings. This synchronization of repeat readings will — 
accomplish two important things: (1) The readings will be at about the same — 
4 on the meter measuring screw to avoid errors from screw irregularities; a 
and (2) soil and rock beneath the stations will be in similar states es with respect. 
to tidal loading. Typical background noise levels (from cultural and natural — 
: sources) are in the 5-.gal-10-.gal range, although during earthquakes, the “‘noise’’ 
4 level can easily exceed the total tidal variation. For a period of several hours 
4] following large earthquakes, obtaining gravity readings in a survey net will be | 
impossible. With the ¢ capacitive readout capability on the Model-D, it is easy a 
to monitor the background n¢ noise level. If required, the | station reading c can be 
obtained as the mean of several readings over a 10-min n period. — a  quekit a 
In spite of the and time- dependent interfi fering eff fects, 
repeated microgravity surveys show promise of being a valuable tool for studying | 
elevation changes, or mass change effects associated with reservoir impounding, _ 
or both. A reservoir representing a total post-impoundment load of the order 
- 10° kg (1.1 x 10° tons or 8.1 x 10° acte- ft) or greater is a candidate for 
‘microgravimetric Studies for determining elevation. changes. Elevation | changes” 
associated with impoundment of Libby Reservoir are probably | between what 
should easily be detectable by microgravimetry and what can be detected only 
by very careful surveys using microgravimetric techniques. 
"Of great importance, from the standpoint of understanding crustal deflection 
and possible earthquake- inducing mechansims, is the need for a study of the © 


this s source to +1 ugal « or ess. Also, readings c can n be synchronized to ‘the cet 


_ reservoir filling and the arrival of pore pressure changes a as a function of depth 
and distance from the reservoir. Considering the nature of the geologic system > 
involved, it is difficult to imagine that the vertical deformations could occur 

“instantaneously” as the -Teservoir load ‘is applied, even though the final 
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__ deformations agree very closely with purely elastic animes calculations. 
_ Repeat microgravity surveys coupled with precise level surveys can play an 
ss... is interesting to note, as a conclusion to this ass assessment, that the length 
of the traverse line is not a limiting factor in assessing the accuracy of a 
microgravity survey. The accuracy of a microgravity survey does not degrade 
 daectly as a function of survey line length. Also, microgravity surveys do not 
depend on on line of sight- -connected survey stations. However, for even the most 
=: se and accurate level survey, the accuracy degrades at least as fast as 
a,L'” in which L = the survey line a in kilometers; and a, depends 


3 pete of this level survey degradation is given in Ref. 14. "Thus, for long 
survey lines, and in rugged terrain where level surveys become more difficult, 

microgravimetry should be considered for elevation change surveys. POT 
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provided by the United States Army Corps of Engineers during the writer’s oe 
—Rerenences 
1. Beck, J. L., ‘“‘Weight-Induced Stresses and the Recent Seismicity at Lake Oroville, 
be California,” Bulletin of the Seismology Society of America, Vol. 66, 1976, pp. 1121-1131. 
2. Butler, D. K., “Assessment of Microgravimetric Site 


3. Butler, D. K., “‘Microgravimetric Techniques for Geotechnical Applications,’’ Miscel- 
 laneous Paper GS- 80-13, United States Army Engineer Waterways Experiment Station, 
Vicksburg, 

4. Gough, D. L., and Gough, W. I., “Stress and Deflection in the Lithosphere near 

Kariba, 1,’ Geophysics Journal ad the Society, ‘Vol. 21, 


‘5. _ Gough, D. Gough, W. “Load- Induced at Lake Kariba,” 

_ Geophysics Journal of the Royal Astronomical Society, Vol.21,1970. = = = “a 

6. Grammel, R. B., Tarman, D. W., Clover, R. C., Leggewie, R. M., Aronstam, P. + " 

i ~ Kroll, R. C., and Eppink, J., " “Precision Gravity Studies at Cerra Prieto—The — 
Second Year,’’ LBL-9551, Lawrence Berkeley Laboratory, Berkeley, Calif., Feb., 2 


am,”” 


et Green, R. Ww. E., “Seismic Activity Observed at the Hendrik Verwoerd Dam 
International Colloquim on Seismic Effects of Reservoir Impounding, The Royal Society, 
+ Guha, S. K., Gosavi, P. D., Agarwal, B. N. P., Padale, J. G., and Manvadi, S. 
oo “Case Histories of Some Artificial Crustal Disturbance,” * Engineering Geology, 
Amsterdam, The Netherlands, Vol. 8, 1974, pp. 59-77. iit Jas, a 

> . Gupta, H. K., and Rastogi, B. K.., ‘Dams and Earthquakes, Elsevier Scientific Publishi 

; 7 10. Lee, T. “A Method for Computing the Deformation of the Crust Caused by the 
tes Filling of Large Lakes * Bulletin of the Seismological Society of America, Vol. 62, 


Lambert, A., and Beaumont, , “Nano Variations in Gravity Due to Seasonal 
Movemeats, Movements, for the Gravitational Deflection of Tectonic 
Movements,”’ Journal of Geophysical Research, Vol. 82, 1977, pp. a9t-S08. ca 


4 
| 
| = 
3 


12. McConnell, R. K., D. Winter, P. J., “An Evaluation of 
be LaCoste- -Romberg Model-D Microgravimeter,”’ presented at the 1974, Seventh Meeting _ 
of the International Gravity Commission, held at Paris, France. 2 a 
_ 13. Neumann, R.., “Gravity Measurement Prospecting Applied to the Location of the 

‘oe Underground Cavities,’ ” Compangnie Generale de Geophysique, Massy, France, 1973. 

14. Savage, J. C., Lisowski, M., Prescott, W. H., and Church, J. P., “‘Geodetic Measure- 
os ments of Deformation Associated with the Oroville, California, Earthquake,” sean! ee 

of Geophysical Research, Vol. 82, 1977, pp. 1667-1671. 
5. Singh, S. K., Mena, M., and Monges, J., ‘Secular Variations of Gravity in Mexico?” 4 
Geophysical Research Letters, Vol. 6, 1979, pp. 469-472, 

‘16. -Westergaard, H. M., and Adkins, A. W., “Deformation of the Earth’s Surface Duc 

_ to Weight of Boulder Reservoir,” Technical Memorandum No. 422, United States — 


of Reclamation, 1934. 


ef mea 


AER, 


| 
wets 


| 


ENGINEERING 


= Note.—This paper is part of the Journal of the Geotechnical Engineering Sas 

“March, 1981. of the American Society of Civil @ASCE, Vol. GT, 


EOE 
16086 
=z 
“3 
| 


JAMAYOL 


Discussions may be submitted on any Proceedings paper or technical note published in in any oe 
Journal or on any paper presented at any Specialty Conference or other meeting, the Proceedings e 
_ of which have been published by ASCE. Discussion of a paper/technical note is open = 
anyone who has significant comments or questions regarding the content of the paper/ technical = 
- note. Discussions are accepted for a period of 4 months following the date of publication — 3 
7 _ of a paper/technical note and they should be sent to the Manager of Technical and Professional 
_ Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period may 
be extended by a written request from a discusser. 
- The original and three copies of the Discussion should be submitted on 8- 1/2-in. (220-mm) . 
= by 11-in. (280-mm) white bond paper, typed double-spaced with wide margins. The length of — 
_a Discussion is restricted to two Journal pages (about four typewritten double-spaced pages 
of manuscript including figures and tables); the editors will delete matter extraneous to the 
~ _ Subject under discussion. If a Discussion is over two pages long it will be returned for shortening. — 
All Discussions will be reviewed by the editors and the Division’s or Council’s Pubieations 7 
‘ Committees. In some cases, Discussions will be returned to discussers for rewriting, or coll Sy 4 
_may be encouraged to submit a paper or technical note rather than a Discussion. amen 
Standards for Discussions are the same as those for Proceedings Papers. A Discussion is 
subject to rejection if it contains matter readily found elsewhere, advocates special interests, 
is carelessly prepared, controverts established fact, is purely speculative, introduces personalities, 
or is foreign to the purposes of the Society. All Discussions should be written in the third 
- person, and the discusser should use the term ‘“‘the writer’’ when referring o himself. The 
author of the original paper / technical note is referred to as “‘the author. 
oir have a specific format. The title of the original paper/technical note appears 
at the top of the first page with a superscript that corresponds to a footnote indicating the 
month, year, author(s), and number of the original paper/technical note. The discusser’s full 
- mame should be indicated below the title (see Discussions herein as an example) together with 
‘page of the manuscript, along with the Proceedings paper number of the original paper / technical 
note, the date and name of the Journal in which it appeared, and the original author’s name. | 
= that the discusser’s identification footnote should follow consecutively from the original _ 
paper/technical note. If the paper/technical note under discussion contained footnote numbers | 
1 and 2, the first Discussion would begin with footnote 3, and a Discussions would 


7 it is suggested that potential discussers request a copy of the ASCE Authors’ Guide to 
the Publications of ASCE for more detailed information on preparation and atuietien of — 
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STRENGTH AND DEFORMABILITY OF Hicmy 4 
Discussion by by J. Laginha Serafim,? F. ASCE 
pay paper by Raphael and Goodman is an excellent contributio nt to 


knowledge of fundamental properties of rock masses with frequent joints and 5 
cracks, viz. deformability and shear strength, 


“core idea of of building up specimens of a very ry jointed rock mass from fitting a 


‘el a post- aA mete were determined. The writer has insisted in various ae 
7 Le occasions that triaxial tests on cylinders of materials with oriented weaknesses 
_ cannot be interpreted correctly by Mohr’s theory (including an — 
- a presentation at the 1962 Rock Mechanics Ss Meeting, held at Salzburg, Austria); 


et 


neither are such tests valid if the surfaces of minimum strength cut through 
the end steel plates. In fact, for layered, aniso-tropic or oriented fissured solids, " 
=. _ the surfaces of rupture will go through the surfaces of minimum strength and 

j : _ not along the direction indicated by the tangent of the supposed envelope and 
_ the Mohr’s circle. On the other side, when the planes of rupture have to cut 
through | the steel loading plates, the modulus of elasticity and maximum strength 
will be affected by them (18), unless some special precautions are > taken (use 
i of brush loading plates, etc.). The post-peak strength indicated by the tests 
pitied by Fig. 3 and summarised in Fig. 4 of the original paper possibly | has 
not allowed to estimate the residual strength. This because of the previous 
- : reason and possibly because the authors were not sure that the final siding 
occurred along joints. Slip took place, as they admit, in new fracture surfaces. 
They certainly were careful in ‘designating this minimum strength the ‘“‘post- peak = 4 
strength’ and not the residual strength of the rock mass, which corresponds _ 
: to the final frictional strength along the joints. For this reason it is rather surprising = 
that the authors have not conducted friction tests in joints. Nor have they 
classified the degree of roughness of the joint surfaces exposed by the cores, 


from which a reasonable value of the peak strength could also possibly 


_ Ina recent paper by Singh (19), the residual strength in a fracture « developed 
‘ - during the triaxial tests of an intact rock specimen was also determined by i 
_ successive increases of the loads after the initial rupture, 
; _ Borehole jack tests (using the Goodman jack) were made for determining 
ss the deformability « of the rock mass. The writer believes that the Problem | = | 
- very close adjustment of the jack surface to the boreholes surfaces was 


1979, by Jerome M. Raphael and Richard E. Goodman (Proc. Paper 14988). 
— of Civ. Engrg., Coimbra Univ., Portugal; also President COBA-International — 
Consultants, Av Marqués de Tomar 96, 1000 Listen 1, Vin) 
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an imperfect can about erroneous results. 
This study was certainly a good opportunity to compare the deformability er 
a of the rock mass obtained by various methods. Unfortunately, large jack tests — 


_ were not performed at the site to compare results with the other determinations on 


rock mass, for instance, according to Bieniawski’s (17) geomechanics dade 

_ tion system, RMR, and introduce new points and | curves in his diagrams. Neither = 

the RQD, nor the percentage of core recovery, nor the complete statistics of 
number of joints per | unit length of the hole and the opening of the ~ joints iy 
were in indicated, which would make this a very good opportunity to add points a 

to the correlation between the RMR and the in situ modulus of elasticity. Maybe - 

the authors can still make the determinations and, in this case, they will give 

an additional value to their important paper. It must be mentioned that the — 
7 possibility of transforming results of dilatometer or borehole jack tests into 5% 

results of in situ large flat jacks is an an important a: aspect | of the ooo lll 

& 


17. Bieniawski, T., “‘The Geomechanics Classification in Rock Engineering Applica- 
i, Rcisdite of the 4th International Congress on Rock Mechanics, Vol. Il, 
18. Serafim, J. L., and Guerreiro, M. Q., “Tests and Observed Behaviour of Rock 
Foundations of Dams,”’’ Proceedings of the 3rd on 

19, Singh, D. P., “A Study of Frictional Properties of Rock,’ * Proceedings of the 4th 
International Congress on Rock Mechanics, Vol. 1, 1979, pp. 301-305. 


Closure by Jerome M. Raphael," F. ASCE ‘on 
and Richard E. Goodman,° M. ASCE 
a nm e authors wish to thank J. L. Serafim for his complimentary and constructive _ 
discussion. The question of multi-stage testing on rupture surfaces in triaxial 
compression experiments is indeed complicated because the failure planes may : 
tend to intersect the end caps. In the experiments in question, so many possible 
failure surfaces exist that it is reasonable to assume that a failure — 
orientation can be found close to the orientation of least strength in a position 
_ within the specimen that does not necessitate intersection with the end caps. 
This would not be true in a general multistage test but is assumed to be true 
: this specific case in which so many incipient fracture directions coexist “tin 
Large block tests were not conducted at the site because money was not : 
available and the results were needed quickly. Direct shear tests on joints were were — 
a. “Prof. of Civ. Engrg., Univ. of California at Berkeley, Berkeley, Calif. enn 
 *Prof. of Geological Engrg., Div. of Engrg. Dept., Univ. of California at Berkeley, 


—, 


not run because samples of individual joints | of ‘satisfactory quality | for s 
From estimates of the ‘rock mass properties, the Bieniawski Rock Mass Ratin g 
= (20) is determined to be from 30-35, qualifying the rock mass as ‘‘Poor Rock,’ = 
a result we already appreciated before starting the tests. The euvenpetiieg ; 
rock mass deformability cannot be obtained from Bieniawski’ s equation: E in 
Gigapascals) equal to twice the rock mass rating minus ; 100. The reason is 
simply that this linear fit to Bieniawski’s data points intersects zero pane Al ay 
at a rock mass rating of 50 while the applicable rock mass rating is less than 7 
i° 50. If a curved fit to Bieniawski’s data points is provided, as shown in Goodman y 
— QD), the point with RMR equal to 33 and E equal to 6 GPa, which corresponds e 
to the results in our paper, will be seen to plot close to the extrapolated curve. 
Much more testing in the low RMR region will be required before such correlations i 


Aprenvix.—Rererences by Bq. 4 bev m 
20. Bieniawski, Z. T., ‘‘Determining Rock Mass from Case 

Histories,” International Journal of Rock Mechanics and Mining Sciences, Vol. 15, 7 


21 ‘Goodman, R . , Introduction ‘to Rock pieslenten John Wiley & Sons, New York, 


Discussion by Yoshiharu Moriwaki, A. M. ASCE. 


and Maurice S. Ss. Power,’ M. _ASCE 


way of evaluating the pore water pressures 
' =e under undrained cyclic loading conditions. The assessment of cyclical- 

. 2 ly-induced pore water pressures is becoming increasingly important because 
of recent emphasis on pore water pressure generation / dissipation and —e 
quake behavior associated with earth structures (12). The writers have also 
used a similar approach in a recent project to evaluate the results of excellent 

: _ undrained cyclic triaxial tests on modeled Oroville gravel by Banerjee, et al. 
(9). Asummary of the results and some comments on the ae of boundary 


_ January, 1980, by Edward G. Prater (Proc. Paper 15100), ts 
*Geotechnical Consultant, Civ. Systems Inc., 2450 Washington | Ave., 
Woodward-Clyde Consultants, San n Francisco, Calif tof) 
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DISCUSSION 


, 2, and 3, but ‘suited for the results of cyclic 


max 


are shown in Fig. 4; also shown in Fig. 4 are the curves computed using Eq. 
_ 4combined with $’ values reported by Marachi, et al. (11), for the same material 
under various initial confining pressures. It is clear from Fig. 4 that Eq. 4 
(which is s based on simple Mohr- Coulomb co considerations) i is sufficiently accurate 
to describe the trends in the observed effects of K. and 
_ The writers feel that Eq. 4 is probably useful in assessing an approximate _ 
value of the maximum residual pore water pressure that can develop in an 
_ undrained cyclic triaxial text specimen. The emphasis on the residual pore water — 7 
p 4 pressure is important because the transient Pore water pressure can easily be 
higher than that indicated by Eq. 4 because an isotropic stress state can exist 
in a specimen during cyclic loading. The emphasis on the triaxial test is important 
7 because in a conventional triaxial test, the specimen boundary is kept at constant 
ambient o, ando, throughout the test. 
= However, in those tests or situations where the specimen or soil boundary — 
is kept from permanent deformations, say in ‘the lateral direction, Eq. 4 is 
no longer applicable since the total lateral stress may change during cyclic loading. 
This is shown in Fig. 5 which shows results from undrained cyclic torsional 
— tests on hollow cylindrical pense by Yasuda (10). Fig. 5(a@) shows a case 
_ where the total lateral stress, o,, was kept constant but the lateral boundary 
ne was allowed to deform permanently during cyclic loading (similar to cyclic triaxial 
tests); the cyclically-induced pore water pressure is limited and probably controlled 
by a relationship somewhat similar to Eqs. 1-4. Fig. 5(b) shows a case where 
permanent lateral displacement was not allowed; o, increases and eventually 
becomes equal to oa, and the pore water Pressure 2 also increases and becomes — - 
la For horizontal soil deposits where the lateral boundary condition is probably J 
closer to that of no permanent lateral displacement because of lateral symmetry, 
it is likely that the pore water pressure can increase to the value equal to 
o”, because the total lateral stress can change during cyclic loading. Therefore, 
the author’s Eq. 3 is considered to be not appropriate for determining the maximum 
residual pore water pressure that can develop in horizontal soil deposits. ? 
‘4 .—Rerenences cals Taw stad at of po mvt aa. 


My Banerjee, N. G., Seed, H. B., and Chan, Cc. K “Cyclic Behaviour of Dense 
Coarse-Grained Materials in Relation to the Seismic Stability of Dams, Earthquake 
“ie Research Center Report No. UCB/ EERC-79/ 13, University of California, , 
™ ‘Ishihara, K., Determination of Liquefaction Parameters,’’ Dynamical Methods 
i) ad in Soil and Rock Mechanics, Proceedings, Karlsruhe, G. Gudehus, ed., Vol. 2 Plastic 
aa % and Long-term Effects in Soils, 1978, pp. 189-193. 


min Marachi, D. N., Chan, C. K., Seed, H. B., and Dunca an, J. M., pen 
_ Deformation Characteristics of Rockfill Materials,”’ Report No. TE ‘69-5, University 4 
Of California, to State of California Department of Water Resources, Berkeley, Calif. = l 


equation simila 

q 1 

| 
j 


= Seed, H. B., ‘*Considerations in the Earthquake- Resistant Design of Earth and Rockfill 


iy The sole discussion contribution is from Moriwaki and Power. The main point 

& concerns the importance of sample boundary conditions in relation to 
horizontal soil deposits. For this field condition lateral restraint under strong 44 
shearing action requires that the associated inelastic strains occur in a vertical 

7 direction only. In the triaxial test, the sample is subjected to constant stress 
boundary conditions. The sample deforms as a cylinder without constraining &. 

any of the principal strains to maintain a zero value. Thus quite different stress 
4 paths are involved in the two cases. This is an important consideration and ale 
the writer is grateful to the discussers for introducing it. They | also refer | in a 
“ this connection to the cyclic torsional tests on hollow cylindrical samples of by 

Yasuda reported by Ishihara (10). The writer wishes to consider these tests 
in greater detail, but beforehand some more — porns in the discussion 


is not modified | so > the use of coy] similar is inappropriate. In ‘Ref. 15) 
writer has also expressed the relation in terms of wollen ao 
‘The writer had also observed (p. 115) that the induced pore pressures can 
transiently exceed the residual value, | even in the case Of dense gravels. This 


writer, however, would like: to stress again ‘the i importance of considering ‘changes 
_ in the grain matrix. As indicated in Fig. 1 if the structure collapses, the initial — 
anisotropic stress state is lost and the effective stress path rapidly moves to 
3 the origin in accordance with a truly liquid state. In the paragraph concerned 
1 with loose horizontal soil deposits (p. 113) the aim was not to find the residual — 
- pore pressure (although unfortunately the author retained the same designation, 
= ), but simply to point out that the failure envelope may be reached during | ~ 
a cycle before Au reaches the value p’. Thus the discussers final remark is 
“not directly relevant to the writer’ "s intention. This, however, brings out the 4 
point as to when it is appropriate, generally, to use Eq. 1. The condition is 2 
J obvious: if the state of total stress after the seismic disturbance is the same z 
4 


as before. Not only in horizontal deposits, but also in any earth structure like 
adam, if the stresses are redistributed the induced pore pressure will change 
accordingly. The question remains to be investigated, however, as to what extent 
4 the static stresses are changed for loose and dense deposits, respectively. 
It is encouraging that more data ; are becoming available for gravel samples, 


accurate description of the friction angle for the gravel tested. The strength 


4 see Fig. 4. The writer would like to take this opportunity — to present a ‘amore 


envelope is, in fact, emnyes with the following empirical relation for ies = 2.28 “a 
kg/cm’ (Dr = 85% (Dr = 85%) 


“Research Engr., Federal Inst. of Tech. Zurich, Inst. of Foundetion Ears and Soll 
Mechanics, ETH-Honggerg, 8093 Zurich-Switzerland. 
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degrees in the range 04, < 8 kg/cm’. mor 
Fer stress paths involving large increments of stress the dependence of ” 
" - a3. should be recognized. For simplicity’s sake, ce, Eq. 1 assumes that the 
strength relationship is a straight line. sh 
Finally, the writer would like to address himself to the torsional | og on 
1 hollow cylindrical samples. The discussers state that for the test without lateral 
confinement [Fig. 5(a)] the induced pore pressure is ‘‘probably controlled 
a relationship somewhat similar to Eqs. 1-4.’ To check this it is, of course, 
necessary to know the complete state of sues in the sample. As pointed out — 
by Vaid and Finn (17) nonuniform stresses can be imposed in this type of © 
test during the consolidation phase. However, assuming uniform. stresses for 


am anisotropic stress state circumferential stress % = For 


-->-0 
(a) 


= 0.43 k g/cm? 


FIG. 6.—Torsional Test (Compare with Fig. 5) with Stress-Path Representation meena 


_ Evaluation of M using Residual Pore Pressure at Lean Condition » (Zero Cyclic Stress) 
_ Showing Influence of Total Stress Redistribution _ 


 Yasuda’s test this would give p, = 0.667 kg/cm’, g, = 0.5 kg/cm? as shown 
in Fig. 6. From Fig. 5(a) the induced pore pressure equals approx. 0.43 kg/cm’. 
Assuming the failure envelope is reached after the period of shaking the line 
OA may be drawn to give M = 2.14, in which ’ = = sin [3 M/(6 + MO = 32 
the medium dense Fuji river sand (Dr = 55%) this value is obviously 
_ far too high. In another investigation (14) using the triaxial apparatus, an average 
_ value of M for compression and extension failure is around 1.25 at this value — 
_ of Dr. The only way for this value to be achieved is to allow the circumferential - 
_ stress to reach a ‘‘full plasticity’’ condition with o, = a, the total stress p = 0.83 


kg/cm?’ moving to point C (Fig. 6 giving p= = 0. 40 40 kg /cm?. ’ The corresponding 


"DISCUSSION 
| 
from 
hicdty higher 
as 


1888 
of M = 1.25, (line OB), i.e., = = 31°. for reasons o of ‘static : 
- equilibrium, if the internal one external radial pressures are equal, the stress 
component o, must equal o,, (13). This suggests that the stress conditions in ~ 
the specimen are far from being uniform. In the case of the second test [Fig. 
-5(b)] radial strain is prevented | and the he shearing action cases the state of stress _ 
to be transformed to a hydrostatic one (point D). The initial deviator stress — 
is eliminated s0 that the specimen is now sheared in the reverse direction - 
__‘There is an obvious need to more closely examine the state of stress and _ 
"i influence of end restraint in torsional tests on hollow cylinders. The proposed a 
“ equation (Eq. 1) should also be verified for other tests, e. eB.» the true triaxial = 
test. 


Ishibashi, ond Sherif, M. A., ‘‘Soil Liquefaction by Torsional Simple Device,” 
Journal of the Geotechnical Engineering Division, ASCE, Vol. 100, No. GT8, Proc. : 
= Paper 10752, Aug., 1974, pp. 871-888. 

14. Ishihara, K., and Okada, S., “Yielding of Overconsolidated Sand and Liquefaction _ 

Model under Cyclic Stresses,’ Soils and Foundations, Vol. 18, No. Mar., 1978, 
1S. Studer, J. , Zingg, N., and Prater, E. G., ‘‘Investigation of the Cyclic Stress-Strain - 
: _ Characteristics of Gravel Material,’’ Proceedings of the 7th World Conference on — | 
Earthquake Engineering, Vol. 3, Sept. 1980, pp. 355-361. 

16. Townsend, F. C., and Mulilis, J. P., “‘Laboratory Strength of Sands under Static 


a and Cyclic Loadings,”’ Research Report S-76-2, No. 6, United States Army Engineer 


a Vaid, Y. P., and Liam Finn, W. B., ‘‘Static Shear and Liquefaction Potential,”’ Journal 
of the Geotechnical Engineering Division, ASCE, Vol. 105, No. GT10, Proc. wpe 
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OF SAN Francisco Bay Mup 


Discussion sion by T Thomas Cc Anderson’ and Robert G. Lukas,* — ASCE 
The work of Clough and Denby is to be commended for the advancement 


of the state of the art for techniques and use of the self-boring pressuremeter 


_ in soft clay. The writers wish to discuss the use of s,/p ratios from the laboratory 2 
UU and vane tests and the comparisons with similar ratios derived son} CU a 


ie “January, 1980, by G. Wayne Clough and Gordon M. Denby (Proc. Paper 15148). _ cane” 
*Sr. Geotechnical Engr., Engrs. Testing Labs., Inc., P.O. Box 21387, Phoenix, “Ariz. 
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. ‘Since the soft clay at both : pveat is slightly overconsolidated, it would ap appear a 
_ to be more appropriate to relate the undrained shear strength (s,) from the — 

laboratory UU and vane tests to the preconsolidation pressure (p.), instead — 

of the present effective overburden pressure (p). This is especially applicable 

for comparisons of this strength ratio with the results of CU triaxial tests, 


= 


situ u stress levels and, thus, were on the \ virgin consolidation curve prior 
‘In 1957, Skempton established the following empirical correlation for normally 
in which PI = the plasticity index, as a percentage. Anderson and Lukas ( > 
showed that this correlation is also valid for ‘slightly to overconsoli-_ 


developed for the sites from the laboratory and fieldtests. 

_ For Site I, the indicated s,/p, ratio from Fig. 3 is 0.275 + 0.01 for the entire 
| depth range of Bay Mud, including the desiccated zone above 5.5 m. The average . 
g for the deposit is 45%, which results in an s,/p ratio of 0.28 from Eq. 7 


. These two ratios are in excellent agreement, ‘but considerably lower than . ; 
s,,/p tatio of 0.35 from the CU triaxial tests. __ | 
mr For Site II, the undrained shear strength values show that the strength i increases 

with depth. The indicated s,/p. ratios from laboratory UU and vane tests are 
- - roughly 0.27 and 0.34, respectively. The average PI for the upper and lower | 
depth zones is about 35%, although no PI values are given for the higher water 
f content, middle zone. For a PI of 35%, as s,/p ratio of 0.24 is computed from 
: «Ea. 1. However, the plasticity index of the middle zone is undoubtedly higher 
Ss probably would result in an average, computed s,/p ratio similar to Site 
I. The indicated strength ratio from the UU tests and the computed ratio are 
( : thus in apparent accord, but much lower than the s,/p ratio of 0.36 to 0.38 — 
‘< from the CU triaxial tests and also the ratio determined from the laboratory — 
The real question is s which of these strength ratios most accurately predicts 
7 the field strength ratio. This can probably only be evaluated by a test embankment, a 
_ which loads the soil to failure, and then back-calculating the field strength by __ 
refined stability analyses. Data by Duncan and Buchignani (16) for a failure 
of an underwater slope within normally consolidated Bay Mud i indicated that el 
‘ undrained strengths from UU tests on trimmed 1.4-in. diam specimens (s, /P. 
ratio of 0.30 + 0.01) overestimated the actual field strength ratio of 0.24 to _ 
0.25, while the strengths from UU tests on untrimmed 2.8-in. diam samples 
, /p. tatio of 0.25 + 0.01) would have closely predicted the field strength ratio. — 
_ Anisotropically consolidated- undrained triaxial test, which resulted in a s, ip a 
ength by approx 40%. — 
For an n average PL of 29%, | as, s,/p ratio of 0.22 is s computed from Eq. 1 and 
“this ratio would have slightly underestimated the strength in the field. ai 
_ Ladd (18,19) and Wu, et al. (21), presented data from case studies of = 
_ embankments over soft ground showing that direct simple shear (DSS) tests. 
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resulted in a closer prediction of the field : i ad Pp atic than CU triaxial tests, 
which overestimated the strength ratio by 25% or more. Are DSS test results z 
available for the Bay Mud at these two sites? Gardner (17) indicated an undrained - , 
_ Strength ratio of 0.25 for a DSS test on San Francisco Bay Mud with a PI ig 
of 45%, while a strength ratio of 0.35 was determined from a triaxial compression 
test. This DSS test strength ratio is in close accord with the s, /p ratios computed 
fom Eq. 1 and also the s,/p. ratios forthe twosites. 
_ Based on the foregoing comments, the writers conclude that the strength 

“ratios p predicted by Eq. | and the DSS test are in reasonable agreement with 
_ the indicated s,/p. ratios for the two sites. These three undrained strength i 

ratios are felt to be most indicative of the actual field strength ratio. As a — 

result, the writers believe that the strength ratios determined from the CU triaxial 
tests considerably overestimate the field strength. Since the strength values 

_ derived from the pressuremeter tests are higher than the CU triaxial tests results, = 


the pressuremeter results also appear to greatly overstate the in- 1 strength. 


‘15. Anderson, T. C., and Lukas, R. G., “‘Preconsolidation Pressure Predicted Using s,/p __ 
. q aa Ratio,”” Proceedings of ASTM { STP Sympesien on Shear Strength of Soils, June, — 


16. Duncan, J. M., and Buchignani, A. L. , “Failure: of Underwater Slope i in San n Francisco é 
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+: SM9, Proc. Paper 10019, Sept., 1973, pp. 687-703. 
17. Gardner, W. S., ‘‘Soil Property Characterization in Geotechnical Engineering Practice,” _ 
Association of ‘Drilled Shaft Contractors Magazine, Jan., 1980, p. 28-72. 
& Ladd, C. C., “‘Test Embankment on Sensitive Clay,” Performance of Earth and 
Earth- -Supported Structures, ASCE, Vol. 1, Part 1, 1972, pp. 101-128. 
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improved. It is sufficiently different to warrant separate 
__ A special protective sheath for the membrane has been developed. It is a 
composite of stainless steel and rubber that is axially stiff to resist friction 
forces on insertion and withdrawal but radially, it is very ‘compliant. From 
the appearance of early models of this device it is known as a Chinese Lantern. 
It consists of a series of narrow, thin, stainless steel strips which butt together | 
around the periphery of the instrument. These strips are partly bonded on their — 
| inside to a thin rubber sleeve which serves only to exclude dirt from getting Ps 
an _ between the strips during the expansion phase. The sleeve seals no pressure 7 
Cambridge Inst., Little Eversden, CB3 7HE, England. 
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anda adds ‘only a a amount t to the membrane 
_ These bonded lanterns have proved very reliable and easy to use and they 
_ provide such good protection that the membrane puncturing described by Clough % 
and Denby is no longer a problem. One membrane now normally lasts for 
4 a Concurrent with the development of the Chinese Lantern, a membrane more > 
resistant to abrasion has been developed. It is made of a polyurethane polymer, _ 
s Du Pont Adiprene L.42, by injection moulding techniques. Examples of this — 
membrane only twice as thick as that used by Clough and Denby have proved _ 
_ capable of use without Chinese Lantern protection in stiff clay of 700 kN m - 
5 undrained shear strength. att gamer A Yee us ll 
q The San Francisco Bay Mud is insensitive to effects due to the rate of loading. — : 
_ However the fact remains that many materials are rate-sensitive. A unit has 
been developed to carry out tests at a constant rate of strain, completely = 
- automatically. A soft-start control system is incorporated which, by limiting 
- the rate of pressure tise to a preset maximum, permits the point at which — mf 
the membrane lifts off the instrument and, thus, the in-situ stress to be accurately 
zi discovered. Ten rates of strain are provided both up and down, together with 
the ability to hold the strain constant at any value indefinitely. . Thus not — - 
can unloading and reloading loops be carried out, giving a reload modulus, 
f but also the changes in stress caused by a given expansion can be monitored — 


as they take place by observing the pressure being supplied to maintain a set - 
expansion. Carrying out the whole test automatically also frees the operator _ 
_ monitoring the progress of the test and also frees the test from the effects 
Recent work on in situ stresses behind a retaining wall has shown that the 
ability to to measure the radial strain at each individual strain arm position is a 
very important. Previously only the total strain registered at the three arms» 
_ together was read out. The smooth curve of total strain versus total pressure — 
“can conceal distinct changes in the slope of the pressure-strain graph if the — 
- changes occur at different values at each strain arm. The modification to read 
individual arms will, with care used when Laateried the instrument, provide 
oa Finally, Clough and Denby mention that advancing of the instrument is _ 
7 trial and error. Recent work by, among others, Clark (2) unpublished communica- 
tion, has shown that it is possible to be a great deal more certain in the initial — 
- setting of the cutter. This, coupled with changes in the cutter geometry, have 
_ meant that it is no longer a procedure of trial and error to determine the correct 
cutter setting. It is also possible, even in stiff clay, to attain driving rates up 
Closure by G. Wayne Clough’ and Gordon M. Denby 


“The w writers appreciate the interest asia in the paper by Dalton, , and Anderson 
; = Lukas through their thorough discussions. Dalton provides a review of { 
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"recent improvements w which have been made to the Cambridge type e self-boring 
_ pressuremeter. This is very appropriate since the device used by the writers 
y is one of the early models. Clearly the reliability of the probe and the sensitivity — 
of the system have been enhanced. However, even with the new modifications, _ 
_ potential users should be aware that problems will still occur in the case of 


conditions such ; as encountered by the authors at China Basin (Site Il). Here, = 
shell layers existed in the clay which clogged the probe, and likely disturbed 
the clay as shell fragments were pushed or dragged by the edges of the cutting 
shoe. Further, ground currents interfered with the recording system. Careful 
_ consideration needs to be given to potential problems like these before embarking a 
on a self-boring pressuremeter test program. 
_ Anderson and Lukas are concerned with defining an undrained ‘‘field’’ strength, 
and suggest an empirical equation for this purpose based on plasticity index 
and preconsolidation load. They conclude that it is more desirable to define = 
an undrained ‘‘field’’ strength by their equation or undrained direct simple shear 
tests than by CU triaxial compression or self-boring pressuremeter 4 
writers ai are in in agreement that it would be nice to have a simple empirical | or s 


this implies that there is one strength for a soil regardless of degree of anisotropy, y 
stress path, or importance of creep effects. In fact, this is not the case, | and 
the comparisons made by Anderson and Lukas between results from different 
types of tests ignore this. For example, it is well known that a CU triaxial 7 
“4 test typically yields a higher undrained shear strength than does a direct simple © 


"shear test due to effects of anisotropy. Of course, this does not invalidate 
that CU test; rather, the results from this test provide one point on a strength © 4 
_ curve, which, for design purposes, should be delineated in terms of stress path, __ 
anisotropy, and creep effects. Data obtained by the self-boring pressuremeter _ 
should be viewed in the same e light. P1 Prevost (11,12) has provided useful theoretical — 
4 findings i in this regard. ent wo wilt 
a In their discussion, Anderson and Lukas utilize data from the writer’s cle 
as well as that from the paper by Duncan and Buchignani (16) to develop an 
estimate of the field shear strength of San Francisco Bay Mud. The writers _ 
- feel that the approach | used in this exercise | is not ot soundly basec based, for a number — 


ua 1. Data from UU triaxial tests and laboratory vane tests are used to justify — 
the results obtained by the discussors’ empirical equation without question as 
to whether the laboratory tests are appropriate to define a field strength. eat 
2. Data from Bay Mud sites are compared without due consideration given 
_ that the minerology and sensitivity of Bay Mud varies $ from site to site and 


such as creep strength loss or strong peak-residual stress-strain response, which 


affect the field strength, while others arenot. = = 
_ 3. There is no acknowledgement that Duncan and Buchignani (16) had to’ 
reduce their UU test results by 30% for sustained loading effects to obtain — 
a field strength value consistent with the observed field behavior. Making this 


allowance leads to strength value below that obtained by the Anderson and 


= sands to o high plasticity clays. Thus, some Bay Muds are subject to factors 
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Francisco » Bay Mud. To the writers’ knowledge, this information is 
available. However, Duncan (4) has reported a series of carefully performed | 

_ active and passive loading tests which yielded S,/p ratios of 0.37 and 0.28 
for the Bay Mud at Hamilton Air Force Base. Since the DSS test commonly 

_ yields a strength about halfway between that defined by active and compression ~ 
loading, the DSS strength could be expected to be about 0.33. The point made 

¥ _by Anderson and Lukas concerning a need for DSS data however, is well placed, a 
and the ser senior writer of the paper intends to undertake these tests in the near ¥ 
"future. Jon to time reguited wauk tank and che 

s In sum, the writers remain of the opinion that the self-boring pressuremeter 

offers promise as a means to define the properties of a soft clay soil. In particular, — 
q should be reiterated that the device is useful because it not only can determine a 
_ the shear strength of a soil, but also the stress-strain and pore pressure response 
and the lateral stress in the ground. Additional research is needed in order 

' to more fully understand the effects of important key parameters such as rate : 

of loading, insertion technique, and to che megnds 
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Discussion by R Robert E. Bigham,” M. ASCE sce ashe 


authors presented a fascinating of failure of a large tank, 


writer controlled loading of several tanks on soft and 
was initially suspicious of the increase in maximum settlement at about a 60.7-ft_ 
_ (18.5-m) load, and decided to make his own analysis of the data of Fig. 4. 
Bell later provided a larger chart similar to Fig. 4 showing actual data points 
upon request by the writer. Analysis of those data yielded the same conclusion 
as was reached based on the smaller chart, and this discussion _Tefers ‘only | 
to ‘the data shown on Fig. 4, as extracted by the writer. Any error in the 
writer’s extraction of data is solely his responsibility, 
_ The writer agrees with the authors that there was no soil bearing capacity t 
failure at Tank T-270, either during loading or later. 
_ ‘The writer analyzes data from tank loading differently than as sometimes =) 
seen in the literature, and the method may be of interest to the authors and 
so “a sFebruary, 1980, by Roy A. Bell and Jun Iwakiri (Proc Paper 15219). el LS 
Engr., , Spencer J. Buchanan & Assoc., Inc., , Consulting Engrs., , P.O. Box 672, 
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writer’s method is based on the theory that, so long | as loads do not 
approach the then current ultimate bearing capacity too closely, the change ; 
in settlement within a standard time after a change in load will bear a constant — 7 
proportion to the change in the load, and that the postulated constant proportion 
is characteristic of the combination of site and structure. emt 
_ With small changes in load, this can be the principal “control for loading 
a tank, and a chart of load versus settlement is all that is required to a 
it to loads following the initial load. The initial load can be selected by a simple 
estimate of ultimate bearing capacity and a large, e.g., four or five, safety — 
oa eet of water (meters “Ae 
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FIG. 16.—Load Settlement Chart Tank 7-270 
vA disproportionately large (based on history at the particular site) increase > 
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Settlement in millimeters 
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an indication that loading i is approaching ‘ultimate bearing capacity n more closely 
than under the previous load. With small, e.g., 1-ft min, up to perhaps 2% 
or 3% of existing load, changes in load, this is a warning and not a bearing © 
"capacity failure. It can even be observed as a trend developing during several 


_ Fig. 16 shows the writer’s analysis of data extracted from Fig. 4. The 4 
points are ; the ‘Settlement for the last day | of a particular load and the settlement 


at constant load. It can be seen that the slopes of the segments — 


| 

by the chart or other evidence. When the chart gives a warning, loading : 
-_ stopped, and a wait for consolidation ensues. At such a time, a chart of time » 


ibe sand in the wc 

a end points are all acs the same, even for ‘the last i increase, and even though _ 


the 60.7-ft (18.5-m) loading changed the load from less to more than the 
_ ‘To the writer’s way of thinking, this goes further to show that there was no 

‘ bearing capacity failure during loading or later at Tank T-270 

__ This method is equally applicable to cases of large settlement, as at Tank 
_ J-270, or small settlement, or small tanks, and has been a valuable aid to the 
writer’s judgment. It does require many (sometimes twice daily) measurements 
: «of settlement. If changes in settlement are as large as 0.01 ft (3 mm), it is 


possible to rig cheap reliable water levels that reduce the time required Zz 
asiit by Roy Bell,“ 
pia 4 The discussion by Bigham presents still another empirical method of evaluating — 
7 tank settlement data that could be a useful supplement to other methods of 
~ tank settlement data interpretation. The fact that it emphasizes conditions during — 
periods of adding water load could provide an early warning signal of impending cf 
by but since it ignors rates of settlement during periods of constant | 
load, other means of settlement evaluation should continue to be used in addition 
\ Errata. —The following corrections should be made to the original paper: 


Page 167, Fig. 13, the formula for ©: : Should read ot 


; 8 = arctan — instead of 6 = arctan 


> z, sin a, 


Page | 167, Fig. 4, the formu formula fi for 8’ ¥: '; Should read 8’ = 


instead of 8’ = 


$1.54 =1273 
014 = 0.13 minstead of n= = 


Further: corrections to ) Figs. 8 13, and 14 can be obtained from the first author. — 


“Sr. Engr., Dames & Moore, 500 Sansome St., San Francisco, Calif. 941 =e 
*Mgr., Civ. Engrg. Dept., and Construction Co., ‘Ltd. Yokohama, 
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Discussion by N. K. Jain’ and Gopal 


would be expected and hardly general enough for use in design. In the writers 
~ opinions, analysis of data in nondimensional form would have been more useful. _ 
Test piles have been pushed into a medium dense dry sand deposit using 7 
a hydraulic jack. It would be of interest to know the initial relative density 
of sand and the change on account of driving. A surcharge | pressure of 100 
KN/ m? (14.5 psi) was applied on the sand pressure. However, the authors 
_ have not indicated any reasons for choosing this value and if any movement 
_ of pile was observed due to the application of surcharge. 
On the basis of observed pile top movement, the authors have identified x 
two distinct type of behavior with the dividing line between 0. 25 Qc/ 0-0. 30 ; 


attributed to reduction in | elastic ‘compression and redistribution of “residual — 4 
_ stresses. In the writers’ opinion, the surcharge pressure value is likely to influence 
the dividing line. Therefore, be relation 


Closure bys Sin-Fi Fatt Chan,* M. M. ASCE at and Thomas H H. er a oat 


di We thank Jain and Ranjan for their contribution. The pr purpose . of this esearch - 
_ work was to examine the general mechanics of long piles when subjected to 

_ slow repeated loadings. This work is not intended to be used directly for full- scale 

design but more to enable full- scale testing programs to be drawn up; consequently, — 
we believe that many of the variables and factors isolated in this work are 


a The piles were driven with the surcharge pressure in place and, . 
it is very difficult to know precisely what the relative density change was due | 
- to pile insertion. The value of 100 kN/m’ of surcharge pressure was chosen > 
4 “February, 1980, by Sin-Fatt Chan and Thomas H. Hanna (Proc. Paper 15222). _ 
- Research Scholar, Dept. of Civ. Engrg., Univ. of Roorkee, Roorkee 247672, India. | 
“Prof., and Head, Geotechnical Engrg., Dept. of Civ. ie Univ., of Ri Roorkee, Roorkee | 
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Me 
as a convenient value to raise ‘the stress level of the sand in the vicinity of b 
| ir pile. Its value was dictated in part by structural considerations of : movement — 4 
4 of the sand container and the patterns of behavior established are controlled be 
_ by many variables, including load level and load amplitude. Also the initial — 
residual stress system locked in the pile as a result of installation is extremely _ 
“a ‘important. It is quite clear to us that much more work has to be done before — 
full understanding of pile behavior under slow repeated loading is available. 
. The suggestion that the dividing line between distinct types of behavior is related _ 
_ to the value of the applied surcharge pressure is interesting. To. date this has r: 
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